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I.  INTRODUCTION 


General  Considerations 


During  the  thirty  years  since  B.  A.  Toms  first  published  his 
work  on  the  subject,  it  has  been  well  established  that  consid¬ 
erably  reduced  resistance  to  the  turbulent  shear  flows  of  liquids 
may  be  achieved  through  the  addition  of  small  quantities  of  high 
molecular  weight  polymers  into  the  fluid.  This  phenomenon  of 
drag  reduction  is  of  great  importance  and  has  far-reaching  impli¬ 
cations  in  the  reduction  of  power  required  for  the  pumping  of 
fluids  and  the  transport  of  bodies  through  liquids. 

Early  drag  reduction  efforts  were  primarily  concerned  with 
well-developed  pipe  flows,  the  area  in  which  the  first  observa¬ 
tion  was  made.  Until  recently,  researchers  have  directed  their 
efforts  toward  characterizing  the  manner  in  which  the  boundary 
layer  velocity  profiles  were  affected  in  an  attempt  to  quantify 
the  effect  of  polymer  additives.  Results  similar  to  those 
achieved  in  pipe  flows  have  been  obtained  in  investigations 
dealing  with  external  flows,  such  as  on  a  flat  plate.  The 
magnitude  of  the  turbulent  fluctuations  in  the  boundary  layer 
are  found  by  some  researchers  to  decrease  with  a  thickening  of  the 
laminar  sublayer  and  a  lowering  of  the  wall  shear  stress.  Other 
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researchers,  however,  have  found  no  evidence  of  this  thickened 
sublayer.  A  limit  to  the  amount  of  drag  reduction  achievable  with 
polymer  use  has  been  demonstrated  and  increases  in  the  shear 
stress  have  been  observed  when  excess  polymer  is  used.  Shear 
stress  reduction  for  external  flows  over  bodies  has  been  achieved 
by  the  ejection  of  polymer  solutions.  The  exact  mechanism  of 
shear  stress  reductions,  however,  still  remains  undefined.  Am¬ 
biguities  in  measurement  techniques  due  to  polymer  effects  on 
commonly  used  instrumentation  and  confusing  features  of  varied 
flow  facilities  have  contributed  to  making  past  investigations 
contradictory  and  difficult  to  interpret. 

It  is  evident  that  a  better  understanding  of  the  mechanism 
of  drag  reduction  In  developing  flows  with  polymer  ejection  is 
required  if  advantage  is  to  be  taken  of  the  phenomenon  for 
practical  application.  Most  cases  of  interest  in  external  flows 
would  require  that  a  carried  polymer  supply  be  used  for  ejection 
into  the  boundary  layer.  Presently,  very  little  is  known  about 
the  effect  of  ejection  techniques  on  polymer  distribution  in 
developing  boundary  layers  of  external  flows  and  no  information 
is  available  on  the  effects  of  polymer  on  boundary  layer  trans¬ 


ition  phenomena. 
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Scope  and  Objectives  of  Present  Studv 

The  object  of  the  research  described  herein  was  to:  (1) 
develop  a  flow  facility  suitable  for  the  study  of  polymer  flows; 
(2)  perform  velocity  measurements  using  non- flow-disturbing 
laser  Doppler  anemometry  techniques  with  real-time  data  process¬ 
ing;  (3)  develop  methods  for  the  prediction  of  local  skin  fric¬ 
tion  values  as  a  function  of  local  polymer  wall  concentration; 
(4)  determine  the  effects  of  ejection  techniques  on  the  distri¬ 
bution  of  polymer  in  the  boundary  layer;  and  (5)  determine  the 
diffusion  characteristics  of  a  polymer  when  ejecting  it  into  a 
laminar  boundary  layer  and  also  its  effect  on  transition. 

Approach 


A  flow  facility  was  constructed  under  the  constraints  of 
limited  funds  and  practicality  of  construction.  The  basic 
tunnel  structure  was  designed  with  a  width-to-height  aspect  ratio 
of  5:1  to  minimize  the  effects  of  side  wall  boundary  layer 
growth  and  corner  flows  in  the  rectangularly-shaped  channel. 

Flow  delivery  was  by  constant  head  gravity  feed  to  minimize 
vibration  effects.  The  polymer  delivery  system  provided  minimum 
mechanical  degradation  of  polymer  solutions  while  allowing 
variations  in  polymer  concentrations  and  ejection  velocities. 
Transition  from. laminar  to  turbulent  flow  was  achieved  at  midplate 


position  allowing  ejection  of  polymer  solutions  into  a  developing 
laminar  boundary  at  the  plate  leading  edge  and  transitioning  to 
turbulent  flow  along  the  plate.  Boundary  layer  velocity  profiles 
through  transition  to  turbulent  flow  were  taken  with  laser  Doppler 
anemometry  techniques,  non-disruptive  of  the  flow  nor  affected  by 
polymer  characteristics. 

Using  boundary  layer  sampling  techniques,  polymer  diffusion 
characteristics  were  measured  at  five  axial  wall  locations  along 
the  plate,  as  well  as  at  eight  positions  normal  to  the  plate  at 
each  wall  station.  The  ejected  polymer  concentrations  were  con¬ 
taminated  with  a  fluorescent  dye,  which  was  later  used  to  analyze 
the  concentrations  of  the  boundary  layer  samples.  Thus  it  was 
possible  to  measure  the  change  in  wall  concentration  and  concen¬ 
tration  profile  in  the  boundary  layer  during  laminar,  transitional, 
and  turbulent  flow. 

Analytical  methods  were  developed  to  predict  skin  friction 
coefficients  as  a  function  of  local  polymer  wall  concentrations 
for  comparison  with  experimentally  determined  values. 


II.  LITERATURE  REVIEW 


Historical  Background 


B.A.  Toms  (19^9)  published  the  first  data  on  the  friction 
reducing  effects  of  polymers  and  showed  chat,  for  conditions  of 
turbulent  flow  with  constant  pressure  gradient,  the  average 
velocity  of  monochlorobenzene,  flowing  in  a  pipe,  increased 
substantially  with  the  addition  of  dilute  quantities  of  polymethyl 
methacrylate.  Friction  reduction  of  up  to  50%  was  achieved 
(as  compared  with  the  pure  solvent)  by  the  addition  of  .25% 
by  weight  of  polymer  to  the  solvent. 

Oldroyd  (1949)  attempted  to  explain  this  phenomenon 
(presently  referred  to  by  many  as  the  "Toms  effect")  as  a  wall 
effect  due  to  exclusion  of  large  polymer  molecules  from  a  region 
near  the  wall  due  to  their  size.  As  a  result  of  studies  on  the 
flow  of  gasoline  thickened  with  aluminum  soaps  at  Eagewood  Arsenal 
during  World  War  II  (but  not  published  at  that  time),  Mysel  (1949) 
applied  for  and  received  a  patent  in  1949.  He  observed  that  in 
tubulent  flow  the  pressure  loss  per  unit  length  of  pipe  was  much 
lower  for  thickened  gasoline  than  that  of  the  pure  solvent. 

The  importance  of  drag  reduction  seemed  to  have  been  over¬ 
looked  until  the  late  1950s  when  Shaver  and  Merril  (1959)  and 
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Dodge  and  Metzner  (1959)  published  results  indicating  the  low 
friction  factors  of  non-Newtonian  solutions,  such  as  sodium 
carboxymethycellulose ,  in  water.  Similar  observations  of  reduced 
friction  were  made  in  the  oil  well  industry  with  the  use  of  the 
plant  derivative.  Guar  Gum.  Solutions  of  the  product  were  used 
to  suspend  sand  in  the  high-pressure,  sand-water  mixtures  utilized 
in  oil  well  fracturing  operations. 

These  experiences  led  to  U.S.  Navy  studies  in  the  early  1960's 
on  possible  military  applications  of  friction-reducing  effects. 
Among  Navy  reseachers,  Hoyt  (1963)  made  the  first  significant 
contribution  followed  by  Fabula  (1964)  who  discovered  the 
spectacular  friction  reducing  ability  of  polyethylene  oxide,  the 
most  effective  friction-reducing  material  known. 

Crawford  (1962)  and  Savins  (1961)  led  the  early  efforts  of 
industrial  research  centers.  The  years  that  have  followed  these 
early  efforts  have  produced  an  ever-increasing  number  of  inter¬ 
national  papers,  reports,  and  conferences  on  the  subject. 

Excellent  summaries  of  research  efforts  and  results 
available  in  the  field  of  drag  reduction  have  been  presented  by 
Lumley  (1969)  and  Hoyt  (1972).  Palyvos  (1974),  however,  prepared 
the  most  extensive  and  detailed  review  of  research  efforts  to  date. 
This  review  reveals  that  after  years  of  active  research  and  the 
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combined  efforts  of  hydrodynamicists,  polymer  chemists  and 
rheologists  in  many  countries,  there  is  yet  no  adequate 
theoretical  explanation  for  the  mechanism  of  drag  reduction,  nor 
is  there  a  satisfactory  correlation  of  the  viscoelastic  properties 
of  dilute  polymer  solutions  with  their  friction  reduction  effect. 

Literature  pertaining  to  the  current  study  will  be  reviewed 
under  four  separate  categories: 

1.  Turbulent  boundary  layer  theory 

2.  Homogeneous  polymer  flow 

3.  Polymer  ejection  studies 

4.  Experimental  facilities  • 

Turbulent  Boundary  Layer  Theory 

The  concept  of  transition  from  laminar  tc  turbulent  flow  was 
first  demonstrated  in  1880  by  Osborne  Reynolds  in  his  classic 
pipe  flow  experiment.  The  next  major  step  did  not  occur  until 
1904  when,  by  demonstrating  the  existence  of  a  thin  "boundary 
layer"  in  fluid  flows,  L.  Prandtl  allowed  the  reconciliation  of 
viscous  flow  with  classical  frictionless  hydrodynamic  relations. 

In  1914,  by  showing  that  boundary  layers  could  be  either  laminar 
or  turbulent,  Prandtl  freed  early  investigations  from  the  limi¬ 
tation  of  considering  only  laminar  boundary  layers.  The 


3 


introduction  of  the  "Prandtl  mixing  length  theory"  in  1925  was  a 
major  contribution  toward  the  understanding  of  the  development  of 
the  boundary  layer  velocity  profile.  Since  that  time  there  have 
been  an  ever-increasing  number  of  contributions  to  the  under¬ 
standing  of  the  flow  of  fluids. 

Many  attempts  have  been  made  to  predict  the  conditions  that 
govern  the  transition  of  laminar  to  turbulent  flow  as  well  as  to 
predict  velocity  profiles  and  wall  shear.  H.  Schlichting's  classic 
text,  Boundary  Layer  Theory  (1968),  remains  the  most  comprehensive 
treatment  of  the  subject  to  date.  In  this  work,  laminar  and 
turbulent  boundary  layers  in  both  compressible  and  incompressible 
flows  over  a  wide  range  of  application  are  addressed.  However, 

t 

analyses  of  the  turbulent  boundary  layer  problem  still  rely  heavily 
on  a  combination  of  dimensional  analysis,  empirical  data  and  flow 
visualization  techniques.  Also,  during  1968,  an  attempt  was  made 
to  introduce  order  into  turbulent  boundary  layer  research.  A 
conference  of  the  world’s  leading  boundary  layer  researchers  was 
called  at  Stanford  University.  Kline  (1968)  edited  the  proceed¬ 
ings  of  that  conference,  in  which  some  29  methods  and  approaches 
to  the  analysis  of  the  turbulent  boundary  layer  were  presented 
and  discussed. 

Viscous  Fluid  Flow,  a  text  by  F.M.  White  (1974),  discusses 
many  of  the  more  significant  methods  of  analysis.  Two  major 


approaches  evolved:  (1)  integral  methods  averaged  across  the 
boundary  layer,  and  (2)  finite  difference  techniques  which  attempt 
to  solve  the  full  partial  differential  equations  of  the  boundary 
layer . 

In  the  last  analysis,  it  remains  that  the  empirical  dimen¬ 
sional  analysis  approach  is  still  the  most  useful  for  engineering 
purposes.  The  scientific  world  is  indebted  to  the  physical  insight 
of  Prandtl  and  Von  Karman  for  the  formulation  of  the  turbulent 
velocity  profile  concept.  Consideration  of  regions  of  influence 
where  the  relative  importance  of  viscous  shear  or  turbulent  shear 
dominate  led  to  the  formation  of  the  three-layer  concept  for  the 
velocity  distribution  in  a  turbulent  boundary  layer.  The  layers 
are: 

Inner  Layer:  Viscous  shear  dominates 

Outer  Layer:  Turbulent  shear  dominates 

Overlap  Layer:  Viscous  and  turbulent  shear  equally 
important . 

The  mean  velocity  distribution  in  a  two-dimensional  turbulent 
boundary  layer,  u(y) ,  depends  upon  four  local  parameters:  the 
local  wall  shear  stress,  tw;  the  fluid  density,  p;  the  fluid  vis¬ 
cosity,  p;  and  the  boundary  layer  thickness,  5.  For  the  inner 
law,  Prandtl  (1933)  deduced  that  the  mean  velocity  depended  on  the 
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wall  shear  stress,  the  fluid  properties,  and  the  distance  y  from 
the  wall  which  led  to  the  following  functional  relation  for  the 
inner  law: 


u  -  f(V  p,  u,  y) 


(1) 


For  the  outer  layer,  von  Karman  deduced  that  the  wall  acts  only  as 
a  source  of  retardation  reducing  the  local  velocity,  u,  below  the 
freestream  velocity,  U£,  in  a  manner  independent  of  the  viscosity, 
U.  The  outer  law  or  velocity  defect  relation,  as  it  is  sometimes 
called,  becomes 


U  -  u  -  f Ct  ,  p,  y,  w) . 
e  w 


C2) 


Coles  (1954)  performed  a  dimensional  analysis  on  these 
* 

relationships,  where  V  is  a  characteristic  velocity  called  the 
wall  shear  velocity,  and  is  defined  as 


(3) 


Sear  the  wall,  Coles  (1954)  found  that  the  velocity  profile 
is  unaffected  by  the  boundary  layer  thickness,  therefore 


+  u 

Inner  Law  U  ■ 


f.Z* 


) 


(4) 
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Far  from  the  wall,  viscosity  does  cot  affect  the  velocity 
distribution.  (The  viscosity  effects  are  concentrated  near  the 
wall.) 

(„,) 

Outer  Law  *■  \  '  ■  g(y/5)  (5) 

V 

In  the  overlap  region  between  the  inner  and  outer  layers, 
the  commonly  known  relation  for  the  "law  of  the  wall"  (shown  in 
figure  1)  may  be  obtained  by  equating  equations  (4)  and  (5),  such 
that 


(6) 


The  constants  X  and  3  have  been  determined  by  the  data  of 
Nikuradse  (1930)  to  be  X  *  .40  and  3  =*  5.5. 


Coles  (1954)  presents  a  correlation  of  the  dimensionless 
* 

velocity  profile,  u/V  ,  with  dimensionless  distance  from  the 
wall,  yvjv,  up  to  a  value  of  about  300.  Neglecting  separating 
flows,  the  data  collapse  into  regions  governed  by  the  inner  and 
logarithmic  outer  laws  with  the  following  limits: 

Inner  Law: 

* 

V 


W 

»  -  ■ 

V 


o  < 


vV 


<  10 


(7) 
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Logarithmic  Law: 

*  * 

-*  -  2.5  ln^  +5.5,  35  <  ^  <  300.  (8) 

7  v  v 

* 

+  vV 

For  values  of  y  =  >  300,  the  viscosity  is  of  negligible 

importance  in  the  determination  of  the  shape  of  the  boundary 
layer.  The  data  will  then  correlate  with  the  velocity  defect  law 
(shown  in  figure  2)  in  such  a  way  that 
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In  a  subsequent  paper,  Coles  (1955)  postulates  that  a  two- 
dimensional  boundary  layer  of  an  incompressible  fluid  may  be 
represented  by  a  linear  combination  of  two  functions: 


The  function  f(yV*/v)  is  given  by  equations  (7)  and  (8).  Coles 
refers  to  the  function  W(y/5)  as  the  "law  of  the  wake."  The 
function,  H,  is  related  to  Clauser's  (1954,  1956)  equilibrium 
parameter,  8,  designed  as  follows: 


(ID 


After  additional  study,  Clauser  concluded  that  the  typically 


fuzzy  thickness,  5,  should  be  replaced  by  the  displacement 


thickness,  <5^,  so  that  the  accepted  parameter  is  now 


3  ” 


(12) 


The  "wake  parameter,"  H,  has  a  reasonable  fit  to  data  when  given 
by  the  relationship 

n  a  0. 8(3  +  0.5)°*  75 .  (13) 


Coles  (1956)  proposed  the  following  curve  fit  for  the  wake 
function,  W,  as 

W  f  =  2sin2  f.f.  (14) 

Then  in  Cole's  notation  for  the  overlap  and  other  layers,  we  nave 
the  following 

u+*  ■^•ln(y+)  +  B  +  £  W  (^)  .  (15) 

It  may  then  be  seen  that,  near  the  wall,  the  velocity 
profile  is  dominated  by  friction  forces  and  the  function 
f(y  •  uT/v)  controls.  Away  from  the  wall,  the  profile  is  dominated 
by  inertia  forces,  with  the  function  w(y/5)  controlling. 


Kline,  Runstadler,  and  Reynolds  (1963)  visually  demonstrated 
the  existence  of  three  regions  of  turbulent  boundary  layers,  using 
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a  combination  of  pulsed  hydrogen  bubble  techniques,  hot  film  probes, 
and  dye  streak  photographs.  Shown  in  figure  3,  the  flow  regions 
of  the  turbulent  boundary  layer  model  are: 

1.  A  laminar  flow  region  having  a  near  regular,  three-dimen¬ 
sional  flow  pattern  for  a  range  of  y'  from  0  <  y(u^./\>)  <  10. 

2.  A  fully  turbulent  region  of  randomly  fluctuating  eddies. 
This  region,  sometimes  called  the  inner  turbulent  region,  exists 
over  a  range  of  y+  from  10  <  y(u^/v)  <  370. 

3.  An  outer  turbulent  region,  or  wake  region,  extending 
from  y+  of  approximately  370  to  the  outer  edge  of  the  boundary 
layer.  In  this  region,  large  scale  turbulent  eddies  are 
found . 


Brady  (1973)  summarized  the  work  of  Kline  et  al.  (1963)  as 
follows : 


Kline  et  al.  found  the  'laminar  sublayer'  to  be  made 
up  of  a  regular  structure  of  low  and  high  velocity  longi¬ 
tudinal  streaks  which  meander  transvers^  and  normal  to  the 
wall.  A  dimensionless  streak  spacing  X  =  \v*/v  *  76.5 
was  found  for  zero  pressure  gradients.  These  streaks 
either  break  up  or  randomly  burst  from  the  sublayer  into 
the  fully  turbulent  region. 

The  fully  turbulent  region  is  one  of  intense  mixing 
and  high  dissipation  of  energy.  Protruding  from  this 
region  are  intermittent  large  eddies  —  visualized  as 
peninsulas  of  turbulence. 


figure  3.  Pictorial  (sketch  of  boundary  layer  flow  (from  Kline  (1963)) 


18 


The  peninsulas  of  turbulence  which  extend  into  the 
third  boundary  layer  zone  —  the  outer  turbulent  region 
—  gives  it  a  characteristic  not  unlike  the  wake  behind 
a  cylinder.  It  has  rather  large  lumps  of  turbulence 
at  intermittent  spacings. 

It  is  clear  that  the  driving  force  for  the  entire 
turbulent  boundary  layer  is  the  generation  of  velocity 
streaks  in  the  sublayer  and  their  subsequent  bursting 
outwards.  Kline  presents  the  results  of  many  other 
investigators,  as  well  as  his  own,  in  support  of  this 
hypothesis . 


Homogeneous  Polymer  Flow 

Since  the  early  work  of  Toms  (1949),  Oldroyd  (1949),  and 
Mysel  (1949) ,  many  high  molecular  weight  polymers  have  been  shown 
to  be  effective  drag  reducers.  Among  the  most  prominent  are: 
polysaccharides  (Guar),  polyethylene  oxide,  polyacrylamides,  and 
sodium  carboxymethyl  cellulose.  For  the  most  part,  investigators 
have  considered  mainly  pipe  flows,  therefore,  this  effort  will 
be  considered  first. 

Pipe  Flow 


Table  1  from  Hoyt  (1972)  shows  the  effectiveness  of  small 
concentrations  of  polymer  on  drag  reduction.  Listed  are  those 
concentrations  required  to  achieve  67%  drag  reduction  in  pipe 


flow  at  Re  =  14  X  10^. 


TABLE  1 


CONCENTRATIONS  (WPPM)  OF  MATERIAL  REQUIRED  TO  ACHIEVE  67% 
DRAG  REDUCTION  IN  PIPE  FLOW  OF  Re  “  14  x  103 
(FROM  HOYT  (1972)) 


Polymer 
Gum  Karaya 
Guar 

Polyacrylamide,  Polyhall-250 
Polyox  WSR-301 


Concentration  (WPPM) 
850 
400 
20 
10 


Hoyt  and  Fabula  (1964)  and  Virk  (1971)  present  data  that  show  a 
maximum  drag  reduction  asymptote.  For  a  smooth  pipe,  this 
asymptote  corresponds  to  80%  of  the  friction  reduction  that  would 
be  attained  if  completely  laminar  flow  were  sustained  at  a  given 
Reynolds  number. 


Meyer  (1966)  and  Elata  et  al.  (1966)  have  shown  that  drag 

reduction  in  pipes  is  due  to  a  thickening  of  the  laminar  sublayer. 

It  was  shown  that  the  constant  B  in  the  law  of  the  wall  equation 

(6)  remained  constant  and  equal  to  the  Newtonian  value  until  a 

* 

critical  threshold  value  of  the  shear  velocity,  V  ,  was  reached. 

* 

Thereafter  the  value  of  3  increased  logarithmically  with  V  : 


20 


B  -  5.5  4-  Z  In 

B  -  5.5  +  /IB, 

where  AB  *  2 

Using  data  from  several  investigators, 
found  that: 


(16) 

(17) 

(18) 


F.M.  White  (1963) 


where  a 


2  -  aCY, 
w 

2.3,  and  y  *  0.5, 


(19) 

(20) 


with  a  maximum  value  of  2  of  approximately  11.  The  data  further 

* 

indicated  that  the  critical  shear  velocity,  V  for  the  onset  of 
drag  reduction  was  .08  ft/sec.  Virk  (1966)  deduced  from  exper¬ 
imental  data  that  the  onset  shear  stress  is  inversely  proportional 
to  the  polymer  molecular  radius  of  gyration.  The  critical  wall 
shear  stress,  t®,  which  must  be  exceeded  for  drag  reduction 
to  occur  is  given  by 

c  ■ .  . . i  .  6  \2  cons t an t  / \ 

t  -  p(0.624  X  10  y/R_)  -  - o -  (21) 

w  G 

XG 

where  R„  is  the  rms  radius  of  gyration  of  the  molecule  obtained 

<3 


from  light  scattering  data. 
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Fabula  et  al.  (1969)  developed  a  criteria  for  the  critical 
wall  shear  stress  at  the  onset  of  drag  reduction  based  on 
molecular  deformation  times  and  a  molecular  characteristic  relaxa¬ 
tion  time.  It  is  thought  that  individual  molecules  are  too 
small  by  several  orders  of  magnitude  to  interfere  with  the 
turbulence  structure,  indicated  by  the  small  value  of  the  ratio 
of  the  polymer  molecule  scale  to  the  scale  of  turbulent  eddies 
at  onset. 

To  circumvent  this  problem  of  length  scales,  Fabula 
postulated  an  interaction  between  the  time  scales  of  the  periodic 
molecular  deformation  in  the  viscous  sublayer,  given  by  y/ZR, 
where  •?  is  the  shear  rate,  and  a  molecular  characteristic 
relaxation  time,  t^.  The  molecular  charcteristic  relaxation 
time,  r  ,  may  be  determined  by  the  Zimm  or  Rouse  theories  which 
relate  a  characteristic  relaxation  time  of  the  solution,  t  ,  to 
the  solvent  viscosity,  us;  the  solution  viscosity,  u;  the  polymer 
molecular  weight,  M;  temperature,  T;  and  the  concentration  of 
polymer,  C.  It  follows  then  that 

_  (22) 

Tl  "  CRT 

where  a  is  a  constant  having  a  value  between  .4  and  .6. 
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The  resulting  criterion  for  the  critical  wall  sheer  stress  is 

given  by 


2n(r/2V  rz  -  i 
or 

c  m  u_  M  uCRT 
w  ”  t z"  a(u-ua)M 


(23) 


Equation  (23)  may  be  written  using  a  relation  between  intrinsic 
viscosity,  molecular  weight,  and  molecular  dimensions  yielding: 


c 

T 

W 


constant 


1.7  X  10 


(24) 


Equation  (24)  of  Fab  Vila  is  similar  to  Virk's  equation  (21). 
The  former  provides  a  better  estimate  of  magnitude  of  onset  shear 
stress  whereas  the  latter  is  considered  more  accurate  when  the 
constant  is  determined  from  test  data. 


In  their  work,  Kowalski  and  3rundrett  (1974)  support 

/ 

Fabula's  postulate  that  the  effectiveness  of  very  dilute  solution 
is  due  to  entanglements  or  "blobs"  of  macromolecules  rather  than 
individual  molecules .  A  formula  was  developed  relating  the  size 
of  the  entangled  molecules  with  the  size  of  a  dissipative  eddy 
and  are  tested  to  predict  the  onset  of  drag  reduction  in  pipe 
flows  of  homogeneous  polymer  solutions. 
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Darby  (1972)  reviews  drag  reduction  theories  comparing 
molecular  hypotheses,  continuum  approaches,  and  conventional 
boundary  layers  length  scales.  Molecular  time  scale  hypotheses 
are  represented  by  the  work  of  Virk,  Fabula,  and  Kowalski.  The 
continuum  approach  of  Seyer  and  Metzner  concludes  that  the 
presence  of  elastic  properties  in  the  dilute  solution  are 
sufficient  criteria  for  drag  reduction.  Elastic  properties  are 
represented  by  a  single  parameter,  the  relaxation  time,  X, 
incorporated  into  a  dimensionless  group,  the  Deborah  number: 


N 


De 


x/tx, 


(25) 


where  t^  is  a  characteristic  time  of  the  system.  Seyer  and 
Metzner  (1969)  chose  this  number  to  be  the  reciprocal  of  the  dissi¬ 
pative  turbulence  frequency, 

a 


a) 


d 


In  3/4 

D  *  Re 


(26) 


Seyer  and  Metzner  further  offer  the  following  phenomena  as  a 
possible  explanation  of  drag  reduction: 


1.  Particulate  effects  due  to  large  molecules  or  agglomerates 


which  may  promote  stability  of  laminar  flow  or  dampen  turbulence 
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2.  The  solid-like  resistance  of  the  fluid  to  sudden  defor¬ 
mation  or  stretching  may  suppress  the  dissipative  turbulent 
frequencies  on  the  generation  of  turbulence 

3.  The  laminar  flow  regime  may  be  stabilized  so  that 
transition  is  delayed  to  higher  Reynolds  numbers  • 

Darby  (1972)  points  out  that,  although  Seyer  and  Metzner 
favor  the  second  concept,  no  distinction  is  made  between  cause 
and  effect.  Dozens  of  different  forms  of  the  Deborah  number 
are  presented,  illustrating  the  problems  encountered  with  the 
present  theoretical  state  of  the  art. 

Boundary  layer  modification  in  terms  of  the  Prandtl  mixing 
length  theory  and  law  of  wall  relation  indicates  a  concept  of 
the  flow  field  divided  into  three  regions: 

1.  A  laminar  sublayer  adjacent  to  the  wall  where  turbulent 
fluctuations  are  absent  and  momentum  transport  is  by  viscous 
shear  only 


2.  A  transition  or  "buffer"  zone  adjacent  to  the  laminar 
sublayer  where  both  viscous  shear  and  turbulent  inertial  fluc¬ 
tuations  contribute  to  the  momentum  transport 


3.  A  turbulent  boundary  layer  or  "core"  outside  the 
buffer  zone  where  transport  is  dominated  by  turbulent  inertial 
fluctuations. 

Good  results  for  friction  factor  correlations  have  been 
achieved  by  Elata,  Lehner,  and  Kahanovity  (1966)  for  Guar  Gum 
solutions  and  Meyer  (1966)  and  Wells  (1965)  for  Polyox. 

Many  authors  have  described  drag  reduction  as  a  "negative 
roughness  effect"  since  polymers  appear  to  thicken  the  sublayer 

•f 

while  maintaining  the  same  slope  of  the  U  versus  In  Y  curve 
in  the  overlap  region.  Nadolink  (1968)  demonstrated  the  existence 
of  the  thickened  sublayer  directly,  using  a  high-speed  motion 
picture  camera  and  a  microscope. 

White  and  McEligot  (1970),  in  their  work  on  transition 
delay  from  laminar  to  turbulent  flow,  found  a  dependence  on 
where  the  critical  onset  shear  stress  is  reached.  If  the  onset 
shear  stress  occurs  in  the  laminar  flow  region,  a  delay  in 
transition  to  turbulent  flow  can  occur. 

Polymer  Injection  Studies 

Practical  application  of  the  phenomena  of  viscous  drag 
reduction  caused  by  dilute  polymer  solutions  Includes  increasing 


speed  or  power  reduction  requirements  for  naval  vessels  and 
increasing  the  flow  through  conduits.  Realization  of  these 
advantages  generally  requires  injection  of  polymer  solutions 
into  a  developing  boundary  layer.  Observation  that  drag  reduc¬ 
tion  is  related  to  changes  in  the  sublayers  implies  that  injecting 
polymer  near  the  wall  affects  considerable  drag  reduction  with 
relatively  small  amounts  of  polymer.  Literature  pertinent  to  this 
area  of  the  field  will  now  be  reviewed. 

Pipe  Flows 


Wells  (1968)  suggested  uniform  injection  through  a  porous 
wall  since  it  raises  the  additive  concentration  to  the  drag 
reducing  level  in  the  wall  region  only.  Using  a  Reynolds-Prandtl 
analogy  to  analyze  the  diffusion  process,  he  calculated  that 
distributed  injection  would  require  40  to  140  times  less  additive 
than  slot  injection  to  maintain  equivalent  drag  reduction.  The 
point  was  made  that  a  continuously  ablating  additive  coating 
would  yield  the  maximum  performance  advantage  if  it  could  be  made 
to  ablate  at  the  optimum  rate. 

Wells  and  Spangler  (1967)  performed  a  series  of  experiments 


in  order  to  determine  whether  the  presence  of  the  additive  only 
in  the  wall  region  could  produce  significant  local  shear  stress 
reduction.  The  experiments  utilized  centerline  injection  and 


circumferential  slot  injection  of  dilute  polymer  into  fully- 
developed  turbulent  flow  of  a  Newtonian  fluid.  It  was  found  that 
the  local  pressure  loss  was  reduced  by  an  amount  comparable  to 
the  flow  of  a  uniform  concentration,  when  the  fluid  was  injected 
at  the  pipe  centerline.  However,  no  reduction  in  local  pressure 
loss  occurred  until  the  fluid  diffused  into  the  wall  region. 

Maus  and  Wilhelm  (1970)  conducted  polymer  injection  tests 
with  a  fully-developed  flow  in  a  1.625-inch  diameter  smooth  acrylic 
pipe.  Five  circumferential  injection  slots  located  6  inches 
apart  were  used.  Each  slot  was  0.050-inch  wide  and  inclined  30° 
to  the  pipe  centerline.  The  effects  of  Reynolds  number,  injection 
rate,  number  of  injection  points,  and  concentration  of  injected 
solution  were  studied.  It  was  found  that  maximum  drag  reduction 
occurred  when  polymer  was  injected  through  the  furthest  upstream 
slot  rather  than  being  distributed  over  the  test  length. 

Walters  and  Wells  (1971)  conducted  tests  using  uniform 
ejection  of  polymer  solution  through  a  stainless  steel  porous  pipe 
section  into  fully-developed  turbulent  pipe  flow.  Concentration 
profiles,  velocity  profiles,  and  wall  shear  stress  data  were 
obtained.  Fluorometric  techniques  were  used  to  obtain  the  concen¬ 
tration  profiles.  Uranine  B  was  selected  as  the  fluorescent  dye 
over  several  others,  as  the  laboratory  apparatus  could  be  easily 
decontaminated.  The  dye  was  mixed  with  the  injection  fluid  and 
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Chen  captured  tracer  sample  concentrations  were  measured  with  a 
fluorometer.  Accuracies  of  the  measurement  instrumentation  of 
two  parts  per  billion  were  reported.  Uniform  polymer  injection  at 
the  wall  appeared  to  somewhat  inhibit  turbulent  diffusion  away 
from  the  wall  resulting  in  high  wall  concentrations.  For  certain 
conditions  of  high  polymer  mass  flux,  a  wall  friction  increase  was 
noted,  possibly  due  to  the  higher  viscosity  in  this  region.  As 
compared  to  water  injection,  a  one  to  two  order  of  magnitude 
reduction  in  total  diffusivity  in  the  ejection  region  was  noticed. 
Downstream  of  the  ejection  section,  an  order  of  magnitude  reduction 
of  total  diffusivity  was  noted  along  with  a  significant  reduction 
in  wall  friction. 

Tullis  and  Ramu  (1973)  studied  the  characteristics  of  mean 
turbulent  flow  in  the  entrance  region  of  a  rough  pipe  for  water 
and  for  polymer  injection  into  a  boundary  layer.  Polvox  WSR-301 
was  injected  through  a  perforated  wall  pipe  section  of  a  12-inch 
diameter,  200-foot-long  steel  pipe  used  for  the  study.  Drag 
reduction  of  up  to  80%  in  the  fully-developed  region  and  90%  in 
the  inlet  region  were  recorded.  Comparison  of  water  and  dilute 
polymer  ejection  showed  the  polymer  concentration  profiles 
developed  slower  than  that  of  the  water  alone,  indicating  lower 
diffusion  of  the  polymer  solution.  The  inlet  length  needed  for 
flow  to  fully  develop  was  found  to  be  greater  for  polymer 
injected  flows  that  for  the  case  of  no  injection.  In  the  fully 
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developed  region,  the  maximum  drag  reduction  measured  was  80%  at 
■  25  WPPM.  Drag  reduction  in  the  inlet  region  appeared  to  be 
independent  of  injection  concentration.  Fluorescent  techniques 
using  rhodamine  WT  dye  were  used  to  measure  the  small  concentra¬ 
tion  of  polymer.  Mean  flow  velocities  between  5  and  50  ft/sec. 

(R^  -  3  X  10^  to  3  X  10^)  were  used.  Injection  concentration 
varied  from  100  to  2400  WPPM  at  rates  of  30  to  220  gal/min 
at  a  location  3.5  pipe  diameters  from  the  pipe  inlet. 

Sellin  (1974)  reported  drag  reduction  of  40-50%  in  large 
scale  tests  at  R^  of  1.5  X  10^.  A  Polyox  WSR-301  solution  was 
added  to  water  flow  in  a  203-mm  diameter  pipeline,  4190  meters 
in  length  at  the  Bristol  Sewage  Treatment  Work  at  Avonmount, 
England.  Polymer  powder,  corresponding  to  a  final  concentration 
of  40-60  WPPM,  was  continuously  mixed  in  a  votex  chamber  and  the 
resulting  slurry  pumped  into  the  pipeline  using  a  gear  pump. 

Open  Channel  Flows 

Latto  and  Shen  (1970)  performed  an  experiment  of  slot 
injection  over  a  flat  plate  suspended  in  a  closed  loop  flume. 
Concentration  of  Polyhall-295  from  200  to  600  WPPM  were  ejected 
at  a  20°  angle.  Using  hot  film  anemometry  it  was  found  that 
momentum  diffusion  was  less  than  that  for  pure  water.  The  flow 
rate  velocity,  and  angle  of  ejection  of  the  polymer  solution  were 
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also  found  to  be  important.  Large  percentage  drag  reductions 
were  reported  by  injecting  high  concentration  aqueous  polymer 
solutions  into  turbulent  boundary  layers  at  low  flow  rates.  It 
was  found  that  tangential  injection  at  the  lowest  possible 
velocity  was  desirable. 

Wetzel  and  Ripkin  (1970)  experientally  studied  the  injection 
of  Polyox  WSR-301  into  a  developing  boundary  layer  in  a  9-foot , 
wide  open  channel.  Pitot  tubes  were  located  at  positions  16,  28, 
and  60  feet  downstream  of  the  injection  slot.  Polymer  injection 
resulted  in  much  fuller  velocity  profiles  than  did  water  injection. 
Several  methods  of  polymer  concentration  determination  were  in¬ 
vestigated  including:  polarography,  the  turbidimetric  method,  and 
the  fluorimetric  method.  The  fluorimetric  method  was  found  to 
give  rapid,  accurate  results.  Concentration  profiles  for  water, 
1000,  2000,  and  3000  WPPM,  were  found  to  be  in  good  agreement  with 
the  curves  developed  by  Morkorin  (1963) . 

A  maximum  drag  reduction  of  25%  was  experienced  over  the  40- 
foot  boundary  layer  length,  which  resulted  in  a  downstream  wall 
concentration  of  30  WPPM.  At  a  distance  of  16  feet  from  the  slot, 
greater  drag  reduction  was  attained  for  the  lower  concentrations 
injected  than  for  the  larger  quantities.  At  distances  further 
downstream  the  reverse  was  true,  with  better  drag  reduction 
attained  with  the  larger  injection  quantities. 
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The  behavior  was  attributed  to  more  complete  mixing. 

Color-dyed  polymer  injection  tests  revealed  that  shortly  after 
injection,  the  flow  organized  into  a  definite  pattern  of  large 
wavering  parallel  streaks.  This  streaking  was  reported  to  be  a 
secondary,  three-dimensional  vortex  motion  superimposed  on  the 
two-dimensional  flow. 

Fabula  and  Burns  (1970)  invoked  the  negative  roughness 

analogy  so  that  the  outer  layer  mean  velocity  similarity  is 

unaffected  by  friction  reduction.  The  similarity  law  of  mixing 

with  polymeric  friction  is  predicted  to  be  the  same  as  without 

polymeric  friction  reduction.  The  velocity  defect  law  of  Coles 

(1956)  as  well  as  his  value  for  the  law  of  the  wake  parameter  II 

0 

of  0.55  were  adopted.  A  relationship  for  calculating  the  local 
additive  concentration  at  downstream  stations  along  the  wall  was 
developed  and  compared  to  preliminary  data  on  mixing  in  an  open 
channel  boundary  layer  injected  with  Poly ox  WSR-301. 

Latto  and  Middleton  (1970)  reported  on  extensive  velocity 
profile  data  taken  with  hot  film  probes  for  flow  of  homogeneous 
solution  of  Polyacrylamine  MRL-402  over  a  flat  plate  suspended  in 
a  recirculating  flume.  Direct  drag  measurements  were  made  over  a 
period  of  time  for  concentrations  of  0,  25,  50,  and  75  WPPM  which 
indicated  no  appreciable  degradation  of  the  polymer.  The  turbulent 
boundary  layer  profiles  appeared  "fuller"  than  corresponding 
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profiles  for  pure  water.  A  friction  velocity  of  *  0.074  ft/sec 
was  found. 

Fruman  and  Tulin  (1974)  performed  a  study  of  diffusion  of  a 
thin  tangential  jet  of  Polyox  WSR-301  solution  injected  into  the 
turbulent  boundary  layer  of  a  flat  plate  vertically  suspended  in 
a  free— surf ace ,  high-speed  water  channel.  Freestream  Reynolds 
numbers  of  the  order  of  3.6  X  10^  were  achieved.  Drag  measurement 
by  reluctance  force  gages  were  taken.  Wall  concentrations  of 
polymer  were  taken  using  a  light-intensity,  dyed-additive  method. 

The  concentration  distribution  along  the  wall  was  found  to  be 
represented  by  two  regions. 

Within  the  first  region,  the  wall  concentration  was  prac¬ 
tically  constant  and  equal  to  the  injected  polymer  concentration. 

In  the  second  region,  the  concentration  varied  inversely  with  the 
distance  from  the  injection  slit.  The  length  of  the  first  zone 
was  found  to  be  15  to  20  times  that  of  water  injection.  This 
extended  initial  zone  appeared  to  be  directly  related  to  the 
thickening  of  the  viscous  sublayer,  the  reduction  of  shear  stress, 
and  the  decrease  of  molecular  diffusivity. 

In  the  region  where  turbulent  diffusivity  was  predominant, 
both  polymer  solution  and  solvent  flows  displayed  similar  behaviors. 
The  distribution  of  the  wall  temperature  over  a  flat  plate  given 
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ov  Seban  (1960)  provided  a  hear  transfer  analogy  for  correlation 
of  data.  According  to  Seban: 


T 
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Whereas ,  Tulin  indicated  o  ■  a  ,  and  therefore 


.  17.01  '‘''(X/S)-0-711. 

Ci  u 


(28) 


Closed  Channel  Flows 


Wu  (1969)  presented  experimental  data  obtained  by  injection 
of  various  polymer  concentrations  along  a  flat  plate,  and  trends 
similar  to  Porah  at  al.  (1971)  were  noted  (lower  initial  concen¬ 
trations  caused  higher  drag  reductions) .  The  results  of  Wu  indi¬ 
cate  poor  mixing  between  the  injected  fluid  and  its  surroundings. 
Photographs  of  submerged  jets  confirm  that  additives  suppress 
turbulent  diffusion.  The  flat  plate  results  also  suggest  that, 
for  efficient  drag  reduction,  the  solution  injected  should  be 
dilute  and  comparable  to  the  discharge  within  the  inner  boundary 
layer. 

Wu  and  Tulin  (1970)  performed  experiments  in  a  circulating 
water  channel  with  ejection  of  various  polymer  concentrations 
along  smooth  and  rough  plates  forming  the  cop  of  the  closed  test 


section.  It  was  recommended  chat,  for  che  nose  effective  drag 

reduction,  che  slot  ejection  angle  should  be  small  with  respect 

to  che  flow  direction  and  che  slot  opening  should  be  comparable 

to  che  thickness  of  che  viscous  sublayer.  Injection  races  were 

also  comparable  Co  the  normal  viscous  sublayer  discharge.  Optimum 

2  3 

additional  concentrations  were  found  co  be  10  -  10  *P?M  for  che 

smooth  place  and  an  order  of  magnitude  larger  for  rough  surfaces, 
where  mixing  due  Co  roughness  causes  Increased  dilution  of  che 
injected  solution. 

Poreh  and  Carmak  (1962)  studied  two-dimensional  turbulent 
mixing  of  ammonia  gas  from  a  line  source  near  a  wall.  They 
envisioned  a  four-zone  diffusion  process  that  is  worthy  of 
description  since  . che  concentration  profiles  are  considered 
representative  of  polymer  diffusion  and  for  this  reason  have  come 
under  much  study.  The  four  tones  as  described  by  Poreh  and 
Carmak  are: 

1.  An  initial  zone  -  This  tone  is  very  close  to  che  source 
in  which  a  large  portion  of  che  diffusion  boundary  layer  is  sub¬ 
merged  in  che  viscous  sublayer.  The  length  of  the  region  is 
determined  by:  the  initial  condition  near  the  source,  relative 
co  the  thickness  of  che  sublayer;  the  injection  velocity;  and 
che  magnitude  of  the  molecular  diffusivity.  The  extent  of 
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this  zone  was  not  determined.  Very  little  reliable  data  were 
obtained  in  this  region  due  to  very  large  velocity  and  concen¬ 
tration  gradients. 

2.  The  intermediate  zone  -  In  this  zone  the  diffusion  plume 
is  submerged  in  the  momentum  boundary  layer  and  its  thickness 
is  large  when  compared  to  that  of  the  sublayer.  The  diffusion 
rate  in  this  region  is  relatively  large  and  the  concentration 
profiles  are  found  to  be  approximately  similar  in  the  sense 
that 


C/CMax  *  f(ZA)’  (29) 

where  X  is  a  characteristic  height  of  the  diffusion  boundary 
layer  defined  as  the  distance  from  the  wall  where  C/C^ax  ■  0.5. 
Measurements  in  air  indicate  that  the  length  of  the  intermediate 
zone  is  20  to  40  boundary  layer  thicknesses  downstream  from  the 
source. 


3.  The  transition  zone  -  In  this  zone,  the  diffusion  rate  is 
slower  due  to  the  lower  level  of  turbulence  in  the  outer  portion 
of  the  boundary  layer  . 


4.  The  final  zone  -  The  growth  of  the  diffusion  boundary 


layer  coincides  with  that  of  the  momentum  boundary  layer  in 
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this  zone.  In  this  stage  the  maximum  concentration  near  the  wall 
is  inversely  proportional  to  the  thickness  of  the  boundary  layer 
and  the  ambient  velocity. 

Morkovin  (1963)  considered  the  experimental  evidence  of 
Poreh  and  Cermak  (1964)  for  the  final  two  zones  of  quasi-similar¬ 
ity  and  concluded  that  they  were  consistent  with  the  concepts  of 
eddy  diffusivity.  He  described  their  data  as  follows: 

Intermediate  zone:  C/C w  -  e"0,693  (y/6  )1,5  (30) 

Final  zone:  C/C  -  e~°'693  (y/5  )2-15.  (33.) 

W  Q 

Figure  (4)  displays  a  plot  of  the  concentration  profiles  in  the 
intermediate  and  final  zones. 

Poreh  and  Hsu  (1971)  indicate  that  the  most  widely  used 
method  to  predict  gross  diffusion  patterns  in  turbulent  flow 
employs  an  eddy  diffusivity  model.  This  model  assumes  that  the 
flux  of  the  diffused  matter  by  turbulent  fluctuation  is 
proportional  to  an  eddy  diffusivity  term,  Dg,  times  the  local 
concentration  gradient : 

0  »  -  D  (3c/3y) .  (32) 

y  e  J 
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igure  i.  Concentration  profiles  ia  the  iatemeciac 
and  final  rones 
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It  is  further  assumed  that  D  is  a  function  of  the  flow  field  and 

e 

its  value  can  be  specified  at  a  pcint  regardless  of  the  position 
of  the  source.  They  point  out  that  chic,  assumption  holds  true 
only  at  distances  from  the  source  that  are  large  compared  to  the 
Lagrangian  integral  scale  of  turbulence.  Measurements  made  by 
Poreh  and  Cermak  (1964)  indicate  that  this  limitation  holds  for 
diffusion  in  turbulent  shear  flows.  Estimates  of  the 
Lagrangian  integral  scale  suggest  that  it  is  of  the  order  of 
boundary  layer  thicknesses. 

Another  method  for  treating  diffusion  patterns  is  based  on 
Batchelor's  (1957)  Lagrangian  similarity  hypothesis  used  to 
predict  the  turbulent  motion  of  particles  in  steady,  self- 
preserving  shear  flows.  Cermak  (1962)  applied  the  Lagrangian 
similarity  hypothesis  to  predict  diffusion  from  a  continuous 
point  and  line  source  with  the  conclusion  that  results  from  the 
application  of  the  Langrangian  similarity  hypothesis  were 
significant  for  the  modeling  of  diffusion. 

Hsu  (1971),  and  Poreh  and  Hsu  (1971)  applied  these  techniques 
to  predict  the  diffusion  boundary  layer  growth  in  the  intermediate, 
transition,  and  final  zones  for  polymer  flow.  The  following 
equations  that  resulted  describe  the  change  of  the  mean  vertical 
position,  y,  and  the  mean  longitudinal  position,  x,  for  an  ensemble 
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of  single  particle  releases,  in  the  logarithmic  portion  of  the 
boundary  layer: 


dv 

dt 


* 


-  bV 


(33) 


where  b  is  Batchelor's  constant,  and 


dx 

dt 


u(y)  . 


(34) 


It  follows,  therefore,  that 


dv 

dx 


(35) 


Poreh  and  Hsu  (1971)  further  conclude  that  the  boundary  concen¬ 
tration  downstream  from  a  continuous  line  source  may  be  given  by 


C 


Max 


-S- 

y  u(y) 


(36) 


where  0  is  the  discharge  of  the  source  per  unit  width. 


Ellison  (1959)  estimated  that  Batchelor's  constant  is  given 
when  b  *  K,  where  K  is  the  Von  Karman  constant,  1/k  *  a  ,  in 
equation  (39).  His  analysis  further  suggests  that  y  (the  mean 
position  of  particles  at  a  given  cross  section,  x)  is  equal  to  Y, 
(the  mean  position  of  single  particle  releases  when  x  »  x) .  The 
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mean  position  of  a  particle  at  any  cross  section,  x,  may  be 
defined  as: 


/”cy  dy 

7  -  - -  •  (37) 

j‘  (c) dy 


By  replacing  y  with  y/o'd  where  is  the  value  of  y  when  C/C^ 
is  .5,  and  C  by: 


‘  £  <y/5d> 

aquation  (37)  may  be  integrated  to  obtain 


(38) 


7 


az  V 


(39) 


Substituting  equation  (39)  into  (35)  results  in  an  expression  for 
the  diffusion  boundary  layer  developing  with  distance  x,  such  that 


a 


i 


dx 


(40) 


Hsu  (1971)  observed  improved  agreement  with  data  when  re¬ 
placing  the  constant  b  with  the  expression 


b  -  K(l-y/o). 


(41) 


Figure  5  presents  a  comparison  of  experimental  data  with 
results  from  aquation  (40)  for  b  »  S,  b  ■  .3R,  and  b  «■  K(l-y/f). 
The  growth  of  the  diffusion  boundary  layer  within  the  momentum 
boundary  through  several  zones  of  diffusion  is  also  shown  in  this 
figure. 

Sampson  (1969)  constructed  a  recirculating  water  tunnel  and 
an  instrumentation  system  for  investigation  of  turbulent, 
non-Newtonian,  flat  plate  boundary  layers.  Mon-Newtonian 
velocity  profile  data  were  obtained  by  laser  velocimeter  at 
Reynolds  numbers  in  the  region  of  10^;  the  1/7 th  power  law  velocity 
profile  was  found  to  correlate  well  with  the  data. 

la  his  studies  with  concentrations  of  100  WPPM  of  the  polymer 
Separan  AP-30,  Rudd  (1972)  found  the  drag  reduction  mechanism 
confined  to  the  viscous  sublayer  close  to  the  wall.  The  polymer 
did  not  appear  to  have  any  significant  effect  on  the  central 
region  of  flow.  A  test  section  of  one-half  inch  square  pipe, 

7  feet  6  inches  long,  was  used  in  a  recirculating  flow  system 
driven  by  a  peristaltic  pump. 

Velocity  measurements  were  taken  with  a  laser  velocimeter. 

It  was  found  that  at  concentrations  of  100  WPPM,  the  polymer 
molecules  produced  quite  a  significant  amount  of  light  scattering, 


requiring  no  additional  seeding. 


figure  Crowlli  ol  tbe  diffusion  boundary  layer  wllliJn  the  momentum  boundary  layer 

(from  Mali  (1971)) 
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Logan  (1972),  using  a  one-half  inch  square  duct  under  gravity 
feed,  confirmed  Rudd's  findings  of  increased  streamwise  turbulence 
intensities  and  decreased  spanwise  intensities. 

Using  the  recirculating  tunnel  built  by  Sampson  (1969),  Kumor 
and  Sylvester  (1973)  conducted  a  series  of  experiments  on  the 
time  effects  of  polymer  drag  reduction  degradation  on  mean 
velocity  profiles.  Sparan  AP-30  was  the  polymer  used.  The 
experiments  show  that  the  velocity  gradient  at  the  wall  increases 
with  time,  while  the  sublayer  thickness  decreases. 

Reischman  (1973)  attempted  to  correct  the  flow  facility  de¬ 
ficiencies  of  previous  investigations.  He  constructed  a  unique 
dimensional  turbulent  flow  channel,  in  which  the  side  walls  were 
slightly  bowed,  to  allow  good  spatial  resolution  at  the  wall  in 
the  center  of  the  channel.  The  flow  channel  was  constructed  of 
1/4- inch  Plexiglas  with  a  cross  section  of  1  inch  by  12  inches 
(an  aspect  ratio  of  12:1)  and  a  length  of  70  inches.  Velocity 
measurements  in  homogeneous  100  WPPM  solution  of  Magnifloc  837-A, 
Sepman  AP  273,  and  Polyox  WSR-301  were  made  using  a  laser  Doppler 
anemometer  measuring  individual  realizations.  For  drag  reduc¬ 
tions  of  40%  or  less  for  the  near  wall  velocity  measurements, 
the  viscous  sublayer  (indicated  by  a  linear  velocity  profile) 
was  found  to  be  of  the  same  extent  as  for  a  comparable  solvent 
flow.  Reischman  found  no  evidence  of  a  thickened  sublayer. 
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Drop  Tests,  Rotating  Disks  and  Cylinders 

Crawford  and  Pruitt  (1963)  were  among  the  first  to  observe 
the  increased  velocity  of  steel  spheres  falling  in  a  solution  of 
Guar  Gum  as  compared  to  their  velocity  in  plain  water.  Ruszczycky 
(1965) ,  using  steel  spheres  of  a  diameter  up  to  1  inch,  performed 
tests  in  concentrated  solutions  of  Guar  Gum  and  Polyox  WSR-301  up 
to  15,000  WPPM.  For  the  1-inch  spheres,  the  maximum  drag  reduction 
observed  was  approximately  28%  at  Guar  Gum  concentration  of  5000 
WPPM  and  a  Polyox  concentration  of  7500  WPPM. 

Lang  and  Patrick  (1966)  found  that  the  drag  of  spheres 
falling  in  Polyox  WSR-301  was  reduced  by  about  70%  while  the 
drag  of  other  shapes,  such  as  cones  and  cylinders,  was  only 
slightly  reduced.  Dyed-wake  photographs  of  falling  spheres  in¬ 
dicated  a  significant  change  in  the  shape  of  the  wake,  with  the 

3 

laminar  separation  point  being  shifted  rearward  (Re  range  6  X  10 
to  2.5  X  105). 

Hoyt  et  al.  (1965,  1964)  showed  that  addition  of  Guar  Gum 
lowers  the  torque  by  as  much  as  70%  for  concentrations  of  300  to 
400  WPPM  for  turbulent  rotating  disk  flow. 

Using  a  multiple  disk  apparatus,  Whitsitt  et  al.  (1969) 
showed  that  rough  disks  gave  higher  torque  reduction  of  about  60% 


45 


in  Guar  Gum  solution  at  250  WPPM  than  did  smooth  disks.  A  poly¬ 
acrylamide  solution  lost  its  torque- reducing  ability  very  quickly. 
However,  the  work  of  Huang  and  Santelli  (1972)  showed  that  poly¬ 
acrylamide  solutions  were  more  resistant  to  shear  degradation 
when  compared  to  polyethylene  oxide. 

Sirmalis  (1974)  conducted  an  extensive  study  of  drag  reduc- 

o  0 

tion  on  axisymmetric  torpedo-shaped  bodies  having  6  and  12 
spherical  tail  cones.  During  the  study,  experiments  covering  a 
range  of  Reynolds  numbers  from  one  to  five  million  in  plain  water 
and  polymer  oceans  of  1  to  60  WPPM  of  Polyox  WSR-301  were  conducted. 
Photographic  studies  revealed  that,  at  very  low  concentrations, 

(2.5  WPPM),  the  fine  structured  turbulence  was  absent  and  only 
large  scale  turbulence  remained.  Computer  routines  tended  to 
overpredict  boundary  layer  thickness  while  indicating  proper  growth 
shapes.  No  shift  in  the  boundary  layer  separation  region  was  noted. 

Testing  in  fresh  water  with  polymer  ejection  from  the  body 
showed  similar  results  but  with  better  visual  representation  at 
the  higher  polymer  concentration.  Dyed  water  ejection  tests 
displayed  the  characteristic  course  and  fine  turbulence  structure. 
Ejection  of  50  WPPM  dyed  polymer  solution  eliminated  most  of  the 
fine  structured  turbulence  leaving  only  the  course  structure. 

Dye  streaking  was  apparent  at  the  higher  concentration  levels. 

The  streak  spacing,  consistent  with  laminar  sublayer  streaks. 
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was  also  in  agreement  with  the  number  of  ejection  holes  in  the 
screen  ejector  in  the  nose  of  the  body.  A  of  702  reduction 

in  skin  friction  was  achieved  with  the  6°  tail  body  in  a  polymer 
ocean  of  20  WPPM. 


Polymer  diffusion  tests  were  conducted  at  a  fixed  rate  of  20.6 
in/sec  and  a  nominal  velocity  of  27  ft/sec.  Wall  samples  and 
boundary  layer  samples  at  .025  inch  and  .055  inch  were  taken 
and  analyzed  by  fluorometric  techniques.  Tests  were  conducted 
with  water  and  concentrations  of  5,  10,  20,  50,  500,  and  1000  WPPM 
of  polymer.  All  polymer  tests  displayed  a  drastically  reduced 
diffusion  resulting  in  an  extended  initial  mixing  zone.  Similarity 
concentration  profiles  were  found  to  have  like  forms,  and  therefore 


C/C 

w 


■0.693 


(y/sd)K3. 


(41) 


Profile  coefficient,  has  an  experimentally  determined 
value  of  0.75.  Sirmalis  indicates  that  the  true  value  is  believed 
to  be  below  this  value,  since  data  were  not  obtained  in  the  crit¬ 
ical  region  of  within  approximately  0.010  inch  of  the  wall.  These 
results  were  considered  significant  in  that  they  imply  substan¬ 
tially  reduced  polymer  quantities  may  be  required  with  proper 
ejection  techniques. 


The  vast  number  of  investigations  concerning  the  phenomena 
of  drag  reduction,  and  the  continuing  attempts  to  conduct  a 
definitive  experiment,  indicate  the  scientific  community's  interest 
in  the  problem  and  the  difficulties  involved.  From  these  studies, 
it  is  known  that  drag  reduction  additives  include  soap  solutions, 
algae,  plant  derivatives,  and  high  molecular  weight  polymers. 
Additionally,  a  maximum  drag  reduction  asymptote  of  approximately 
80%  has  been  shown.  In  spite  of  these  efforts,  the  exact  nature 
of  the  phenomena  of  drag  reduction  remains  elusive  to  the  various 
techniques  applied  to  its  study. 

Flow  visualization  experiments  have  provided  some  insight 
from  a  qualitative  view  point,  while  pressure  drop-flow  rate 
experiments  have  contributed  an  indication  of  the  magnitude  and 
scope  of  the  phenomena.  The  bulk  of  these  efforts  are,  however, 
still  contradictory  and  difficult  to  interpret  due  to  the  confusing 
features  of  the  flow  facilities  and  ambiguities  in  the  measurement 
instrumentation  employed. 

Reischman  (1975)  points  out  that,  in  general,  experimental 
techniques  for  velocity  measurement  may  be  separated  into  four 
groups:  pitot  tubes,  bubble  tracing,  hot-film  probes,  and  laser 


Doppler  anemometers. 
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PitoC  probe  investigations  comprise  the  largest  group  and  are 
subject  to  errors  requiring  complex  corrections.  Smith  et  al. 
(1967)  have  demonstrated  that  measurements  in  identical  flow 
situations  made  by  various  size  pitot  probes  yield  different 
results.  Metzner  and  Astaria  (1967)  attribute  these  errors  to 
the  influence  of  additional  viscoelastic  normal  stress  terms. 
Corrections  for  pitot  probe  techniques  are  generally  omitted, 
however,  as  being  too  complicated.  In  spite  of  this,  Tomita 
(1970)  applied  viscoelastic  corrections  to  his  pitot  probe 
measurements . 

Hot  film  probes,  the  next  largest  group,  depend  upon  the 
heat  transfer  characteristics  of  the  medium.  Frieche  and 
Schmartz  (1969)  have  shown  that  hot-element  sensors  are  difficult 
to  calibrate  because  polymer  additives  alter  the  heat  transfer 
characteristics  of  the  medium  and  collect  on  the  sensors.  It  is 
pointed  out  that  the  calibration  drift  problem  is  particularly 
severe  in  dilute  polymer  solutions  where  the  sensitivity  of  the 
probe  to  velocity  changes  is  lower  than  it  is  in  the  water  alone. 

Donohue  et  al.  (1972)  point  out  that  bubble  tracing 


techniques  are  hampered  by  large  uncertainties  as  well  as  being 
an  extremely  tedious  process. 
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The  fourth  group  comprise  laser  Doppler  anemometer  efforts. 

A  relatively  new  discovery  by  Yah  and  Cummings  (1964) ,  the  laser 
Doppler  anemometer  is  a  linear,  non- interfering  instrument  which 
does  not  depend  on  che  rheological  or  intensive  properties  of  the 
working  fluid.  The  laser  anemometer  was  first  used  in  drag 
reducing  flows  during  four  investigations  almost  simultaneously: 
Chung  and  Graebel  (1969),  Goldstein  et  al.  (1967),  Rudd  (1969), 
and  Shankan  (1969).  The  measurements  of  Shankan  (1969)  have 
insufficient  sample  si2e,  however,  for  accurate  results,  while 
Goldstein  et  al.  (1967)  only  made  measurements  on  the  tube  center- 
line. 


Chung  and  Graebel  (1969)  experimented  in  a  small,  .47-inch 
diameter  pipe  where  che  spatial  resolution  of  their  laser 
anemometer  was  very  poor.  Subsequently,  five  researches  have 
made  acceptable  laser  velocimetry  measurements  in  drag-reducing 
flows.  This  group  consists  of  Rudd  (1969),  Logan  (1972), 
Sampson  (1969),  Kumor  and  Sylvester  (1973),  and  Reischman 
(1975).  Of  these  experiments,  all  were  conducted  in  a  square 
cross-section  apparatus,  except  the  work  of  Reischman  who  used 
a  12:1  aspect  ratio  tunnel.  Rudd  and  Logan  utilized  a  one-half 
inch  square  duct  and  Sampson,  and  Kumor  and  Sylvester,  used  an 
8-inch  square  tunnel  with  a  submerged  off-center  flat  plate. 
White  law  (1973)  has  shown  that  non-symmetric  secondary  flows  are 
significant  in  square  duct  flows. 


Sampson  (1967)  reported  that  pump  vibrations  were  signifi¬ 
cant,  prohibiting  meaningful  tests  at  higher  flow  rates  and  that 
the  effect  of  these  vibrations  on  the  laser  system  were  unknown. 
Additionally,  accommodations  in  both  equipment  and  system 
design  were  necessary  to  minimize  costs  of  the  construction  of 
the  laser  system.  Although  the  electronic  output  was  noisy  with 
a  relatively  poor  signal- to-noise  ratio,  usable  data  were 
obtained  by  use  of  a  storage  oscilloscope. 

Kumor  and  Sylvester  utilize  the  same  facility  as  Sampson  but 

offer  no  water  data  to  allow  verification  of  the  suitability  of 

this  facility. 

Reischman  avoided  the  detrimental  recirculation  aspects  of 
the  previous  researchers  by  using  his  flow  facility  in  a  blow 
down  mode.  In  contrast  to  the  previous  researchers,  Reischman  used 
an  individual  realization-type  laser  Doppler  anemometer,  which 
measures  the  period  for  10  cycles  of  a  Doppler  burst  scattered 
from  a  single  particle.  However,  the  orientation  of  the  laser 
light  beam  optical  system  was  in  an  undesirable  side  scatter 
mode  where  scattered  light  intensity  is  predictably  weak.  The 
singular  results  of  Reischman,  Indicating  that  the  viscous  sub¬ 
layer  does  not  thicken,  were  obtained  in  an  improved  flow  facil¬ 
ity  under  weak  signal  conditions  for  homogeneous  polymer  flow. 
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This  Is  in  direct  conflict  with  the  conclusion  of  all  previous 
investigators. 

Rudd,  Logan,  and  Sampson  used  continuous  wave  anemometers 
employing  a  frequency  tracker  which  converts  Doppler  frequency  to 
an  analog  signal  proportional  to  the  instantaneous  velocity. 

To  date,  all  of  the  work  using  laser  Doppler  anemometers 
has  been  performed  with  homogeneous  polymer  flows,  not  with 
ejection  of  polymer  solutions.  A  goal  of  the  present  work  was 
to  improve  the  deficiencies  of  previous  researchers  and  obtain 
data  on  velocity  profiles,  with  and  without  injection  of  polymer 
solutions  into  the  boundary  layer  on  a  submerged  flat  plate  in 
a  high  aspect  ratio,  rectangular  water  tunnel.  Laser  Doppler 
anemometry  velocity  measurements  in  this  experiment  were  made 
in  the  forward  scatter  mode  where  scattered  light  intensity 
distributions  are  strongest. 


III.  SZPF3.IMENTAL  APPARATUS 


The  experimental  apparatus  in  this  research  consists  of  a 
flow  facility  and  water  tunnel,  a  polymer  ejection  system,  a 
laser  Doppler  anemometer  system  for  velocity  measurement,  a 
fluorometer  to  measure  boundary  layer  concentration  profiles, 
and  a  minicomputer  for  real-time  data  processing.  The  character¬ 
istics  of  each  element  are  discussed  in  detail  in  the  sections 
that  follow. 


Flow  Facility 


The  overall  flow  facility  designed  for  this  experimental 
effort  is  shown  schematically  in  figure  5.  Filtered  water  is  fed 
to  a  15,C00-gallon  capacity  water  tank  having  a  constant  head 
overflow,  providing  5  psi  flow  head  for  the  tunnel.  Various 
filter  cartridges  allow  filtration  of  particles  ranging  in  size 
from  5  microns  to  .5  micron.  A  high  pressure  steam  heat  exchanger 
allows  deaeration  of  the  water  by  heating.  The  flow  rate  of 
water  from  the  storage  tank  to  the  tunnel  is  measured  by  a  Cox 
Instruments  Co.  turbine  flowmeter  and  controlled  by  a  manual  ball 
valve  in  the  exit  line.  A  round  to  rectangular  transition  section 
provides  entrance  to  a  rectangular  diffuser.  The  flow  is  then 
straightened  in  a  settling  chamber  designed  to  form  small  scale 
turbulence  of  uniform  distribution.  The  flow  enters  the  water 
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tunnel  through  a  specially  designed  transition  nozzle.  The 
flat  plate  is  attached  to  the  floor  of  the  tunnel  allowing  wall 
samples  to  be  taken  through  the  bottom  wall  of  the*,  tunnel.  The 
boundary  layer  approaching  the  flat  plate  is  bled  off  through  a 
suction  port  ahead  of  the  leading  edge  of  the  flat  plate. 

Boundary  layer  sampling  probes  are  lowered  through  ports  in  the 
tunnel  ceiling.  Velocity  measurements  are  taken  throughout  the 
entire  internal  volume,  without  disturbing  the  flow  using  the 
laser  Doppler  anemometer  techniques  and  a  three-dimensional 
traversing  system.  The  flow  then  exits  the  tunnel  through  an 
exit  plenum.  Constant  tunnel  exit  pressure  is  provided  by  exiting 
the  flow  into  a  container  with  constant  exit  head  overflow.  The 
boundary  layer  bleed  flow  is  also  exited  into  this  container. 
Polymer  solutions,  for  ejection  over  the  flat  plate,  pass  through 
a  delivery  line  coiled  in  the  exit  flow  container.  The  temper¬ 
ature  of  the  polymer  solution  then  approaches  the  temperature  of 
the  water  flowing  through  the  tunnel.  Figure  7  is  a  photograph 
of  the  actual  test  setup. 

Initial  Attempts 

In  order  to  avoid  the  design  of  a  round  to  rectangular 
transition  diffuser  and  because  axial  length  was  at  a  premium 
in  the  test  setup,  a  “ilatively  large  cross-section  settling 
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chamber  was  designed  and  bullc  of  anodized  aluminum  as  shown 
in  figures  8  and  9.  The  dimensions  were:  width,  3.25  inches; 
height,  6.5  inches;  and  length,  11  inches.  The  top  and  rear  side 
were  made  of  Plexiglas  for  measuring  and  photographing  the  flow. 
Dye  probe  access  ports  were  located  in  the  top  and  the  flow 
entered  through  a  3/4-inch  fitting.  Two  overlapping  perforated 
plates  (0.25-inch  diameter  holes  and  43%  open  area)  were  thought 
to  distribute  the  flow,  followed  by  1/8-inch  3EXC2LL  honeycomb 
(1  inch  long)  for  flow  straightening  and  four  screens  (24-mesh, 
.007-inch  diameter  wire  of  67.42  open  area)  for  the  even  dis¬ 
tribution  of  fine,  small-scale  turbulence.  A  rounded  edge  slot 
formed  the  tunnel  entrance. 

The  perforated  plates  were  totally  ineffectual  as  a  high¬ 
speed  central  core  water  jet  maintained  itself  well  into  the 
settling  chamber.  An  impact  distribution  plate  was  attached  to 
the  inlet  fitting.  This  fix  transformed  the  high-speed  core  flow 
into  violent  turbulent  recirculating  flows.  In  addition,  corner 
vortices  formed  in  the  exit  flow,  as  shown  in  figure  10.  At  this 
juncture  it  became  evident  chat  a  diffuser  would  be  necessary  to 
decelerate  the  flow  and  an  entrance  nozzle  to  the  tunnel  would 
also  be  required  to  eliminate  the  formation  of  comer  vortices. 

A  redesign  of  the  settling  chamber  to  provide  equal  contraction 
ratio  in  the  horizontal  and  vertical  planes  was  also  performed  in 
an  attempt  to  further  minimize  the  formation  of  corner  vortices. 
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Settling  Chamber 

The  second  settling  chamber  design  dimensions  now  remained 
firm  at  a  width  of  8.25  inches,  a  height  of  2.75  inches,  and  a 
length  of  18.5  inches.  The  contraction  ratios  from  settling 
chamber  to  tunnel  were  fixed  at  3.67  in  both  the  horizontal  and 
vertical  planes,  providing  an  overall  contraction  of  13.5:1.  The 
settling  chamber  was  constructed  entirely  of  Plexiglas  and  had  an 
internal  configuration,  shown  (in  order)  in  figure  11,  consisting 
of  a  perforated  plate,  two  screens,  a  honeycomb,  and  four  addition¬ 
al  screens,  spaced  at  1-inch  intervals  with  2  inches  upstream  of 
the  honeycomb.  The  entrance  to  the  tunnel  was  a  well-rounded 
rectangular  hole,  shown  in  figure  12. 

The  configuration  straightened  the  flow;  however,  corner 
vortices  at  the  tunnel  entrance  remained  (figure  13) .  A  standoff 
rectangular  nozzle  was  designed  and  tested  with  no  success,  as 
shown  in  figures  14  and  15.  In  addition,  small  flow  oscillations 
appeared  due  to  intermittent  stall  in  the  diffuser. 

The  corner  vortices  were  finally  eliminated  by  the  design 
and  fabrication  of  a  transition  nozzle  in  the  form  of  a  ninth 
order  polynomial  whose  first  four  derivatives  are  zero  at  each 
end.  The  polynomial  coordinates  were  placed  on  paper  tape  and 
used  to  cut  templates  for  the  vertical  and  horizontal  contours  in 


Figure  13.  Settling  chamber  flow  pattern  photograph 
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a  cape-driven  automatic  milling  machine.  These  templates  were  used 
to  cut  a  mandril  from  a  cast  block  of  epoxy  resin.  This  mandril 
was  placed  in  a  casting  box  of  the  same  internal  dimensions  as  the 
settling  chamber.  The  transition  nozzle  block  was  then  cast  over 
the  mandril.  Removal  of  the  mandril  after  an  appropriate  curing 
time  yielded  a  block  with  proper  internal  flow  passage  configura¬ 
tion  that  fit  into  the  settling  chamber.  The  transition  nozzle 
used  is  shown  in  figures  16  and  17. 

Discussions  on  the  internal  configurations  of  the  settling 
chamber,  with  Drs.  Andrew  Chawat  and  Steven  Barker  at  the 

University  of  California  at  Los  Angeles,  provided  the  key  element 
for  the  removal  of  the  final  obstacle  for  flow  conditioning.  A 
3-inch  length  of  reticulated  Scott  foam  quieted  and  steadied  the 
flow,  removing  any  remaining  flow  disturbances. 

The  final  settling  chamber  design  consisting  of  a  perforated 
plate,  Scott  foam,  honeycomb,  two  screens  and  a  transition  nozzle 
is  shown  in  figure  18  and  schematically  in  figure  19. 

Diffuser 


A  36- inch  long  rectangular  Plexiglas  diffuser  (figure  20)  was 
designed  having  no  appreciable  stall,  according  to  Kline  (1959). 

The  entrance  dimensions  were  0.75  inches  by  0.75  inches  with  rectan- 
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gular  exit  dimensions  of  8.25  Inches  in  width  and  2.75  Inches  in 
height  to  match  the  new  settling  chamber  dimensions.  In  order  to 
minimize  separation  effects  In  going  from  . 6 25- inch  (inside 
diameter)  round  pipe  to  a  .75-inch  square  cross  section,  a  special 
transition  section  was  fabricated.  The  walls  at  the  entrance  of 
the  diffuser  were  thickened  with  extra  layers  of  Plexiglas  and  a 
round  entrance  hole  of  a  diameter  equal  to  the  diagonal  of  the 
entrance  was  drilled.  The  Interfaces  where  the  round  hole  broke 
through  the  interior  walls  were  hand-blended  smooth.  A  conical 
transistlon  section  from  the  .625-lnch  (Inside  diameter)  inlet 
pipe  to  the  rounded  diffuser  Inlet  was  made  and  is  shown  in  figure 
21.  A  screen  was  Incorporated  to  aid  In  maintaining  flow  distri¬ 
bution. 


Tunnel  Design 

The  tunnel  is  constructed  of  .375-inch  thick  Plexiglas  walls 
and  inverted  aluminum  channel  sections  forming  the  top  and  bottom 
of  the  tunnel.  The  interior  dimensions  are:  height,  .75  inches; 
width,  2.21  inches;  and  length,  24.25  inches.  Ten  access  ports 
along  the  top  of  the  tunnel  allow  insertion  of  boundary  layer 
sampling  probes  into  the  tunnel  flow.  The  five  stations  used  In 
this  experiment  are  shown  in  figure  19. 


Figure  21.  Conical  diffuser  Inlet  transition  section  photograph 
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Starting  1.75  inches  from  the  entrance  to  the  tunnel,  an 
anodised  aluminum  flat  plate  is  fixed  to  the  floor  of  the  tunnel. 

A  one-half  inch  suction  slot  at  the  upstream  edge  of  the  flat  plate 
allows  the  incoming  boundary  layer  to  be  bled  off  and  a  new  bound¬ 
ary  to  be  started  in  the  flat  plate.  An  axially  adjustable  lead¬ 
ing  edge  allowed  the  ejection  slot  height  to  be  varied  up  to  0.050 
inches  In  the  horizontal,  x,  direction.  Boundary  layer  vail 
concentration  sampling  holes  (.030  inch  in  diameter)  were  drilled 
at  locations  corresponding  to  the  probe  sampling  planes.  The 
polymer  ejection  slit  was  located  1.75  inches  from  the  leading 
edge  and  slanted  at  an  angle  of  15°  to  the  flaw. 

Eight  different  leading  edges,  shown  in  figure  22,  were  tried 
in  an  attempt  to  obtain  transition  in  the  center  portion  of  the 
plate.  Of  all  the  leading  edges  tried,  the  most  effective  one 
(figure  23)  consisted  of  an  epoxled  sand  surface  (average  grain 
size  0.020  inch)  sloped  at  7.5°. 

Boundary  Layer  Probes 

Boundary  layer  probe  rakes  were  fabricated  and  are  shown 
assembled  in  figure  24  and  as  an  exploded  view  in  figure  25.  The 
mounting  fixture  has  a  slot  machined  in  its  side  to  allow  the  rake 
to  be  retracted  up  into  the  roof  of  the  tunnel  out  of  the  flow 
paasage.  Figure  26  is  a  close-up  of  the  probe  tips  shoving  their 


Figure  24.  Boundary  layer  probe  rake  photograph 


capered  rectangular  shape  of  .005  inch  high  by  .015  inch  vide. 
Figure  27  presents  probe  dimensions  and  table  2  presents  the 
location  of  the  sampling  points  vith  and  vlthout  the  spacers  used 
to  reposition  probes  in  the  vertical,  7,  direction.  The  bound¬ 
ary  la7er  samples  would  pass  through  the  probe  and  along  a  36-inch 
b7  .020-inch  diameter  flexible  tube  to  a  sampling  stand  containing 
six  cuvettes.  Figures  28  and  29  show  this  arrangement  with  reser¬ 
voirs  in  place  of  the  cuvettes  in  the  sampling  stands. 

Polymer  Inlection  System 


The  polTmer  injection  S7stem  was  designed  to  minimize  polTmer 
degradation  during  the  delivery  process.  As  shown  in  figure  6, 
water  is  recirculated  to  a  water  reservoir.  Part  of  the  water 
flow  is  bypassed  through  a  water-calibrated  rotometer  to  the  poly¬ 
mer  reservoir.  The  polymer  reservoir  contains  13  liters  of  polymer 
solution  above  a  separation  piston.  The  water  flow  pushes  the 
piston  upwards.  This  action  in  turn  pushes  the  polymer  mixture 
out  of  the  reservoir  into  the  delivery  line  to  the  tunnel.  The 
delivery  line  is  coiled  in  the  constant  exit  head  container  and 
submerged  in  a  bath  of  tunnel  exit  flow  water  to  equalize  the 
polymer  flow  temperature  with  that  of  the  main  tunnel  flow. 

Figure  30  i3  a  photograph  of  Che  polymer  delivery  system. 
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TABLE  2 

CONCENTRATION  SAMPLING  LOCATIONS  IN  THE  BOUNDARY  LAYER 


Probe 

Sample 


i 


2 


3 


5 


6 


7 


8 


Yc  Yc 

(Spacer  .020-inch)  (Spacer  .103-inch) 


0 


.010 


.030* 


.072 


.092* 


.134 


.154* 


.196 


.216* 


0 


.010 


.072 


.  113* 


.134 


.175 


.196 


.237* 


.299* 


i 


*  Spacer  in  place 
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Figure  29.  Boundary  layer  probe  sampling  system  photograph 


Polymer  deliver 


Boundary  Layer  Concentration  Measurements 

The  methods  used  for  the  measurement  of  concentration  distri¬ 
butions  of  polymers  in  solution  include:  polarographic  (investi¬ 
gated  by  Goren  (1965));  turbidity  (investigated  by  Ruch  et  al. 
(1965));  and  opaqueness  and  fluorometric  techniques.  Wetzel  and 
Ripkln  (1970)  present  a  complete  discussion  of  all  but  the 
opaqueness  method,  which  was  used  by  Fruman  and  Tulin  (1974). 

Fruman  and  Tulin  (1974)  constructed  a  specially-built  light 
absorption  system  utilizing  photocells  and  a  recording  oscillograph. 
The  intensity  of  the  light  source  was  chosen  in  order  to  achieve 
the  maximum  possible  output  readings  between  zero  dye  and  a  dye 
concentration  corresponding  to  practically  zero  light  transmission. 
This  dye  concentration  was  chosen  to  be  1250  WPPM  of  India  ink  in 
tap  water.  The  voltage  output  of  the  photocell  was  40.0  mV  and 
the  equivalent  spot  displacement  on  the  recorder  was  12.0  cm. 

When  a  continuously- running  sample  provided  a  steady  state  output, 
the  reading  was  compared  to  previous  readings  obtained  during  the 
calibration  of  the  system.  The  scatter  in  the  data  was  attributed 
to  stability  problems  with  the  light  Intensity  and  galvanometer 
of  the  recording  oscillograph. 

The  method  selected  for  use  in  this  study  is  the  fluorometric 


method.  This  method  has  been  used  successfully  by  Wetzel  and 
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Ripkin  (1970) ,  Walters  and  Wells  (1971) ,  Tullis  and  Ramu  (1973)  , 
and  Sirmalis  (1974).  The  method  consists  of  injecting  a  tracer 
dye  into  the  fluid  being  analyzed  and  capturing  a  sample  of  the 
fluid  for  analysis  in  a  fluorometer. 

Fluorescence,  the  basic  phenomenon  being  measured,  is  the 
instantaneous  emission  of  light  from  a  molecule  or  atom  that  has. 
absorbed  light.  Consideration  of  the  equation  of  fluorescence  in 
dilute  solution  shovs  that 

F  -  KI  £Cd0  (42) 

Qi  O 

where 

F  is  the  meter  reading, 
m 

R  is  the  amplification  factor  of  the  photodetector, 

IQ  is  the  intensity  of  the  exciting  light, 
e  is  the  extinction  coefficient  of  the  compound 

(a  constant  for  any  given  wavelength  of  I  ) , 

o 

C  is  the  concentration  to  be  measured, 
d  is  the  path  length  (a  constant) , 

0  is  the  quantum  yield  (a  constant) , 
indicates  that  decreasing  concentration,  C,  may  be  compensated  by 
increasing  the  amplification  factor,  K,  or  the  light  intensity,  IQ. 
This  is  done  in  fluorometry,  yielding  larger  scale  readings  for 
lower  concentration.  The  fundamental  principle  of  fluorescence 
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measurement  is  shown  in  figure  31,  a  simplified  schematic 
of  a  fluorometer  taken  from  Turner  (1972).  The  intensity  of  the 
exciting  light  is  controlled  by  a  range  selector  having  four 
apertures  or  sensitivity  ranges.  The  desired  wavelength  of 
exciting  light  is  selected  by  a  primary  monochromator  or  filter 
placed  between  the  light  source  and  the  sample.  Emitted  by  the 
sample,  the  wavelength  of  light  to  be  measured  is  selected  by 
a  second  filter  placed  between  the  sample  and  a  photodetector. 

The  output  of  the  photodetector,  a  current  proportional  to  the 
intensity  of  the  fluorescent  light,  is  amplified  to  give  a  reading 
on  a  meter  or  recorder. 

In  operation,  a  standard  sample  is  placed  in  the  fluorometer 
and  the  sensitivity  or  intensity  of  exciting  light  adjusted  to  a 
desired  reading.  Unknown  samples  and  a  blank  are  then  read.  The 
net  readings  of  the  standard  and  unknown  samples,  with  the  blank 
subtracted,  are  in  the  same  proportion  as  their  concentration. 
Standard  linear  curves  may  be  prepared  and  the  concentrations  of 
samples  read  directly  from  the  curve. 

The  basic  assumption  made  in  all  polymer  solution  studies  is 
that  the  diffusion  of  the  tracer  dye  is  identical  to  that  of  the 
boundary  layer  mixture.  This  is  a  reasonable  assumption  based 
on  data  comparing  various  measurement  methods  by  Wetzel  and  Ripkin 
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(1970).  Tests  of  tvo  widely  used  fluoroscent  trace  dyes,  Rhodimine 
and  Uranine-B,  found  both  to  be  compatible  with  Polyox.  The 
solutions  were  stable  and  their  drag-reducing  properties  were 
unaffected  when  tested  In  a  capillary  tube  blowdown  rheometer 
similar  to  Hoyt  (1966) . 

Selection  of  the  fluorometric  method  for  the  present  experi¬ 
ment  was  based  on  the  accuracy,  ease,  and  rapidity  In  determining 
dye  concentrations  reported  by  previous  authors  and  the  availabil¬ 
ity  of  a  G.K.  Turner  Model  111  Fluorometer,  figure  32.  The  use  of  a 
fluorometer  of  the  type  utilized  by  previous  Investigators  should 
minimize  data  correlation  problems  due  to  apparatus  and  technique 
differences. 

Although  it  is  easy  to  use,  the  Rhodamine  dye  was  rejected 

because  it  is  very  difficult  to  clean  off  the  apparatus.  The 

Uranine-B  dye  was  selected,  however,  since  it  is  readily  adaptable 

9 

to  laboratory  use  and  gave  a  reproducible  reading  to  1  part  in  10 
on  a  weight  basis.  This  range  was  adequate  as  the  dye  concentration 
of  the  injected  solution  could  be  adjusted  to  produce  samples 
along  the  plate  sufficiently  above  the  background  blank  levels  to 
give  reproducible  results.-  Calibration  curves  were  experimentally 
developed  for  the  fluorometer  using  known  concentrations  of 
Uranine-B  dye  in  the  type  water  to  be  used  in  the  experiment. 


The  G.K.  Turner  Model  111  Fluorometer  has  four  sensitivity 
aperture  settings  allowing  calibration  curves  to  be  drawn  over  a 
range  of  about  1  part  per  10^  to  1  part  per  10^  on  a  weight  basis. 
For  concentrations  greater  than  the  lower  value,  controlled 
dilution  prior  to  measurement  was  performed.  The  calibration 
curves  used  in  this  experiment  are  shown  in  figures  33  and  34. 
Sample  sizes  of  approximately  3.5  to  4.5  ml  are  required  for 
measurement. 

The  sample  is  placed  in  a  small  test  tube  called  a  cuvette 
and  inserted  in  the  instrument.  The  measurement  is  automatically 
taken  when  the  door  is  closed.  Cleanliness  must  be  observed  and 
care  taken  to  wipe  the  outside  of  the  cuvette  prior  to  measurement 
as  fingerprints  affect  the  reproducibility  of  the  sample  readings. 
Wiping  the  exterior  of  the  cuvette  with  a  clean  paper  towel  was 
found  to  be  sufficient.  All  samples  were  allowed  to  achieve  room 
temperature  prior  to  reading. 

Table  3  provides  a  list  of  apparatus  used  in  the  concentra¬ 


tion  measurements. 
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FLUIO  CONCENTRATION  URANINE  -3  OYE 


Figure  33.  Calibration  curves  for  model  Ill  fluoromecar 

opening  1  and  10 
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TABLE  3 

CONCENTRATION  MEASUREMENT  APPARATUS 


Component 

Specification 

Fluor ome ter 

G.  K.  Turner  Model  111 

Excitation  lamo 

G.E.  #  G4T4/1 

Primary  filter 

2A  (closet  to  lamp) 

47B 

Together  these  filters  trans¬ 
mit  only  light  at  435  nm 
wave  length. _ 


Secondary  filter 

2A-12 

This  filter  only  passes  light 

at  wavelengths  greater  than 

500  nm. 

Cuvettes 

Matched  cuvettes 

Fisher  Scientific  No. 

14-385-900A 

Dye 

Uranine-B 

Fluoresces  at  470  nm  and 

515  nm 
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Laser  Doppler  Anemometer  System 

The  velocity  measurements  in  this  experiment  were  taken  with 
a  laser  Doppler  anemometer.  (LDA  will  be  used  as  short  form 
terminology.  Operation  of  an  LDA  system  is  based  on  the  prin¬ 
ciple  that  coherent  laser  light  scattered  from  a  particle  in  a 
moving  fluid  will  be  Doppler  shifted  in  frequency.  The  basic 
components  of  an  LDA  system,  shown  in  figure  35,  include  a  laser, 
transmitting  optics,  receiving  optics,  photodetector,  and  signal 
processing  equipment.  Table  4  presents  a  list  of  the  laser 
Doppler  system  components  used  in  this  experiment. 

Among  the  advantages  of  LDA  systems  are  the  following: 

1.  The  flow  is  not  disturbed 

2.  No  flow  calibration  is  required 

3.  There  is  excellent  spacial  resolution  due  to  small 
measuring  volume 

4.  High  frequency  response  occurs. 

Constraints  on  the  use  of  LDA  systems  require  (1)  an  optically 
clear  medium,  (2)  a  window  to  the  flow  (two  windows  are  preferred), 

/ 

and  (3)  particles  in  the  flow.  Because  the  LDA  measures  particle 
velocity  and  not  fluid  velocity,  there  must  be  sufficient  concen¬ 
tration  of  particles  of  the  proper  size  and  density  to  provide  the 
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TABLE  4 

LASER  DOPPLER  ANEMOMETER  SYSTEM  COMPONENTS 


Component 

Specification 

•  Laser : 

•  Transmitting  Optics: 

Lens  I: 

Lens  II : 

Lens  III: 

•  Beam  separation: 

•  Receiving  Optics: 

Lens  I : 

Lens  II: 

Lens  III: 

•  Photodetector: 


Spectra  Physics  model  120  Helium- 
Neon  laser  power  -  5  milliwatts 
beam  diameter  -  .8  millimeters 
(at  1/e^  point)  wavelength  - 
632.3  x  10  meter3 

TSI  Model  910 

Model  917 

Focal  length  103.  mm 

Beam  intersection  angle  25.0 

Model  918 

Focal  length  241.  mm 

Beam  intersection  angle  11.58 

Model  919 

Focal  length  598.  mm 

Beam  intersection  angle  4.71 

50  mm 

TSI  Model  930 
Model  937 

Focal  length  102.  mm 
Model  938 

Focal  length  242.2  mm 
Model  939 

Focal  length  577.2  mm 

TSI  Photodiode  model  960 
3andwidth  100  Mhz 


•  Aperture  diameter:  0.25  dia  mm 

•  Signal  Processor:  TSI  tracker  model  1090 


Frequency  response: 


2kHz  to  50  MHz 
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required  signal  while  closely  following  turbulent  fluctuations  in 
the  flow. 

In  general,  LDA  systems  are  characterized  by  the  arrangement 
of  their  optical  components  as  shown  schematically  in  figure  36. 

In  the  "forward  scatter"  arrangement,  the  receiving  optics 
and  photodetector  are  located  opposite  the  laser  to  receive  light 
scattered  in  the  forward  direction.  This  is  the  most  favorable 
mode  as  the  majority  of  incident  light  is  scattered  by  the 
particles  in  the  forward  direction.  A  typical  polar  scattered 
light  intensity  distribution  predicted  by  the  Mie  scattering 
theory  for  micron  size  particles  is  shown  in  figure  37.  The 
scattered  lighc  intensity  in  the  forward  direction  is  typically  of 

3 

order  10  times  greater  than  the  backscattered  light  Intensity. 
Requiring  two  windows  on  the  flow,  the  forward  scatter  mode  is 
ideal  for  water  tunnel  measurements.  Lower  laser  power  and  par 
tide  concentrations  may  be  used  due  to  increased  light  intensity. 

The  "backscatter"  arrangement  has  the  receiving  optics  and 
photodetector  located  on  the  same  side  of  the  flow  as  the  laser 
light  source.  This  arrangement  requires  only  one  window  on  the 
flow.  However,  higher  laser  power,  shorter  focal  lengths,  in¬ 
creased  particle  concentration  and/or  size  are  required  due  to 
decreased  light  intensity  scattered  in  the  backwards  direction. 


"Reference  beam"  mode  refers  to  a  single  observed  light  beam, 
whereas  the  "dual  beam"  mode  refers  to  the  system  where  the  light 
from  two  crossed  beams  is  observed.  Although  the  reference  beam 
mode  in  forward  scatter  was  developed  first,  the  dual  beam  mode 
will  be  discussed  in  detail  here  as  it  is  the  actual  mode  used 
and  the  easiest  to  conceptualize. 

In  the  dual  beam  mode,  the  laser  beam  is  split  into  two  beams 
of  equal  intensity,  equidistant  from  the  original  beam.  The  beams 
are  then  focused  at  a  crossover  point  where  they  interfere  with 
each  other  to  form  a  fringe  pattern,  shown  schematically  in  figure 
38.  These  fringes  are  caused  by  the  interaction  of  the  light  in 
the  two  beams  cancelling  and  reinforcing  each  other  to  form  a 
standing  electromagnetic  interference  pattern.  The  crossover 
region  where  the  two  beams  are  focused  is  called  the  measuring 
volume.  Figure  39  is  an  actual  photograph  of  a  plane  of  this 
measuring  volume  magnified  600  times,  containing  over  100  fringes. 
A  particle  moving  through  the  measuring  volume  in  the  plane  of 
the  two  beams  then  passes  through  alternating  regions  of  low  light 
intensity  (light  cancelling)  and  regions  of  high  light  intensity 
(light  reinforcing) . 

Light  scattered  from  particles  passing  through  the  measuring 
volume  is  focused  by  the  receiving  optics  onto  a  photodetector 
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Figure  38.  Fringe  pattern  schematic 


Fringe  pattern  photograph 
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whose  output  signal  is  shown  schematically  In  figure  40.  An 
aperture  is  placed  in  front  of  the  photodetector,  blocking  light 
from  locations  other  than  those  from  the  beam  crossover  region. 

Due  to  the  Gaussian  intensity  distribution,  the  amplitude  will 
vary  with  a  maximum  occurring  in  the  center  of  the  measuring  volume. 
The  frequency,  fp,  of  this  signal,  produced  by  the  scattered  light 
pulses  from  a  particle  travelling  through  the  measuring  volume, 
is  linearly  proportional  to  the  velocity  with  which  the  particle 
is  moving  across  the  fringes  (figure  41).  This  frequency  is  de¬ 
termined  by  the  frequency  tracker  of  the  LDA  system.  The  distance 
between  fringes,  d^,  is  a  function  of  the  laser  wavelength  and  the 
angle  of  beam  intersection  only,  such  that 

.  1 

r  ”  2  sin(0/2) ’  (43) 

The  velocity  of  the  particle  may  then  be  given  by 

0  -  f  •  df .  (44) 

Conversion  of  frequency  to  velocity  requires  only  that  the 
laser  wavelength  and  angle  of  beam  intersection  be  known.  A  flow 
direction  ambiguity  of  ISO0  does  exist,  however,  in  that  the  fringe 
system  cannot  distinguish  reverse  flow  from  forward  flow.  The 
problem  may  be  resolved  for  reversing  flow  by  frequency  shifting. 


Figure  40.  Photodetector  output  schematic 
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The  fringes  may  he  considered  to  be  moving  if  one  beam  is  fre¬ 
quency  shifted  relative  to  the  other.  If  the  difference  in 
frequency  between  the  two  beams  is  40  MHz,  then  a  stationary 
particle  in  the  measuring  volume  will  yield  a  40-MHz  signal  from 
the  photodetector.  A  particle  travelling  in  the  same  direction 
as  the  fringes  will  produce  a  lower  frequency  while  a  particle 
moving  in  the  opposite  direction  will  raise  the  frequency. 
Frequency  shifting  provides  a  zero  offset  or  bias  allowing  bi¬ 
directional  measurements.  Techniques  for  frequency  shifting 
include  Bragg  cells  (40  MHz-shift)  and  rotating  diffraction 
gratings  (a  few  MHz  shift). 


The  region  of  the  beam  crossover  point  called  the  measuring 
volume  may  be  considered  shaped  as  shown  in  figure  42.  The  di¬ 
ameter  of  the  measuring  volume,  d  ,  and  the  length  of  measuring 

m 

volume,  1  ,  are  expressed  in  terms  of  the  diameter  of  the  laser 
in 

after  focusing,  d^,  and  the  beam  intersection  angle  0  as 
follows: 


A  Air 

l  “  IT  D  -2  » 
e 

d  -  d* 
m  sin  0/2  ’ 


1 

m 


d 


i 


cos  0/2 


(45) 


(46) 


(47) 
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The  volume  of  the  measuring  volume  is  given  by  equation  (48) 

-2  3 

V  -  ir(D  )  / 6  cos(0/2)  sin(0/2).  (48) 

m  e 

_2 

The  laser  beam  diameter,  Dg  ,  is  defined  as  that  region  where  the 

light  intensity  is  greater  than  1/  2  (about  14%)  of  the  center- 

e 

line  beam  intensity.  Table  5  shows  typical  values  of  these  terms. 

TABLE  5 

MEASURING  VOLUME  PARAMETERS 


Parameter 

Values 

F  (mm) 

103 

241 

598 

0  (°) 

25.0 

11.58 

4.71 

D  -2  (mm) 
e 

.8 

.8 

.8 

d^  (mm) 

.104 

.243 

.602 

d  (mm) 
m 

.106 

.244 

.603 

1  (mm) 

.479 

2.406 

14.657 

d^  (mm) 

1.462  x  10“3 

3.136  x  10~3 

7.700  x  10~3 

nfr 

73.0 

77.0 

78.0 

V  (cm3) 
m 

2.766  x  10"6 

74.600  x  10“6 

2,785.900  x  10"6 

The  position  of  the  measuring  volume,  when  making  measure¬ 
ments  in  a  liquid  flow,  will  be  different  as  compared  to  that  in 
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air  flow.  Figure  43  shows  the  location  of  the  measuring  volume 
relative  to  the  wall  in  terms  of  the  parameters  of  Snell's  law 
equation,  such  that 


Xlnlm  X2n2  * 


(49) 


The  position  of  the  measuring  volume  may  be  determined  with 
the  use  of  equations  (50)  and  (51)  as  follows: 

01  *2 

r-jSin  —  -  n2sin  —  ,  (50) 

h 

*  m  2  cot  ~2  '  (51) 


Photodetector 


The  quantum  efficiency,  a  measure  of  the  photodetector's 

efficiency  to  convert  light  energy  to  electrical  current,  is 

typically  higher  for  photodiodes  than  for  photomultiplier  tubes. 

However,  photomultipler  tubes  make  up  for  this  lack  of  efficiency 

4  7 

with  internal  gains  of  from  10  to  10  .  In  general,  the  selection 
of  a  photodetector  depends  on  signal  intensity,  with  low  scattered 
light  levels  requiring  photomultiplier  tubes;  high  light  levels 
require  photodetectors.  The  photomultiplier  also  has  higher 
frequency  response  than  photodiodes  making  them  more  applicable 
to  frequency  shifted  systems. 


The  signal  processor  used  is  of  the  tracker  type.  Basically 
it  consists  of  a  combination  of  tracking  filter  and  frequency-to- 


voltage  converter,  with  provision  for  signal  validation  and  re¬ 
tention  of  the  last  valid  reading  when  no  signal  is  present.  The 
fundamental  limitation  of  a  tracker  is  the  capture  range.  Should 
the  change  in  velocity  between  one  particle  and  the  next  exceed 
the  capture  range,  the  unit  will  not  track  the  new  value  and  any 
succeeding  values  that  remain  outside  the  capture  range.  However, 
the  tracker  system  will  return  the  last  valid  reading  until  a  new 
reading  occurs  within  the  captive  bandwidth.  Table  6  indicates 
the  range  of  tracking  frequencies  available  with  this  system. 


TABLE  6 

TRACKING  FREQUENCIES 


Frequency 

Tracking 

Range 

Range 

Output 

Capture 

(MHz) 

(FS-10V  DC) 

Window 

High  frequency 

0.5  -  50 

0.2V/MHZ 

5  MHZ 

Mid  frequency 

0.02  -  5 

2.0V/MHZ 

0.5  MHz 

Low  frequency 

0.002  -  0.5 

20V/MHZ 

0.05  MHz 
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Particles 


Particles  in  the  flow  form  the  basic  requirement  for  the 
application  of  LDA  techniques.  The  relative  relationships  between 
mode  of  operation  and  direction  of  maximum  scattered  light  in¬ 
tensity  have  been  discussed  previously.  The  types  of  particles 
and  their  concentration  are  also  related  to  the  mode  of  LDA 
operation  as  well  as  the  fluid  medium  to  be  measured. 

For  gas  flows,  the  diameter  of  particles  should  be  of  the 
order  of  0.1  to  1  micron  while  liquid  flows  generally  required  1 
to  10  micron  particles.  Excellent  results  have  also  been  attained 
with  natural  contaminants  in  water  flows.  If  possible  the  water 
used  should  be  filtered  to  remove  large  particles  providing  a 
uniform  dispersal  of  contaminants  and  lessening  the  chance  of 
large  particles  saturating  the  photodetector  output. 

Mailing  (1971)  has  shown  that  1-micron  sized  particles  are 
able  to  follow  flow  fluctuation  of  up  to  10  kHz  in  air  with  an 
accuracy  of  1.0Z.  Table  7  presents  a  list  of  types  of  seeding 
and  the  medium  in  which  they  are  used. 
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TABLE  7 

TYPES  OF  SEEDING 


Medium 

Seed 

Water 

Polystyrene  beads 

Milk,  cream  powder 

Natural  pigment 

Hydraulic  oil 

Paint  pigment 

Air 

Ammonium  chloride 

Magnesium  oxide 
Silicon  oil  droplets 
Smoke 

Water  vapor 
Dioctyl  phthalate 
Aluminum  oxide 


Burning  gas  Silicon  dioxide 

Magnesium  oxide 
Silicon  oil  droplets 
Dibutyl  phthalate 


116 


Laser  Power 


A  determination  of  laser  power  required  to  attain  a  desired 
signal-to-noise  ratio  (SNR)  requires  precise  information  on  the 
configuration,  characteristics  of  the  fluid,  properties,  and  con¬ 
centration  of  the  scattering  centers  as  well  as  characteristics 
of  the  photodetector.  This  selection  is  quite  complex  and  beyond 
the  scope  of  this  work.  However,  the  following  general  remarks  by 
Fingerson  (1976)  may  be  used  as  a  general  laser  selection  guide. 

In  forward  scatter  with  focal  lengths  up  to  250  mm  for  sub¬ 
sonic  flows  containing  scattering  centers  above  0.1  micron  in 
size,  a  5-mW  laser  is  adequate  and  generally  of  the  He-Ne  type. 

In  the  backscatter  mode  with  short  focal  distances  of  120  mm 
or  less  and  favorable  conditions,  a  15-mV  laser  is  the  minimum. 

For  long  focal  lengths,  especially  in  backscatter,  1-W  lasers  of 
the  Argon-Ion  type  are  generally  used.  For  these  applications,  a 
variable  power  adjustment  feature  is  desirable  in  matching  laser 
power  to  actual  conditions. 

In  general  a  5-mW  or  15-mW  He-Ne  laser  is  adequate  for  most 
laboratory  experiments  particularly  where  forward  scatter  can  be 
used.  For  increased  focal  length,  increased  velocity,  backscatter 
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or  dual  component  velocity  measurements,  a  larger  laser  such  as 
Argon-Ion  is  desirable. 

An  oscilloscope  photograph  of  the  photodiode  output  due  to 
the  reflected  light  from  a  single  particle  travelling  through  the 
measuring  volume  is  shown  in  figure  4.  The  total  Doppler  burst 
consists  of  approximately  100  fringes  indicating  a  particle 
travelling  close  to  the  center  of  the  measuring  volume.  Figure  45 
represents  good  particle  distribution,  providing  nearly  continuous 
signal  output  due  to  sequential  occupancy  of  the  measuring  volume 
by  individual  particles.  Multiple  particles  in  the  measuring 
volume  at  the  same  time  are  to  be  avoided.  This  type  of  output 
is  shown  in  figure  46.  Although  providing  nearly  continuous 
signal,  the  relatively  high  background  noise  level,  caused  by 
varying  phase  and  amplitude,  affects  the  accuracy  of  the  apparent 
frequency  measurement.  The  accepted  term  for  this  is  Doppler 

3 

ambiguity.  A  count  rate  distribution  of  2.7  x  10  particles/ 
second  corresponding  to  figure  45  is  shown  in  figure  47.  Each 
spike  contains  a  Doppler  burst  envelope  for  a  single  particle. 

Errors  in  Laser  Doppler  Average 

Edwards  et  al.  (1971),  and  George  and  Lumley  (1973)  treat 
laser  anemometry  measurement  errors  in  detail.  Limitations  of  a 
theoretical  nature  include: 


▲ 


Figure  47.  Count  rate  distribution  photograph 
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1.  Velocity  fluctuations  in  the  measuring  volume 

2.  Mean  velocity  gradients  across  the  volume 

3.  Effect  of  the  finite  transit  time  of  the  particles 
through  the  volume. 

The  finite  measuring  volume  of  the  laser  anemometry  has  an 
effect  on  mean  flow  measurements  in  steep  velocity  gradients. 
Turbulent  particle  movement  within  the  measuring  volume  can  also 
influence  resulting  turbulence  data.  Since  the  Doppler  signal 
from  a  single  particle  has  a  beginning  and  an  end,  higher  fre¬ 
quency  variations  are  generated  appearing  as  spectral  broadening 
(even  in  laminar  flow) .  This  spectral  broadening  gives  a  min¬ 
imum  turbulence  intensity  which  can  be  measured. 

Doppler  ambiguity  caused  by  phase  and  amplitude  modulation  of 
the  Doppler  signal  is  a  source  of  inherent  noise.  This  noise  is 
not  a  problem  at  low  particle  concentrations  but  becomes  signifi¬ 
cant  where  the  probability  of  multiple  particle  in  the  measuring 
volume  is  high.  However,  interpretation  of  spectral  broadening 
and  its  effects  on  signal  processor  remains  controversial. 

Even  under  the  most  ideal  conditions,  there  are  certain 
errors  associated  with  laser  Doppler  anemometry.  These  errors 
are  frequently  called  "broadening  errors"  since,  in  the  frequency 
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domain,  Chey  would  cause  Che  frequency  to  b.  broadened  into  a 
distribution  with  a  finite  width.  In  a  tracker,  this  broadening 
translates  into  a  voltage  fluctuation  in  the  output.  The  major 
sources  of  broadening  are 


1.  Oj,  -  finite  transit  time  broadening 

2.  -  small  scale  turbulence  fluctuations 

3.  aG  -  mean  velocity  gradient  broadening 

4.  Oj  -  instrument  broadening 


If  it  i3  assumed  that  each  of  the  broadening  components  has  a 
Gaussian  distribution,  the  output  variance  may  be  written  as 


(  2  )  *  (~~2  )+(  2  )  "*”  (  2/+(  2/' 

\Z  h  \E  4  \E  4  \E JG  \E  4 


(52) 


If  a  measurement  in  a  laminar  flow  with  no  velocity  gradient 
is  considered,  then  the  broadening  reduces  to 


(53) 


The  second  term  on  the  right  hand  side,  the  instrument  broadening, 
may  be  determined  by  entering  a  constant  amplitude,  constant  fre¬ 
quency,  voltage  fluctuations  as  the  tracker  input,  and  observing 
the  corresponding  voltage  fluctuation  at  the  output.  This  gives 
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the  level  of  the  Instrument  broadening.  Comparison  with  the 
results  of  the  laminar  flow  measurement  allow  the  finite  transit 
time  broadening  to  be  determined. 

The  transit  time  broadening,  a  ,  comes  about  as  a  result  of 

r 

the  signal  from  any  one  particle  existing  for  only  a  finite  time. 
If  there  is  more  than  one  particle  in  the  control  volume  at  a 
given  time,  then  there  will  be  voltage  and  phase  variations  due  to 
"old"  particles  leaving  the  control  volume  and  "new"  particles 
entering  it.  The  finite  transit  time  broadening  is  a  strong 
function  of  control  volume  geometry  and  many  discussions  of  its 
cause  and  effects  may  be  found  in  the  bibliography. 


and  2  scattering  volume  dimension  in  the  direction  of  the 
velocity  gradient,  and  =  Doppler  frequency. 

The  previous  approximation  shows  the  desirability  of  keeping 
the  smallest  dimension  of  the  control  volume  in  the  direction 
of  the  velocity  gradient  in  order  to  minimize  gradient  effects. 

The  final  source  of  broadening  is  usually  the  parameter  of 
interest.  In  some  measurements  the  turbulence  variation  has  been 
taken  as  the  difference  between  the  variance  measured  in  turbulent 
flow  and  that  measured  in  laminar  flow  under  the  same  conditions . 
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Data  Processing 

A  modes t  real- else  computer  data  system  was  used  la.  con¬ 
junction  with  the  laser  Doppler  anemometer.  The  data  system 
components  shown  in  figure  43  are  listed  in  table  3. 

TART. 7  3 

DATA  SYSTEM  COMPONENTS 


Laser  Doppler  Anemometer 

Display  Keyboard 

Central  Processing  Unit  (CPU) 

Plotter 

Printer 

Data  Interface _ 


TSI  Model  1090-1 
Wang  Model  2220 
Wang  Model  2200-72 
Wang  Model  2272-2 
Wang  Model  223 IV 
NUSC _ 


The  data  system  provides  a  printout  of  seven  variables  and 
functions  of  variables,  a  plot  of  normalized  velocity  ratio  versus 
position,  and  a  cassette  recording  of  four  data  functions.  Table 
9  lists  printer  output. 
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TABLE  9 

PRINTER  OUTPUT  LIST 


Variables 

Variable  Functions 

y  (inches)* 

Height  above  plate 

u  (ft/sec) 

Velocity 

1  * 
u  (ft/sec) 

Velocity  fluctuation 

u/  * 

'UFS 

Normalized  velocity  ratio 

u^/u 

Local  turbulence  level 

u1/UFS 

Relative  turbulence  level 

volts 

Tracker  output 

*  stored  on  tape  cassette. 

The  program  "DATA"  selects  61  data  points  from  each  channel 
of  input  and  computes  their  average  for  computational  use.  The 
four  channels  of  Input  include  velocity  and  rms  of  the  velocity 
fluctuation  from  the  LDA  system  output,  position  y  from  the 
digital  position  indicator,  and  a  counter  channel. 

From  listings  and  laboratory  data,  x  station  and  velocity 
profile  data  were  key  punched  as  input  to  a  comprehensive  data 
processing  program  called  "DATAPAC."  This  program  calculated 
hydrodynamic  parameters,  listed  data,  plotted  data,  and  theoret¬ 
ical  curves  using  a  high-speed  Calcomp  plotter.  Appendix  B 
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presents  a  copy  of  program  DAXAPAC  and  a  sample  output  listing. 
Typical  output  curves  from  one  velocity  profile  are  shown  in 
appendix  C. 


IV.  EXPERIMENTAL  PROCEDURES 

Three  types  of  tests  were  performed  in  this  research.  These 
were:  (1)  tunnel  characterization  tests,  (2)  water  injection 

boundary  layer  sampling  tests,  and  (3)  polymer  injection  boundary 
layer  sampling  tests.  The  preparations  for  these  experiments  will 
be  discussed  first. 


Polymer  Preparation 

The  polymer  used  in  this  experiment  was  polyethylene  oxide 
(Polyox  WSR-301)  of  approximate  molecular  weight  6  x  106.  This 
polymer  was  selected  due  to  its  availability  and  widespread  use 
among  researchers,  thereby  minimizing  data  correlation  problems 

due  to  use  of  polymers  of  different  molecular  weight. 

As  the  experimental  program  called  for  solutions  of  1000  WPPM 
or  less  to  be  tested,  a  master  solution  was  prepared  at  2000  WPPM 
and  diluted  to  the  required  concentrations. 

The  master  solution  was  prepared  using  the  same  procedure  as 
Sirmalis  (1974).  Premeasured  Polyox  powder  was  sifted  onto  the 
surface  of  the  carefully  weighed  water,  which  was  slowly  stirred  by 
a  magnetic  mixing  bar.  After  a  clear  solution  resulted  (1  hour), 


130 


131 


the  solution  was  left  to  stand  for  approximately  40  hours  to 
minimize  the  high  viscoelastic  effects  of  freshly  mixed  solution 
and  to  assure  homogeneity.  The  master  solution  was  kept  in  a 
dark,  cool  place  after  initial  mixing  to  reduce  the  auto-oxidation 
problem  with  polyethers.  When  required,  the  proper  amount  was 
pipetted  out  of  the  master  solution  and  diluted  to  the  required 
concentration.  These  dilute  solutions  were  used  within  24  hours 
of  the  makeup  time. 

The  Uranlne  master  solution  was  prepared  in  a  similar  manner. 
When  a  test  solution  was  required,  the  proper  amount  of  each 
constituent  was  added  to  enough  water  to  make  one  liter  of 
solution.  This  solution  was  poured  back  and  forth  into  the  con¬ 
tainer  ten  times  to  assure  complete  mixing  and  then  let  to  stand 
for  1  hour  to  achieve  homogeneity.  A  series  of  calibration  tests 
with  the  fluorometer  indicated  that  this  procedure  provides 
accurate  results. 


Injection  System  Calibration 


The  polymer  injection  system  was  designed  to  minimize 


mechanical  degradation  of  the  polymer  concentration  and  avoid 
polymer  effects  on  the  flow  measuring  device.  This  was  accomp¬ 
lished  by  using  water,  bypassed  from  a  recirculatory  loop  to 
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displace  polymer  from  a  reservoir.  The  reservoir  contained  a 
free-sliding  piston  with  0-rlng  seals' and  was  thoroughly  cleaned 
after  each  test  to  prevent  alteration  of  the  next  concentration 
to  be  tested.  The  water  flow  passed  through  a  Fisher  Scientific 
Company  variable  area  flowmeter  with  tube  number  448-225.  The 
float  diameter  was  1/8  Inch  and  made  of  stainless  steel. 

3 

The  calibration  curve  for  the  flow  range  of  20,  40,  and  80  cm 
per  minute  is  shown  in  figure  49.  The  accuracy  of  the  curve 
applied  to  the  delivery  system  was  confirmed  by  the  capture  of 
polymer  solution  in  a  graduated  cylinder.  There  were  no  oscil¬ 
lations  of  the  float  in  the  variable  area  tube  indicating  the 
steadiness  of  delivery. 

Boundary  Layer  Sampling 

Five  boundary  layer  sampling  stations  were  used  in  these 
experiments,  located  at  axial  distances  of  3.75,  8.5,  12.5,  16.5, 
and  20.5  inches  from  the  leading  edge  of  the  plate.  All  probe 
rakes  were  retracted  up  into  the  ceiling  of  the  tunnel  out  of 
the  flow  at  the  start  of  the  test.  The  rake  furthest  downstream 
was  lowered  to  the  plate  first  and  a  set  of  samples  captured  in  / 

the  cuvettes.  The  wall  sample  was  not  taken  through  the  hole 
in  the  plate  at  this  time  to  prevent  probe  interferences.  When 
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Figure  49.  Polymer  injection  flowmeter  calibration  curve 
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a  sufficient  sample  volume  was  captured,  at  approximately  3  to 
10  minutes,  the  sample  lines  were  transferred  to  a  reservoir 
and  the  probe  rake  raised  to  position  number  2.  The  new  position 
was  determined  by  a  gage  block  inserted  between  a  collar  on  the 
rake  body  and  the  top  of  the  tunnel  fitting.  Using  the  reservoir 
and  placing  the  probes  in  the  next  position  allowed  the  lines 
to  purge  themselves  while  the  next  set  of  cuvettes  was  being 
readied.  When  the  probe  was  raised  to  position  number  2,  the 
wall  sample  was  taken.  After  all  the  samples  were  taken  at  the 
station,  the  rake  was  retracted  up  into  the  ceiling  of  the  tunnel, 
out  of  the  flow.  The  upstream  probe  stations  were  used,  in  order, 
until  all  the  samples  were  taken.  All  samples  were  expelled 
under  internal  tunnel  pressure.  No  suction  was  used  to  avoid 
interference  effects  between  probes  during  sampling.  Volume 
rates  of  sampling  were  based  on  the  slowest  velocity  in  the 
boundary  layer  that  allowed  all  the  samples  to  be  taken  within 
the  capability  of  the  delivery  system.  Table  10  presents  the 
non-dimensional  locations  of  each  of  the  stations  at  which  con¬ 
centration  measurements  were  taken. 

Fluorometer  Calibration 


The  fluorometer  used  in  these  tests  is  a  G.K.  Turner  Model 
111  Fluorometer.  Uranine-B  dye  was  selected  since  it  was  readily 
adaptable  to  laboratory  use,  compatible  with  polymer  flows,  and 


•  M 


9 

gave  reproducible  readings  Co  abouc  1  part  in  10  on  a  weight 
basis.  This  range  was  more  than  sufficient  for  the  experiment  to 


TABLE  10 

NON-DIMENSIONAL  LOCATIONS  OF  CONCENTRATION  MEASUREMENTS  STATIONS 


Station 

X 

(Inches) 

X/L 

(L*22.5  Inches) 

X/S 

(S».020  Inches) 

X/S 

(S-.005  Inches) 

0 

1.75 

.078 

87.5 

350.0 

1 

3.75 

.167 

187.5 

750.0 

2 

8.50 

.378 

425.0 

1700.0 

3 

12.50 

.556 

625.0 

2500.0 

4 

16.50 

.773 

825.0 

3300.0 

5 

20.50 

.911 

1025.0 

4100.0 

be  performed  as  the  dye  concentration  in  the  solution  to  be 
ejected  could  be  adjusted,  thus  producing  samples  along  tne  plate 
sufficiently  above  background  levels  to  give  reproducible  results. 
The  Model  111  Fluorometer  has  four  sensitivity  settings  allowing 
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calibration  curves  to  be  drawn  over  a  range  of  about  1  part  per 
6  9 

10  to  1  part  per  10  ,  on  a  weight  basis.  For  concentrations 
greater  than  the  lower  value,  dilution  prior  to  measurement  was 
performed.  The  calibration  charts  for  this  experiment  are  shown 
in  figures  33  and  34.  Approximately  3.5  to  4.5  ml  of  sample  is 
required  for  measurement.  The  sample,  contained  in  a  cuvette 
(a  small  test  tube),  is  placed  in  the  instrument  and  the  measure¬ 
ment  is  taken.  Cleanliness  was  found  to  be  a  key  factor  in  the 
technique  as  fingerprints,  body  oils,  and  trace  contamination 
could  affect  readings.  All  cuvettes  were  carefully  wiped  with 
soft  laboratory  paper  towels  and  handled  only  by  the  upper  edge. 
All  samples  were  allowed  to  achieve  room  temperature  prior  to 
reading. 


Wall  Location 


The  location  of  the  plate  surface  was  determined  with  use  of 
the  laser  beam.  The  beam  of  diameter  d  *  .0096  inches  was 

SI 

lowered  until  two  visual  effects  were  observed.  First,  surface 
imperfections  were  observed,  highlighted  by  the  laser  beam  as  it 
approached  the  plate.  Secondly,  with  the  aid  of  a  viewing  tele¬ 
scope,  observation  of  the  measuring  volume  at  the  beam  crossover 
point  showed  the  reflection  of  the  intersected  beams  from  ths 
plate  surface  as  a  second  image.  When  the  second  image  blended 
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with  the  original  image,  the  beams  were  considered  on  the  surface. 
The  center  of  the  measuring  volume  was  then  about  0.0048  inch  off 
the  surface  of  the  plate.  The  beams  were  lowered  to  that  dis¬ 
tance  and  the  vertical  height  digital  indicator  set  at  zero.  The 
measurements  were  all  considered  to  be  made  at  a  height  correspond¬ 
ing  to  the  center  of  the  measuring  volume  above  the  surface  of 
the  place. 


Tunnel  Characterization 


To  determine  the  suitability  of  the  water  tunnel  for  the 
experiment,  a  velocity  field  map  of  the  tunnel  was  made  using  the 
LDA  system.  (The  two-dimensional  character  of  the  tunnel  is 
shown  in  figure  50.)  Additional  tests  downstream  and  at  higher 
flow  rates  also  showed  that  corner  flow  effects,  although  present, 
were  removed  from  the  center  of  the  tunnel. 

Velocity  profiles  and  turbulence  intensity  distributions  were 
taken  and  are  presented  in  comparison  with  the  work  of  Laufer 
(1951)  as  figures  51  and  52.  These  curves  are  in  excellent 
agreement . 

A  critical  requirement  for  this  experiment  was  that  transi¬ 
tion  be  effected  on  the  plate,  preferably  at  its  mid-length  to 
allow  observation  of  laminar,  transitional,  and  turbulent  flow 
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boundary  layers.  The  incoming  velocity  profile  is  shown  in  figure 
53.  Velocity  profiles  at  axial  distances  of  3.75,  14,  and  18 
inches,  presented  in  figure  58,  are  also  laminar,  transitional, 
and  turbulent.  Ejection  of  a  dye  onto  the  surface  of  the  plate 
showed  turbulent  bursting  occurring  at  X  *  6.00  inches. 

Application  of  a  laser  Doppler  technique  to  locate  transition 
on  a  flat  plate  as  developed  by  P.  E.  Gibson  of  the  Naval 
Underwater  Systems  Center,  Newport,  Rhode  Island,  is  shown  in 
figure  54.  From  laminar  and  turbulent  velocity  profiles,  the 
height  above  the  plate  which  yielded  the  greatest  corresponding 
velocity  difference  was  determined.  The  laser  measuring  volume 
was  then  set  at  that  height  and  transversed  the  length  of  the 
tunnel,  producing  a  plot  of  axial  velocity  distribution.  This 
technique  confirmed  that  transition  was  occurring  on  the  plate. 

Test  Series 


The  test  program  consisted  of  two  parts:  (1)  water  ejection 

tests,  and  (2)  polymer  ejection  tests.  All  tests  were  run  at  a 

constant  tunnel  flow  of  8  gal/min  and  an  average  velocity,  V&  , 

of  2.2  ft/sec.  Volumetric  flow  rates,  Q^,  of  20,  40,  and  80 
3 

cm  /min  were  selected  to  be  injected  through  two  slot  openings, 


s,  of  0.020  and  0.005  inch.  Table  11  lists  the  injection  parameter 
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used.  At  a  given  velocity  of  ejection  determined  by  the  combi- 
[  nation  of  slot  area  and  volumetric  ejection  flow  rate,  various 

|  concentrations  of  polymer  were  ejected.  Table  12  presents  the 

[  concentration  flux  matrix  tested. 
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V.  ANALYTICAL  CONSIDERATIONS 


Boundary  Laver  Model 

This  section  presents  an  analytical  approach  for  developing 
turbulent  flow  over  a  submerged  flat  plate  with  polymer  injection. 

A  model  is  developed  for  predicting  growth  of  the  boundary  layer, 
diffusion  of  polymers,  and  prediction  of  skin  friction  coefficients 
and  transition  as  a  function  of  local  wall  polymer  concentrations. 

In  the  analysis  that  follows,  the  velocity  profile,  continuity 
equation,  momentum  equation,  and  an  expression  for  polymer  wall 
concentration  from  a  regression  analysis  which,  in  turn,  was 
obtained  from  experimental  data  are  combined  with  an  effective 
Reynolds  number  analogy.  This  combination  yields  expressions  for 
the  growth  of  the  boundary  layer  and  the  skin  friction  coefficient 
as  a  function  of  the  local  polymer  wall  concentration. 

Velocity  Profile  Relation 

Meyer  (1966)  provided  a  AB  correction  for  the  law  of  the  wall 
accounting  for  polymer  effects 

u+  -  7  In  y+  +  B  +  AB.  (56) 
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An  alternate  method  of  accounting  for  the  polymer  effect  would 

be  by  adjustment  of  the  mixing  length  constant,  according  to  Virk 

(1971) .  In  his  work  on  maximum  drag  reduction  along  the  ultimate 

asymptote  line,  he  noted  that  a  change  in  K  by  a  factor  of  up  to 

5  Is  possible  and  In  full  agreement  with  data.  However,  from 

his  studies,  Meyer  concluded  that  drag  reduction  in  pipes  is  due 

to  thickening  of  the  laminar  sublayer.  The  mixing  length  constant, 

K,  which  equals  approximately  0.4,  is  apparently  unaffected  by 

dilute  concentrations  of  additive.  Meyer  further  correlated  the 

constant  B  and  showed  it  to  remain  constant  at  the  Newtonian  value 

* 

of  5.5  until  the  friction  velocity  reached  a  threshold  value,  Vq  , 

* 

after  which  B  increased  logarithmically  with  V  , 

B  *  Z  ln(V*/V  *).  (57) 

o 

The  term  Z  is  a  dimensionless  constant  dependent  on  the  type  of 
polymer  additive  and  the  concentration.  A  constitutive  relation 
for  Z  is  suggested  by  White  (1968)  from  simple  curve  fit  expres¬ 
sions  as 


Z 


««v>\ 


(58) 


where  a  *  2.3  and  y  *  0.5. 


(59) 
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Skin  Friction  Relations 


The  boundary  layer  continuity,  and  momentum  equations  for  flow 
over  a  submerged  flat  plate  with  zero  pressure  gradient  as  applied 
to  the  polymer  flow  case  studied  are  given  as 


3(pu)  ,  3(pv)  ,  n 

3x  3y  u 


(60) 


and 


pu 


3r 

3y‘ 


(61) 


Upon  changing  independent  variables  x  and  y  to  x,  y+,  and 
Cw,  the  x  derivatives  must  be  handled  by  the  chain  rule,  since 
the  parameters  y+  and  Cw  are  functions  of  x  in  the  law  of  the  wall. 
Thus,  the  substitute  is 


(62) 


As  may  be  noted  in  equation  (62) ,  concentration  derivatives  are 
added.  The  fluid  properties  of  density,  p,  and  viscosity,  u, 
are  assumed  constant  and  equal  to  those  of  the  solvent  water  in 


this  analysis. 
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The  resultant  boundary  layer  equation,  after  considerable 
algebraic  manipulation,  is 


Appendix  A  presents  a  complete  development  of  the  equations 
related  to  this  method. 

With  elimination  of  the  polymer  terms,  equation  (63)  compares 
favorably  with  the  classical  skin  friction  relations  of  Prandtl- 
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Schllchtlng  and  S chul tz-Grunow ,  as  well  as  Che  simplified  relation 
of  F.  White.  This  comparison  is  presented  in  figure  55  •. 

Prior  to  the  solution  of  equation  (63)  and  (65)  ,  it  is  neces¬ 
sary  to  determine  a  means  of  calculating  C  •  A  regression  analysis 

w 

was  performed  on  all  of  the  data  for  the  slowest  injection  rate  of 
20  cc/m±n  and  injection  concentrations  of  c  *  100,  200,  400,  500, 
and  800  WPPM.  Equation  (66)  represents  the  results  of  that 
analysis  which  successfully  correlates  88%  of  the  data  obtained, 
so  that 


^  -  .834  3  +  2. 30  3  3  0**)  -  7.385(|)2  +  4.2855(^)3.  (66) 

Equation  (66)  is  valid  for  X/L  >  .27.  For  X/L  <  .27,  equation 
(67)  holds  true.  It  follows  then  that 

C 

-  1.0.  (67) 

Gi 


Effective  Reynolds  Number  Analogy 

White  (1968)  rearranged  the  friction  factor  relation  derived 
by  Meyer  (1966)  and  evolved  the  following  effective  Reynolds 
number  analogy.  The  friction  factor  in  the  presence  of  polymer 
may  be  considered  equal  to  the  Xewtonian  friction  factor  evaluated 
at  an  effective  Reynolds  number  given  by 
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Thus,  the  polymer  has  the  effect  of  increasing  the  Newtonian 
Reynolds  number.  Re,  and  thereby  yielding  the  lower  polymer  solu¬ 
tion  friction  factor.  For  constant  freestream  velocity  and 
assuming  constant  polymer  ocean  flow,  skin  friction  versus 
Reynolds  number  curves  may  be  obtained  as  follows : 


1.  Calculate  Reynolds  number  at  given  X, 


Re  -  ■=*  .  (69) 

xv 

2.  Calculate  initial  skin  friction  value  using  White's 
simplified  relation. 


.455 


[ In ( . 06  Rex)]‘ 


3.  Calculate  initial  shear  velocity, 


(70) 


(71) 


4.  Calculate  effective  polymer  Reynolds  number, 


(72) 


* 


M 
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where 

2  -  aCY  for  C  <  C*  (73) 

w  w  —  w  v  ' 

2  -  <x(C*)Y  for  C  >  C* 
w  w  w 


5.  Calculate  polymer  skin  friction, 


.455 

[ln(. 06  Re  )]2 
P 


6.  Calculate  shear  velocity. 


(74) 


(75) 


Equations  (72)  through  (75)  should  be  repeated  until  the  shear 
velocity  does  not  change  by  more  than  0.5%.  This  iterative  pro¬ 
cedure  yields  skin  friction  values  as  a  function  of  Reynolds 
numbers.  A  new  value  of  x  along  the  plate  should  be  selected  and 
the  procedure  should  be  repeated.  Figure  56  presents  curves  of 
skin  friction  versus  Reynolds  number  for  concentrations  of  5,  10, 
15,  20, and  25  WPPM  as  well  as  turbulent  and  laminar  flow  skin 
friction  curves  given  by  equations  (70)  and  (76)  for  water  flow 
only, 


(76) 


i 
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Using  the  concept  of  the  effective  Reynolds  number  analogy  and 
the  variation  in  the  wall  concentration  of  polymer  determined 
from  the  regression  analysis  performed  on  the  experimental  data, 
it  is  possible  to  numerically  integrate  the  right  hand  side  of 
equation  (63)  and  compare  it  with  the  resulting  Reynolds  number 
obtained  from  integration  of  the  left  hand  side  of  equation  (63). 

The  computational  procedure  is  presented  in  appendix  F  as 
a  listing  of  computer  program  SKINFRIC.  The  model  allows  for 
an  initial'  laminar  flow  over  the  plate,  changing  to  turbulent 
flow  at  x  ■  x(.  The  integrations  are  performed  from  initial 
values  of  x  where  injection  occurs,  an  initial  wall  concentration 
equal  to  the  injected  concentration  of  polymer  and  an  initial 
value  of  X  determined  from  laminar  flow  conditions  given  by 
equation  (76).  This  is  based  on  the  assumption  that  small  concen¬ 
trations  of  polymers  do  not  affect  laminar  flow.  The  work  of 
Hoyt  and  Fabula  (1964)  showed  that  non-Newtonian  is  a  misnomer 
applied  to  low  concentrations  of  polymers.  They  found  that  these 
solutions  were  of  constant  viscosity  and  greater  magnitude  than 
that  of  the  solvent.  White  (1968)  determined  the  limiting  value 
of  2  to  be  approximately  11.0  by  single  curve  fitting  of  data. 
Application  of  the  maximum  drag  reduction  concentration  for  Polyox 
of  C*  *  30  WPPM  as  determined  by  Hoyt  and  Fabula  (1964)  ,  and 
shown  in  figure  57,  allowed  the  limiting  value  of  2  to  be  taken 


as  12.6. 
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Figura  57.  Percent  drag  reduction  as  a  function 
of  polymer  concentration 


It  is  further  assumed  that  the  effective  Reynolds  number 
analogy  is  a  solution  to  equation  (63) .  Values  of  X  are  then 
used  in  the  numerical  integration  of  the  right  hand  side  of  equa¬ 
tion  (63).  Values  of  the  polymer  vail  concentration,  equation 
(66)  ,  from  the  regression  analysis  are  also  used.  The  integration 
proceeds  from  XQ  ■  0  at  the  point  of  injection  of  the  polymer 
and  neglects  leading  edge  effects.  Laminar  flow  equations  give 
the  flow  over  the  initial  region  of  flat  plate  until  the  preset 
transition  point  is  reached  at  x  ■  after  which  turbulent 
effective  Reynolds  analogy  parameters  apply.  Comparison  of  the 
difference  between  integrand  1  and  integrand  2  on  the  right 
hand  side  of  equation  (63)  with  the  left  hand  side  Reynolds  number 
is  made  at  the  end  of  the  plate.  If  they  are  not  within  1 Z,  the 
location  of  the  point  of  transition  is  adjusted  and  the  procedure 
repeated  until  the  results  are  within  the  set  tolerance  limits. 

The  transition  region  is  defined  as  starting  at  the  transition 
point  and  ending  at  the  distance  where  the  Reynolds  number  equals 
350,000.  In  this  zone  a  linear  intermittency  factor  is  assumed 
and  an  appropriately  proportioned  friction  factor  calculated. 

Typical  plots  of  skin  friction  coefficient  versus  Reynolds 
number  for  various  polymer  concentrations  injected  at  »  20  cc/min 
are  shown  in  figures  92  through  96.  Boundary  layer  growth  over 
the  length  of  the  plate  is  presented  in  figures  85  through  88  and 
compared  with  experimental  data. 


VI.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


General  Remarks 


The  experimental  program  was  conducted  at  a  constant  tunnel 
flow  of  8  gallons  per  minute.  Injection  Into  the  boundary 
layer  was  performed  through  two  different  size  slots  of  width, 
s  -  .020  Inch,  and  s  *  .005  Inch,  at  Injection  rates  of  20,  40, 
and  80  cc/mln.  These  six  injection  velocity  ratios  were  tested 
with  dyed  water  injection.  Polymer  concentrations  of  100,  200, 

400,  500,  and  800  WPPM  were  Injected  at  20,  40,  and  80  cc/min 
through  the  s  ■  .020  inch  slot  and  at  80  cc/mln.  through  the  s  * 
.005-inch  slot,  producing  injection  velocity  to  average  freestream 
velocity  ratios  of  0.0174,  0.0348,  0.0696  and  0.2795,  respectively. 

Velocity  profiles  were  taken  at  nine  stations  along  the  plate 
corresponding  to  axial  distances  of  x  *  3.75,  6,  8,  10,  14,  15, 

16,  17,  and  18  inches.  Concentration  measurements  were  taken  at 
five  stations  along  the  plate  corresponding  to  axial  distances  of 
x  ■  3.75,  3.5,  12.5,  16.5,  and  20.5  Inches.  At  each  concentration 
measuring  station,  nine  samples  were  taken  throughout  the  boundary 
layer  including  a  wall  sample.  Table  2  indicates  the  relative 
heights  in  the  boundary  layer  at  which  the  concentration  samples 
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were  captured  according  to  which  spacer  was  used  In  conjunction 
with  the  sampling  probe. 

In  total,  243  velocity  profiles  and  130  concentration 
profiles  were  measured.  Space  limitations  prohibit  the  complete 
listing  of  the  data  file,'  however,  copies  of  the  data  file 
are  available  from  the  author  at  the  Naval  Underwater  Systems 
Center,  Newport,  Rhode  Island.  The  focus  of  the  discussion 
will  concentrate  on  the  20  cc/min  injection  rate  through  the  .020- 
inch  slot  which  yielded  an  injection  velocity  ratio  of  0.0174. 

This  slowest  injection  velocity  was  selected  as  providing  minimum 
boundary  layer  disturbance  at  the  point  of  injection.  Previous 
investigations  such  as  Sirmalis  (1976),  Fruman  and  Tulin  (1976), 

Wu  (1968),  and  Wells  (1968)  have  recommended  minimization  of  the 
injection  velocity.  The  hydrodynamic  parameters  for  the  tests 
discussed  are  presented  in  appendix  D  with  velocity  profiles 
presented  in  appendix  E. 
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Water  Flow  Boundary  Layer  Characterization 

Velocity  profiles  for  the  laminar  flow  region  are  presented 
in  figure  58  and  compared  with  the  classical  Blasius  profile. 

The  slight  variation  from  the  classical  form  may  be  attributed 
to  the  effects  of  the  small  pressure  gradient  at  the  leading  edge 
of  the  plate  and  due  to  its  shape  and  roughness.  Displacements  among 
the  data  themselves  are  due  to  the  effects  of  injection  velocity, 
with  minimum  injection  velocity  approaching  the  theoretical  profile. 

Downstream,  at  a  station  x  ■  18  Inches,  the  velocity  profiles 
follow  a  l/6th  power  law  turbulent  flow  shape  as  shown  in  figure 
59.  Patel  (1968)  has  shown  that  in  developing  flows  the  velocity 
profiles  go  through  various  power  law  shaped  profiles  until  a 
l/7th  power  law  profile  is  achieved  for  fully-developed  turbulent 
flow.  The  measured  profiles  are  in  agreement  with  Patel's  findings 
as  fully-developed  turbulent  flow  was  not  totally  achieved  due  to 
the  length  of  the  plate  used. 

The  transitional  nature  of  the  velocity  profiles  from  laminar 
to  turbulent  flow  for  the  case  of  20  cc/min  injection  at  a 
velocity  of  injection,  Vi  ■  .0383  ft/sec,  is  shown  in  figure  60. 
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NORMALIZED  VELOCITY  U/UFS 
Figure  58.  Normalized  laminar  velocity  profiles 
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The  turbulent  Intensity  distributions  for  these  stations  are 
presented  in  figure  61.  As  shown  by  Schubauer  and  Klebanoff  (1955) , 
the  velocity  becomes  greater  closer  to  the  wall  as  the  flow  becomes 
turbulent.  Transition  to  turbulent  flow  is  generally  observed  by 
a  sudden,  rapid  growth  in  the  height  of  the  boundary  layer.  The 
growth  of  the  boundary  layer  for  water  injection  is  shown  in  figure 
62  and  reveals  this  characteristic. 

An  additional  criterion  to  characterize  turbulent  flows  is 
that  the  velocity  distribution  in  the  wall  region  should  follow 
the  well  known  law  of  the  wall,  such  that 

u+  »  A  In  y+  +  B, 

where  A  and  B  are  universal  constants  having  the  value  of  2.5  and 
5.5,  respectively.  Patel  (1968)  indicates  that,  for  developing 
flow,  the  value  of  B  decreases  until  the  value  of  5.5  is  achieved 
for  fully-developed  turbulent  flow.  Application  of  this  concept 
to  the  channel  flows  under  discussion  revealed  the  variation  in 
the  value  of  B  as  shown  in  figure  63  as  a  function  of  length 
Reynolds  number.  A  representative  value  of  B  *  7.2  was  used  in 
this  study  for  the  turbulent  flow  region  which  is  not  fully  developed. 
Figure  64  presents  a  plot  of  the  law  of  the  wall  in  inner  variables 
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Figure  63.  Variation  in  value  of  B  as  a  function  of 
Reynolds  number  for  developing  channel  flow 
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y  and  u  for  B  *  5.5  and  7.2.  A  comparison  Is  also  made  with 
experimental  data  showing  excellent  agreement  with  the  value  of 
B  «  7.2.  The  value  of  B  ■  7.2  is  therefore  used  in  the  expression 
for  the  law  of  the  wall  to  account  for  the  developing  velocity 
profile  aspect  of  the  flow. 

Polymer  Flow  Boundary  Layer  Characterization 

Developing  velocity  profiles  for  injected  polymer  concentra¬ 
tions  of  ■  100,  200,  400,  500,  and  800  WPPM  are  presented  in 
figures  65  through  69.  The  profiles  are  from  stations  x  * 

3.75,  10.0,  and  18.0  inches  representing  laminar,  transitional, 
and  turbulent  flow  regions.  The  characteristic  increasing  fullness 
of  the  velocity  profiles  over  the  length  of  the  plate  is  evident  in 
each  test  sequence. 

Turbulence  intensity  distributions  for  these  profiles  are 
plotted  in  figures  70  through  74.  The  freestream  turbulence 
intensity  was  .3%  for  the  tests  where  water  and  polymer  concen¬ 
trations  equal  to  100,  200,  and  400  WPPM  were  injected.  For  the 
injection  of  polymer  at  500  and  800  WPPM,  the  freestream  turbulence 
level  was  .4"  and  .5 Z,  respectively. 
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Figure  69.  Developing  velocity  profiles  for  injected  polymer 
concentration  »  800  WPPM 


Figure  70.  Turbulence  intensity  distribution  for  injected  polyme 
concentration  C.  •  100  WPPM 
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Figure  73.  Turbulence  intensity  distribution  for  injected  polymer 
concentration  =*  500  WPPM 
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Figure  74.  Turbulence  Intensity  distribution  for  injected  polymer 
concentration  “  800  WPPM 
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The  velocity  profiles  at  %  m  3.75  inches  are  all  laminar  in 
shape.  At  x  ■  10  inches,  the  velocity  profiles  appear  less  turbu¬ 
lent  than  the  water  case.  This  characteristic  is  also  shown  in 
the  lower  turbulence  intensity  distributions  at  this  station  for 
all  concentrations  injected.  The  polymer  causes  a  reduction  in 
the  peak  value  of  turbulence  intensity  and  shifts  the  position 
further  away  from  the  wall.  The  entanglement  or  "hand  holding" 
of  the  long  chain  polymers  may  act  as  a  woven  blanket  to  dampen 
or  suppress  turbulence. 

A  comparison  of  velocity  profiles  at  x  *  18  inches,  normalized 
with  respect  to  boundary  layer  thickness  and  freestream  velocity 
compared  with  a  l/5th  power  law  velocity  profile  in  figure  75 
shows  good  agreement  for  the  »  400  WPPM  injection  case.  The  C * 
500  WPPM  injection  data  also  plotted  on  this  figure  indicate  a 
lesser  developed  profile  and  a  possible  thicker  sublayer. 

When  2000  WPPM  was  injected  through  the  .005-inch  wide  slot 
at  a  velocity  ratio  U/Uavg  *  .280,  the  freestream  turbulence 
intensity  was  3.0%.  Figure  76  presents  a  comparison  of  the 
velocity  profile  development  for  the  2000  WPPM  injection  case  with 
the  comparable  water  injection  case.  These  data,  taken  at  x  a  3.75 
inches  (2  inches  downstream  from  the  injection  slot) ,  indicate 
a  disturbance  effect  due  to  the  higher  injection  velocity  for 
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Figure  76.  Developing  velocity  profiles  for  injected  polymer 
concentration  ■  2000  WPPM 
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the  narrow  slot.  The  tendency  of  the  thick  2000  WPPM  injection 
concentration  to  stay  together,  coupled  with  the  high  injection 
velocity,  contribute  to  the  disturbance  of  the  boundary  layer. 

Here  again,  at  x  *  18  inches,  the  polymer  velocity  profile 
compares  well  with  the  water  turbulent  velocity  profile  in  the 
outer  region  of  the  boundary  layer.  Closer  to  the  wall,  the 
polymer  velocity  profile  departs  from  the  water  profile  and  i3 
less  full.  This  anomaly  in  the  near  wall  region  may  be  due  to 
a  thickening  of  the  laminar  sublayer. 

Figure  77  presents  a  comparison  of  the  turbulence  intensity 
distributions  at  x  ■  18  inches  for  the  2000  WPPM  polymer  injection 
case  and  a  similar  water  case  at  the  same  injection  velocity  ratio 
((V^/Vavg  ■  .280)  and  slot  width  (s  ■  .005  inches).  The  charac¬ 
teristic  lowering  of  the  peak  value  of  the  turbulence  intensity 
and  shifting  of  its  location  away  from  the  wall  are  also  shown. 

The  lowering  and  redistribution  of  the  primary  turbulent 
activity  into  a  broader  band  found  in  this  study  is  supported 
by  the  work  of  Reischman  (1973).  However,  due  to  experimental 
differences  only  qualitative  comparisons  can  be  made  with  the 
work  of  previous  researchers.  Rudd's  data,  typical  of  Logan's 
and  Rumor  and  Sylvester's  measurements,  are  presented  in  figure 
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Figure  77.  Turbulence  intensity  distribution  for  injected  polymer 
concentration  ■  2000  WPPM 
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78  taken  from  Reischman  (1973).  However,  Reischman  used  a 
rectangular  channel  while  the  others  conducted  their  experiments 
in  a  square  duct. 

A  controversy  exists  concerning  the  so-called  "sublayer 
thickening"  which  is  represented  by  the  velocity  profile  remaining 
linear  to  a  y+  value  of  15  to  20  instead  of  the  solvent  value  of 
8  to  9.  There  have  been  five  investigations  (including  this  one) 
that  have  made  velocity  measurements  at  y+  <  10.  Table  13  from 
Reischman  (1973)  indicates  the  extent  of  the  linear  curve  and  the 
lowest  y+  location  at  which  velocity  measurements  were  made  in 
each  prior  investigation.  Data  from  this  investigation  have  also 
been  included.  Reischman  alone  concluded  that  the  sublayer  does 

■f  +  + 

not  thicken  based  on  the  maximum  value  of  y  for  which  u  ■  y  . 
Previous  results  achieving  larger  values  of  y+  which  indicated  a 
thicker  sublayer  were  held  in  serious  question  due  to  experimental 
facilities  and  technique.  This  problem  was  discussed  in  section 
II  which  reviewed  experimental  facilities.  It  has  been  found 
that  in  the  rectangular  channel  used  in  this  research  linear 
values  of  y+  up  to  15  have  been  achieved  with  Polyox  WSR-301. 
Figures  79  through  83  present  velocity  profiles  plotted  in  terms  of 
inner  variables  and  compared  to  the  law  of  the  wall.  It  is 
recommended  that  future  investigations  concentrate  their  efforts 
in  the  region  of  y+  from  7  to  50  to  provide  additional  data  at 
the  interface  of  the  sublayer  and  the  buffer  region  or  elastic 
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Figure  78.  Turbulence  intensity  distribution 


Comparison  of  the  law  of  the  wall  with  experimental  data 
or  Injected  polymer  concentration  C  -  100  WPPM 


or  injected  polymer  concentration  C  **  200  WPPM 


Comparison  of  the  law  of  the  wall  with  experimental  data 
or  Injected  polymer  concentration  C  ■  400  WPPM 
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TABLE  13 


COMPARISON  OF  NEAR  WALL  LASER  ANEMOMETER 
MEASUREMENTS  IN  DRAG  REDUCING  FLOWS 


Author 

Maximum  y+  + 
for  Which  u  *y 

Lowest  y+ 
Measurements 

Lazar 

10 

10 

Rudd 

19 

8 

Kumor  & 

16 

3 

Sylvester 

Reischman 

9 

3 

Present  Work 


15 


7 
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sublayer  which  joins  the  viscous  sublayer  and  the  elevated  log 
region.  Figure  84  indicates  the  large  variation  in  the  velocity 
profile  in  this  region.  Profiles  in  the  outer  region  tend 
toward  the  turbulent  while  becoming  more  laminar  close  to  the 
wall. 

Comparison  of  Analytical  Predictions  With  Experimental  Results 

Typical  boundary  layer  growth  patterns  are  presented  in 
figure  85  through  88  for  various  concentrations  injected  at  a 
constant  rate  of  V  /V avg  ■  .0174.  There  is  excellent  agreement 
between  the  experimental  data  and  the  theoretical  profile  predicted 
by  the  mathematical  flow  model  utilizing  the  effective  Reynolds 
number  analogy.  The  gradual  increase  in  the  boundary  layer  height 
for  the  500  WPPM  injection  case  may  be  explained  in  terms  of  the 
water  temperature  and  associated  solvent  viscosity.  For  this 
case,  the  water  temperature  was  the  coldest  at  57°F.  The  resulting 
higher  kinematic  viscosity  resulted  in  lower  axial  Reynolds  number 
affecting  transition  and  the  growth  of  the  boundary  layer. 

Skin  friction  calculations  as  a  function  of  local  polymer 
wall  concentration  were  based  on  three  techniques:  (1)  slope 
of  the  velocity  profile,  (2)  Clauser  technique,  and  (3)  flow 
model.  The  results,  indicating  good  agreement,  are  presented  in 
tables  14  through  19.  Skin  friction  computation  from  the  slope 
of  the  velocity  profile  at  the  wall  is  relatively  straightforward 
and  will  not  be  discussed  in  detail.  This  method  is  sensitive, 
however,  to  the  value  of  the  velocity  near  the  wall.  Velocity 
measurements  must  be  made  in  the  vicinity  of  the  laminar  sublayer 
and  are  susceptible  to  velocity  gradient  biasing  in  the  laser 
measuring  volume. 
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NORMALIZED  VELOCITY  U/UFS 

Figure  84.  Comparison  of  normalized  velocity  profiles  with 
l/5th  power  law  at  station  x  *  18  inches 
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Comparison  of  experimental  boundary  layer  growth  with  analytically  predicted 
values  for  injected  polymer  concentration  C  =  100  WPPM 


>arison  of  experimental  boundary  layer  growth  with  analytically  predicted 
values  for  injected  polymer  concentration  C  =  200  WPPM 
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TABLE  14 


COMPARISON  OF 

SKIN  FRICTION 

COEFFICIENTS 

FOR  WATER  CASE 

x (Inches) 

16 

17 

18 

R 

354252 

376793 

402064 

ex 

Cf  (Model) 

.00458 

.00453 

.00447 

Cf  (Clauser) 

.00450 

. 00440 

.00425 

Cf  (Velocity 
t  Profile) 

.00497 

.00478 

.00461 

C 

0 

0 

0 

w 

TABLE  15 

COMPARISON  OF  SKIN  FRICTION  COEFFICIENTS  FOR 
C1  -  100  WPPM  OF  POLYMER 


x  (inches) 

16 

17 

18 

R 

ex 

330,812 

352,991 

375,215 

Cf  (Model) 

.00347 

.00380 

.00426 

Cf  (Clauser) 

.00325 

.00335 

.00350 

Cf  (Velocity 
Profile) 

.00405 

.00406 

.00424 

10 


4.6 


3 
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TABLE  16 

COMPARISON  OF  SKIN  FRICTION  COEFFICIENTS  FOR 
Cf  -  200  WPPM  OF  POLYMER 


x  (inches) 

16 

17 

18 

R 

345601 

346215 

365801 

ex 

Cf  (Model) 

.00327 

.00345 

.00359 

C^  (Clauser) 

.00300 

.00340 

.0040 

Cf  (Velocity 
r  Profile) 

.00356 

.00395 

.00420 

C 

15 

5 

1 

V 

TABLE  17 

COMPARISON 

OF  SKIN 

FRICTION  COEFFICIENTS 

FOR 

Cf  -  400 

WPPM  OF  POLYMER 

x  (inches) 

16 

17 

18 

R 

322,476 

343262 

365  ,325 

ex 

Cf  (Model) 

.00337 

.00354 

.00350 

C^  (Clauser) 

.00220* 

.00265* 

.00280* 

C.  (Velocity 
x  Profile) 

.00345 

.00354 

.00374 

C 

80 

60 

35 

w 

*  Clauser  skin  friction  values  only  valid  for  values  of  Cf  ^  200 
WPPM  polymer. 
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Clauser  Skin  Friction  Techniques 

The  Clauser  technique,  presented  In  1954,  and  discussed  by 
C.C.  Lin  (1959),  bears  some  discussion.  Utilizing  the  law  of  the 
wall  and  expressions  (77) ,  (78) ,  and  (79) ,  Clauser  obtained  a 
family  of  curves  of  U/UFS  versus  log  (Re^)  with  the  skin 

friction,  C^,  as  a  parameter,  so  that 

U 

Ur  -  —  (77) 


and 


This  family  of  curves  for  water  flow  utilizing  the  value 
of  B  •  7.2  is  shown  in  figure  89.  Application  of  this  figure 
to  a  determination  of  C^  merely  requires  the  placing  of  a 
measured  velocity  distribution  thereon  and  reading  off  the 
value  of  C£,  interpolating  where  necessary.  Adapting  this 
technique  to  polymer  flows  requires  that  the  appropriate  form  of 
the  law  of  the  wall  be  used,  resulting  in 
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Figure  89.  Plot  of  Clauser  family  of  curves  for  water  flow 
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where  3  *  a  CY.  (81) 

w 

Figure  90  demonstrates  the  application  of  this  technique  to 
polymer  flows.  A  new  family  of  curves  must  be  obtained  for  each 
polymer  concentration  considered. 

Designed  for  fully-developed  turbulent  flow,  application  here 
has  been  attempted  for  developing  flow.  Additionally  the  AB  shift 
for  polymer  flow  is  assumed  to  be  of  the  form  given  by  equation 
(81).  For  this  reason,  the  skin  friction  values  are  presented 
for  comparison  purposes  only.  Additional  experimentation  with 
fully-turbulent  flows  over  a  range  of  polymer  concentrations  should 
be  performed  to  further  develop  this  technique  since  the  present 
data  indicate  much  promise. 

Figures  91  through  96  present  variations  in  the  value  of  the 
local  skin  friction  along  the  plate  for  various  injection  polymer 
concentrations  as  predicted  by  the  flow  model.  Also  plotted  on  these 
figures  are  values  of  skin  friction  determined  from  data  taken  at 
the  plate.  The  comparative  agreement  is  excellent.  The  lower 
values  of  determined  by  the  modified  Clauser  technique  may  be 
attributed  to  the  developing  aspect  of  the  velocity  profiles  and 
their  near-wall  anomaly  due  to  the  polymer.  The  flow  model 
prediction  compares  very  well  with  the  velocity  profile  data  in 
all  cases.  During  transition,  a  linear  intermittency  function 
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gure  90.  Plot  of  Clauser  family  of  curves  for  injected  polymer 
concentration  *  100  WPPM 
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Figure  91.  Variation  of  local  skin  friction  along  the  plat 
compared  with  experimental  results  for  water  flow 
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Figure  92.  Variation  of  local  skin  friction  along  the  plate 
compared  with  experimental  results  for  injected 
polymer  concentration  -  100  WP PM 
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Figure  94.  Variation  of  local  skin  friction  along  the  plate 
compared  with  experimental  results  for  injected 
polymer  concentration  C  =  400  WPPM 
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Figure  95.  Variation  local  akin  friction  along  the  plate 
compared  with  experimental  results  for  injected 
polymer  concentration  C  ■  500  WPPM 
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Figure  96.  Variation  of  local  skin  friction  along  the  plate 
compared  with  experimental  results  for  injected 
polymer  concentration  C  -  800  WPPM 


is  assumed  with  completion  of  transition  to  turbulent  flow  at 
Re^  *  350,000.  The  physical  limitation  of  the  test  setup  and 
the  limited  length  of  the  plate  precluded  the  taking  of  additional 
data  at  larger  length  Reynolds  numbers.  More  research  in  this 
region  with  fully- developed  turbulent  flow  is  recommended  to 
complete  the  development  of  the  flow  model. 

The  effective  Reynolds  number  analogy  values  of  the  local 

skin  friction  agree  very  well  with  the  values  predicted  by  the 

analytical  model.  Figure  97,  plotted  for  the  conditions  of  the 

200  WPPM  injection  case,  typically  represents  the  other  injection 

conditions.  Curves  of  versus  Rex  are  plotted  with  C  as  a 

parameter.  To_ obtain  a  value  for  the  skin  friction,  the  plot 

should  be  entered  with  the  value  of  Re  while  the  value  of  C. 

x  f 

is  taken  off  at  the  appropriate  value  of  the  polymer  wall  concentration. 
Utilization  of  this  technique  to  aid  in  the  solution  of  equation  (63) 
has  demonstrated  its  suitability  for  application  to  polymer  flows. 

Polymer  Diffusion  Patterns 

Typical  polymer  wall  concentration  distributions  normalized 
with  injection  concentration  and  plotted  against  normalized  axial 
length  are  shown  in  figures  98  through  99.  The  polymer  flows 
injected  into  the  boundary  layer  showed  suppressed  diffusion 
characteristics  in  the  transition  region  0.3  <  X/L  <  0.7. 


Effective  Reynolds  number  analogy  values  of  local 
riction  verses  Reynolds  number  for  injected 
solymer  concentration  C,  “  200  WPPM 


NORMALIZED  WALL  CONCENTRATION  C/Cw  n  NORMALIZED  WALL  CONCENTRATION  C/Cw 
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98.  Normalized  polymer  wall  concentration  distribution 


Figure  99.  Normalized  polymer  wall  concentration  distribution 
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Higher  concentrations  of  injected  polymer  resulted  in  higher 
plate  trailing  edge  concentrations.  The  500  WPPM  injection  case 
showed  the  nighest  wall  concentrations  for  the  largest  distance  along 
the  plate.  It  is  believed  that  this  distribution  was  influenced  by  a 
Reynolds  number  effect  due  to  the  lower  water  temperature  of  this 
run.  The  lower  temperature  produced  a  higher  value  of  the  kinematic 
viscosity.  The  resulting  lower  Reynolds  number  had  an  additional 
delaying  effect  on  transition. 

The  polymer  appears  to  affect  transition  by  delaying  its 
onset,  dampening  the  turbulent  transport  of  material  away  from 
the  wall  region,  and  stretching  out  the  transition  region. 

Figures  100  and  101  show  the  effect  of  polymers  on  transition. 

The  maximum  velocity  difference  between  laminar  and  turbulent 
velocity  profiles  occured  at  a  height  of  y  =  .020  inch  above 
the  plate.  Axial  velocity  distribution  profiles  at  this  height 
along  the  plate  show  the  characteristic  increase  in  velocity  at 
transition  to  turbulent  flow  for  the  water  injection  case.  Injec¬ 
tion  of  polymer  showed  no  such  characteristic  increase  in  velocity 
indicating  delayed  transition. 

Diffusion  patterns  normal  to  the  plate  were  taken  at  five 
stations:  1  =  3.75  inches,  2  =  8.5  inches,  3  =  12.5  inches,  4  =  16.5 
inches  and  5  =»  20.5  inches.  No  concentration  measurements  were 
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Figure  100.  Axial  velocity  profile  at  y  »  .020  inches  for  injected 
polymer  concentration  *  500  WPPM  compared  with  water 
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Figure  101.  Axial  velocity  profile  at  y  *  .020  inches  for  injected 
polymer  concentration  ■  2000  WPPM  compared  with  water 
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generally  obtainable  above  y  ■  .010  inch  to  .015  inch  above  the 
plate  for  stations  1  and  2,  indicating  very  little  to  no  diffusion 
in  these  laminar  regions.  Figures  102  through  107,  therefore, 
compare  the  three  concentration  measurement  stations  along  the 
plate  at  x  ■  12.5,  16.5  and  20.5  inches.  The  injected  polymer 
concentrations  of  «  100  and  200  WPPM  Indicate  a  continuous 
gradual  decay  in  concentration  at  stations  3  and  4,  with  rela¬ 
tively  high  concentration  levels  in  the  near  wall  region.  The 
water  injected  case  at  these  stations  shows  rapid  concentration 
decay  in  the  near  wall  region.  At  station  4,  the  diffusion  pattern 
above  y  «*  .015  inch  takes  on  a  linear  decay  characteristic 
found  in  the  turbulent  flow  outer  region  at  station  5.  The 
concentration  distributions  normal  to  the  wall  at  station  5  for 
polymer  injection  of  *  100  and  200  WPPM  are  similar  to  the  water 
case  at  the  same  station.  This  is  due  to  the  very  low  wall 
concentration  at  this  station. 

The  reduced  diffusion  characteristics  of  polymer  flows  may 

be  seen  by  plotting  the  growth  of  the  diffusion  boundary  layer 

indicated  by  a  characteristic  height,  X  .  This  characteristic 

c 

height ,  X  ,  is  defined  as  that  height  at  which  the  concentration 
c 

is  equal  to  one-half  of  the  wall  concentration.  Increasing 
concentration  of  injected  polymer  reveals  a  reduced  diffusion 
characteristic  away  from  the  wall  and  a  thicker  layer  of  higher 
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103.  Normalized  polymer  diffusion  pattern  for  injected 
polymer  concentration  *  100  WPPM 
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normalized  concentration  c/cw 

Figure  104.  Normalized  polymer  diffusion  pattern  for  injected 
polymer  concentration  *  200  WPPM 
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Figure  105.  Normalized  polymer  diffusion  pattern  for  injected 
polymer  concentration  *  400  WPPM 
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concentration  in  the  near  wall  region  as  shown  In  figure  108. 

The  quantity  of  polymer  or  mass  flow  rate  Injected  Into  the 
boundary  layer  may  be  considered  to  be  given  by  the  flow  rate 
of  Injection  times  the  concentration.  Figure  109  Is  a  plot  of 
diffusion  boundary  layer  height  versus  axial  length  Injected  at 
various  flow  rates  and  concentrations  In  combination  to  produce 
a  constant  quantity  of  polymer  Injected  Into  the  boundary  layer 
of  QiCi  '•=  3,000.  The  lowest  rate  of  injection  at  the  highest 
concentration  produces  the  slowest  growth  of  the  diffusion  boundary 
layer.  The  effect  of  Increasing  the  quantity  of  polymer  Injected 
Into  the  boundary  layer  is  shown  in  figure  110  where  the  value 
of  ■  16,000.  Again,  the  lowest  rate  of  injection  at  the 

highest  concentration  produces  the  slowest  growth  of  the  diffusion 
boundary  layer.  The  increased  amount  of  polymer  injected  Into 
the  boundary  layer  greatly  suppressed  the  diffusion  of  polymer 
away  from  the  wall.  The  effect  of  increasing  injection  concen¬ 
tration  flux  into  the  boundary  layer  may  be  seen  in  figure  111, 
which  may  be  considered  to  present  a  measure  of  the  injection 
"displacement"  thickness.  At  constant  velocity  of  injection, 
increased  injection  concentration  flux  greatly  suppressed  the  growth 
of  the  diffusion  boundary  layer.  The  lowest  injection  velocity 
produced  the  lowest  diffusion  boundary  layer.  In  the  final  zone, 
diffusion  was  found  to  vary  with  an  exponent  of  .6  as  shown  in 
figure  112.  Sirmalis  had  predicted  that  the  exponent  for  similarity 
profiles  would  be  lower  than  his  value  of  .75  as  he  was  unable  to 
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obtain  data  within  .01  Inch  of  the  wall.  Similarity  profiles  for  the 
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Figure  109.  Diffusion  boundary  layer  growth  for  C^C 
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Figure  112.  Diffusion  boundary  layer  growth  in  the  final  zone 
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Figure  113.  Diffusion  boundary  layer  growth  in  the 
transition  region 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


General 

The  object  of  this  research  was  threefold:  (1)  to  develop 
an  experimental  facility  and  techniques  for  the  study  of  the 
effects  of  drag-reducing  polymers  on  submerged  flat  plate  flow, 

(2)  to  qualitatively  define  the  effects  of  polymers  on  boundary 
layers  in  external  flows,  and  (3)  to  verify  predictive  techniques  by 
experiment.  A  flat  plate  water  tunnel  facility  was  constructed 
and  laser  Doppler  anemometry  techniques  were  developed  for  boundary 
layer  velocity  measurements.  The  flow  facility  is  capable  of 
variable  injection  velocity  and  concentration  of  polymer  with 
boundary  layer  sampling  at  several  heights  throughout  the  boundary 
layer  and  at  several  axial  stations.  The  injection  process  was 
unique  because  injection  occurred  very  close  to  the  leading 
edge  into  laminar  flow,  then  developing  through  transitional 
to  turbulent  flow.  Application  of  laser  Doppler  techniques 
provided  non-disruptive  velocity  measurements  in  a  low  turbulent 
level  facility  for  a  wide  range  of  polymer  concentration  and 
velocity  of  injection. 
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Analytical  Predictive  Methods 


The  case  of  developing  flow  with  polymer  injection  has 
received  little  attention  for  external  flows.  Predictive  methods 
are,  therefore,  not  well  developed.  Past  works  on  the  diffusion  of 
polymer  boundary  layers  have  centered  on  the  far  downstream 
region  in  turbulent  flow  where  the  polymer  attained  the  same 
diffusion  characteristics  as  the  medium. 

An  analytical  model  was  developed  for  comparison  with  experi¬ 
mental  data.  The  model  is  based  on  integral  boundary  layer  analysis 
combined  with  velocity  profile  relation  for  flat  plate  flow  at 
zero  pressure  gradient  and  variable  polymer  wall  concentration. 
Accurate  prediction  of  local  skin  friction  and  boundary  layer 
height  as  a  function  of  variable  local  wall  polymer  concentration 
has  been  demonstrated  over  laminar,  transitional,  and  turbulent 
flow  regions  by  comparison  with  experimental  data. 

Polymer  Boundary  Layer  Characteristics 

Polymer  additives  appeared  to  have  no  effect  in  the  laminar 
flow  region.  However,  increasing  polymer  concentrations  delayed 
and  stretched  out  the  transition  to  turbulent  flow.  The  effect 


of  polymer  additives  is  predominant  in  the  region  of  maximum 
turbulence  production. 


233 


Streamwise  turbulence  Intensity  distributions  are  lower  In  peak 
value  and  displaced  further  from  the  wall  over  a  broadened  buffer 
region  near  the  wall  than  found  in  comparable  solvent  flows. 
Near-wall  velocity  measurements  have  shown  that  the  laminar 
sublayer  thickens  with  the  introduction  of  polymer.  A  thicker 
velocity  profile  sublayer  reduces  the  velocity  gradient  at  and 
near  the  wall.  The  wall  shear  stress,  proportional  to  the  velocity 
gradient  at  the  wall,  also  decreases,  thus  reducing  the  frictional 
drag  of  the  surface.  Although  thickening  the  sublayer,  polymers 
cause  the  overall  boundary  layer  to  become  thinner.  Skin  friction 
prediction  by  Clauser's  technique  has  been  shown  applicable,  with 
modification,  to  polymer  flows. 

f 

Diffusion  Characteristics 


Various  polymer  concentrations  were  injected  into  a  laminar 
flat  plate  boundary  layer  at  various  injection  rates.  Transition 
to  turbulent  flow  was  delayed  and  extended.  Kowalski's  (1974) 
explanation  that  polymer  molecules  seem  to  act  as  barriers 
restricting  communication  across  the  boundary  layer,  forcing  the 
liquid  to  flow  in  a  semi- laminar  manner  is  most  appropriate  here. 
The  growth  of  the  diffusion  boundary  layer,  indicated  by  a 
characteristic  height,  decreased  with  increasing  concentration 
and  quantity  of  polymer  injected  into  the  boundary  layer. 
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The  increased  concentration  of  polymer  showed  reduced  diffusion 
in  Che  transition  region  resulting  in  an  extended  initial  mixing 
zone.  This  effect  is  also  supported  by  the  work  of  Sirmalis  (1976). 

Conclusion  Summary 

An  experimental  facility  has  been  developed  that  allows 
laboratory  measurements  of  polymer  drag  reduction  concentration 
profile  distribution  for  various  injector  geometries,  velocities 
of  injection,  and  injection  concentrations.  Hydrodynamic  properties 
of  the  boundary  layers  were  measured  using  laser  Doppler  anemometry 
techniques.  Experiments  indicated  thickening  sublayers  and  overall 
thinning  of  the  boundary  layer  with  polymer  concentration.  Turbu¬ 
lence  intensity  distributions  were  lower  in  peak  value  and  distri¬ 
buted  over  a  broad  region.  Diffusion  of  polymer  injected  into 
the  laminar  flow  at  the  leading  edge  of  the  flat  plate  delayed 
transition  and  spread  the  transition  region  out.  Very  little 
diffusion  occurred  in  the  laminar  region.  Diffusion  in  the 
turbulent  region  indicated  by  a  characteristic  height  of  the  diffusion 
boundary  layer  decreased  with  increasing  concentration  and  quantity 
of  polymer.  Predictive  methods  were  developed  building  on  reliable 
theories  and  techniques  for  the  prediction  of  local  skin  friction 
as  a  function  of  polymer  wall  concentration.  Comparison  with 
experiment  has  been  performed  and  good  agreement  achieved. 
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Recommendations 


1.  Additional  testing  should  be  performed  at  boundary  layer 
heights  between  values  of  y+  equal  to  7  to  50.  This  region  would 
supply  additional  data  to  define  the  increasing  height  of  the 
sublayer  with  polymer  concentration,  as  well  as  the  polymer- induced 
AB  shift.  Low  concentrations  of  injected  polymer  below  50  WPPM 
should  also  be  tested. 

2.  Application  of  the  modified  Clauser  skin  friction  calcula¬ 
tion  technique  to  polymer  flows  should  be  continued  and  will 
require  additional  data  in  the  above-mentioned  region. 

3.  Application  of  traversing  laser  Doppler  techniques  should 
be  extended  to  include  rotation  of  the  measuring  volume  by 

90°  to  obtain  distribution  of  the  v'  turbulence  intensity  dis¬ 
tribution  allowing  computation  of  Reynolds  stress. 

4.  Reduction  of  flow  noise  indicates  that  frequency  spectrum 
analysis  should  be  performed  to  determine  which  frequencies  are 
attenuated  by  various  polymer  concentrations. 

5.  The  experimental  facility  should  be  lengthened  to  allow 
a  more  gradual  and  natural  transition  to  turbulent  flow. 
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6.  Multiple  horizontal  injection  stations  along  the  plate 
should  be  tested  with  lower  concentration  of  injection  to  verify 
the  disentanglement  theory.  Optimization  of  the  injection 
system  should  sd.so  include  pulsed  low  concentration  injection 

to  determine  persistence  and  possible  time  scale  correlation. 

7.  Modification  of  the  tunnel  passage  will  allow  the  study 

of  polymer  effects  in  boundary  layers  with  adverse  pressure  gradients. 
The  capability  for  acoustic  modulation  of  the  laser  beams  should 
be  added  to  the  test  setup  to  measure  flow  reversals  using  non- 
disruptive  laser  Doppler  anemometer  techniques. 
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Appendix  A 

DEVELOPMENT  OF  SKIN  FRICTION  RELATIONS 


1 


238 


The  law  of  Che  wall  incorporating  a  polymer  term,  according 
to  Meyer  (1966) ,  is  given  by 

u+ "  i^)iny+ + b + a  (m) 

where  a  •  2.3,  y  ■  .5. 

The  conservation  of  mass  and  Che  momentum  equations  for  flow 
over  a  flat  plate  are  given  by 

3_  (pu)  +  3_  (pv)  -  0  (A2) 

3x  3y 

and 

For  a  flat  plate,  it  may  be  assumed  that  the  pressure  gradient 
is  negligible,  such  that 


3P  .  n  (A4) 

3x 


Therefore,  equation  (A3)  becomes 
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„  /au\  ,  _  /3u\  3t 
\»*)  \3y)  "  3y* 


(A5) 


Ic  may  further  be  assumed  that  y  ■  0,  p  •  £,  and  v  «■  t  m  ^  are 


all  constants. 


<A6) 


Equation  (A2)  becomes 


3u  3v 
3x  3y 


0. 


CA7) 


An  expression  for  v  becomes 


3v  m  _  3u 
3y  3x 


dv 


7. 

v(y)  -  »  -  -  |® 


(A8) 


Changing  to  law  of  the  wall  variables  we  have 


+ 

u  ” 


(A9) 
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u  ■  V  u+ 

(A10) 

3u  .  3(V*u) 

3x  3x 

(All) 

* 

+  yV 

y  ■  * — 

7  V 

(A12) 

* 

dy+  -  •“  dy 

(A13) 

dy  *  ^7-  dy+ 

V 

(A14) 

,* .  v*w .  vx 

(A15) 

Combining  equations  (A8j  All,  A14) , 

equation  (A8)  becomes 

c-)”* 

(A16) 

Changing  equation  (A5)  into  law  of 

the  wall  variables  using 

equations  (A9)  thru  (A16)  we  have 


*  +  a  *  + 
P(V  U  )g^(V  U  )  + 


*  + 
IS1  a  ) 
3y 


i 

I 

< 


3x 

3y  • 


(A17) 
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Substituting  equation  (A20)  into  equation  (A18) , 


+.  f  +  3V*  .  _»  3y+/3u+\  .  „»  9CW  /3u+\l 

,(vu)iu  ■sr*v  IT^JJ 

.  + 

I  /  3  ,  *  +.  v  .  +  3  +s  3r 

[i  3j(,“>7d!rj3?(7“)  'I?' 


~P 


(A21) 


Obtaining  an  expression  for  ~  we  have 


JJu  ,  3(V*u+) 
3y  3y 


*  + 

+  3V  .  *  3u 

“  1T  +  7  IT  ’ 


(A22) 


However. 


*  *,  . 
V  «  V  (x) 


(A23) 


and  hence 


3V 

3y 


(A24) 


Equation  (A22)  then  becomes 


»(vV).  ,  v*  ini 

3y  3y 


(A25) 


Changing  3y  to  3y  we  have 
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but 


3<V*u) 

ay 


CA26) 


0 


and  therefore 

3(vV) 

3y 


Equation  (A21)  becomes 


P 


*  + 
(V  u+) 


av 

3x 


-P 


(A27) 


As  v 


t ,  equation  (A6)  and  V 


V*(x) 


equation  (A23)  may 


then  be  removed  from  within  the  integral  yielding 


P(V 


'  +.  +  3v 

u  >  u  a T 

m 


-P 


ii 

3y* 


(A28) 


After  substitution  of  equation  (A20)  with  equation  (A28)  we  have 
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Substituting  into  equation  (A29)  we  have 


*  u 

Dividing  then  by  pV  and  remembering  that  v  ■  — 
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(A30) 


+ 

Multiplying  by  dy  ,  expanding  terms  and  intergrating  we  have 


From  equation  (A6)  and  (A22)  up0n  rearringing  equation  (A31) 

and  integrating  through  the  boundary  layer  such  that  t(y+)  ■  0,  we  have 


+ 


+ 
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r  y+ 


■/j  i  <•*(«-■>  •  •'(«-*) 


v  V 


Defining  groups  of  terms  we  have 


+ 


Therefore  equation  (A32)  becomes 


3V 

3x 


•(G)  + 


* 

V  H 


* 

V  J  - 


N 


1 

~  p  V 


Recalling  the  law  of  the  wall  equation,  (A33)  becomes 
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(A32) 


(A33) 


(A34) 


(A35) 


(A36) 


(A37) 


(A38) 
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G  ■  f0(k +  53  * 2-3  7 *  (7*))  d*+ 


t 

J  |  J2(lny+)2  +  (5.5) 2  +  (2.3)2C^*141h  +  2(5.5)  £  lny+ 


+  2(2. 3)|  lny+  C^57  In  +  2(5.5)  (2.3)C^57  In  /_  I  dy+.  (A39) 

V0  VQ*  ) 


”2  f  Clny+)2dy+  +  (5.5)V  +  (2. 3) 2  cj*14ln2  (~\  y+ 
K  0  \Vrt/ 


+  2^1  /  lny+  (d/)  +  ^  ln  fLA  Jr 

*■  \V0/  0 

A 

+  2(5.5)(2.3)C^57ln  (^\  y+ 


(A40) 


Note:  (ln  x)2  dx  »  x(ln  x)2  -  2[x(ln  x-1)]  *  x[ln2x-2  ln  x  +2] 

*  x[ln2x-2(ln  x-1)]  (A41) 


hence 


G  -  -2  {y+[ln2y+  -  2 (In  y+  -  1)]|  +  (5.5)2  y+ 
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G 


+ 


,,  „,2  1.14  T  2 

(2.3)  C  In 

w 


+ 

7 


+  [y+(in  y+  -  D3  +  ^y1  ^  ln 

Vvo  ' 


* 

+  2(5.5)  (2.3)  C^57ln  y+ 


y+|i2  [ln2y+  -  2 (ln  y+  -  1)]  +  |(ln  y+  -  1)  js.5  +  2.3 
R 


-  1) 

(A42) 


+ 


12(5.5)  +  2.3 


(A43) 


Evaluating  term  H,  we  have 


H 


(A44) 


3y+  3_  ,  V*  £  /^L\ 

3x  "  3x  v  v  \3x  j 


(A45) 


7  ln  y+  +  5.5  +  2.3 

'  IN. 


(A46) 
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H 


■  i£ln7+  +  5-5  +  2-3!TI57Wr 


v  lr]g+j  dy+-  (A47> 


Now  since 


+  V 

y  ■  y  — 

M 


+ 

L  ,  Z_. 

V  * 

V 


Substitute  into  equation  (A47) 

+ 

i 

H 


(A48) 


Ky 


dy+  (A49) 


y+f*l  + 

/  K  ln  y  +  5' 
-/n  L 


5  + 

0  7 


1 _  3V_ 

V*K  3x 


dy 


1  3V  r  1 1  ,  +  .57  /v  \ 

/[k  11  y  +  5-5  +  2-2S  i"  (-*) 


+ 

*  y  r 


.57  .  /V  \|  + 

dy 


*  .  “I 


ffi-y  i  a*  7+ -  „  +  ,.3  +  2.xi’7*  fa' 


V  K 


0  'J 


(A50) 
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Evaluating  term  J: 


J 


3C, 


3x 


'W  ,  + 
dy 


(A51) 


*1 L£k  lny+  +  5-5  +  2*3Cw57ln(7^] 


‘  . 57)  (2 . 3)  C^* 43 dy+ 


I 

(.57)(2.3)^-43ln^!i  jf  (±ln  7+  +  5.5  +  2.3C^71»  ^  <ty+ 

7+  ^ln  y+  -  1)  +  5.5 


(A53) 


•< 


Term  2  becomes 
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1  dV*  + 
K  dx  7  * 


(A 59) 


Now  evaluate  term  N  : 
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where 


(A65) 


,  13V 

K  3x 


k(1°T+ 


1)  +  5.5  +  2. 


(A66) 


4 


(A67) 


5 


1 

K 


<2.3)(.57)<£* 


(A68) 


And  we  then  have 


•T 

kjrf0  i  <■***  -  »  * 5-5  +  2-3Ci5h°(fy  ^ 


1  iL 

K  3x 


J  (y+ln  y+  -  2y+  +  5.5  y+  +  2.3C^5?ln y+j 
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* 

3V_ 
3x  . 


(G) 


In2  y+  -  2(lny+  -  1) 


+  i  <lny+ 


i) 


5.5  +  2.3C^57 


In 


..57. 


..57 


+  2 . 3C^  lnj  — -y j 1 2 ( 5 . 5 )  +  2.3C~'  It 


V*(H) 


+  +  ^  (lny+  -  1) 


+  5.5 


+  2.3 


* 


1  3V 
K  3x 


1)  +  5.5  +  2.3C^57 


V*(J) 


(.57)  (2.3)C^'43ln ^  y+  K(lny+  "  L)  H 


(5.5)2|(A73) 


(A7A) 


5.5 
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1  V*  3CW  -.43 

+  !  •5I2<2-3)  <‘57>  ^  ^ 


(9't 


(A76) 


Combining  into  equation  (A37)  we  have 
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1 


1 


+  7+  £-jWV  +  ^5-5^2’  C2.3)(.57)“-43ln  L_ 

(K  TT 


+  (2.3) (.57) 


Now 


1 

U 


(A77) 


(A78) 


since 


Tw  "  Cf(Js  Pu2> 


(A79) 


and  letting  ^ 

V 


(A80) 


Recalling  equation  (A15) 


* 

V 


P 
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-4 


Hpu2c, 


resulting  in 


ft 


v  -  u  '-f-  . 


from  equation  (A80) 


x  -  — l_  -vf 

fTn  C£ 


c*  ’  ?• 


Therefore  equation  (A79)  becomes 


(A81) 


(A82) 


(A83) 


1 
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Substituting  equations  (A78)  and  (A85)  into  equation  (A77)  and 
combining  terms  we  have 


* .  r 


|  (In  y+  -  1)  +  2(5.5) 


+  2.3  C 


"’"61 


(+  |(ln  y+  -  1)(5.5  -  £) 


i2(ln  y+)2  +  (5.5)2| 


3C 

-y+x-£  (-57)(2 


•3)c»‘43"0 


1  +  -  (In  y+  -  1) 
k 


+  5.5 


1 

k 


(A86) 


The  value  of  y+  In  equation  (A86)  is  the  value  of  y  at  the  edge 
of  the  boundary  layer. 

Therefore  at  the  edge  of  the  boundary  layer  we  have 


u  -  U 


(A87) 
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and 


<$ 


+ 


hence 


+ 

u 


(A88) 


Recalling  equation  (A80) 


X 


U 

* 

V 


X 


+ 


at  y 


(A89) 


Whence  equation  (Al)  becomes  at 


+ 

y 


5 


+ 


X  ■  •p  ln6  +  +  B  +  a  In/— 

k  w  \v0/ 


(A90) 


where  a  *  2.3  and  y  •  .5 


and 


In  5 


KX  -  KB  - 


K3 


CA91) 
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where 

8-acJ 


and 


K(X-B) 

e 


Equation  (A86)  then  becomes 


5+|31n^-*[^|<ln  5+- 


+  2B  +  Z  In' 


*  1 


0  -I 


+  |(ln  <5+  -  1)(B  -  i)  +  ~2(ln  5+)2  +  B2  j 


^  .te  „  2., 

r»  »0 


1  +  K  (In  «+  -1  )  +  | 
IT  K 


1  \  3C, 


_W 

3x 


Multiplying  by  dx  and  integrating,  we  have 
X(X)  , 


R  (X)  -  R  (X  )  -  f  5+  {  Bln  — . 
e  e  o  J  * 

X(X  )  1  V0 


f(ln<5+  -  1)  +  2B  + 


+  |(ln6+  -  1)(B  -  +  — 2 (In 6  +  )2  + 

1C 


CA92) 


(A93) 


.  (A94) 
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H 


For  injection  of  polymer  at  or  close  to  the  leading  edge  we  may 

assume  X  *  0.  hence  Re(X  )  ■  0  and  C  (X  )  *  C. .  Equation  (A95) 
o’  o  w  o  i 

is  numerirxLly  integrated  in  the  flow  model  described  in  chapter  V. 


APPENDIX  B 

COMPUTER  PROGRAM  DATAPAC  AND  SAMPLE  LISTING 
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3V6?LATOi» 


COMMON/ lOVt*/ 15  NO.  I  PASS 


COMMON  CXtOI.UI.SU, SLOT 
laaMQM  mi T*.nti.Tii ..iri.uy.MTfiu^i 


WCmlliM,  *■■  *!T- 


22(31.21  .TPk.USU31.23 
common -cm  m-2) _ 


CIMCNSZON  *CCM(31),»Y(31> 
-JICIU-UUMOA _ 


PCMINO  1 

■  ip  iiPiss.ga.n  sa  io-io- 
!PASS*1 

ititL  ll^Caui&.ICMiUJ- 


*210  lC.NPk.OT.NCM,  (MClSe(X)  .Xal.lll 
-ii  r  as  mat  1 32131  _ 


1210  2«.UP**TU*>,3K*t,19> 
2B.-S0*M»T(1C14) 


*210  IPMT.lUNITIOd, 2J, I«l, KOIAS) 

-TEAa-Ig.’IT.^Ukf.TA  IX. 2 ) .  I»1Uk3I^UJ _ 


*210  30 . ( YPHJS1 ( C.ZI , X»1 *<UA1LI 
—  *EAQ-  OO-UUPUiS  ( I.-2 1  .  l»U<lUk 
10  P0«M»T(6Fia.5l 


*210  50. (XAX(JK.K) ,JK*l,9) .NXCHA*(K> 

anxuK  *ju  .u*px»il.au!  caasuiu- 


-o 

- 20 

—  20 

.  —  20. 


*210  50 • (TITL2 ( JK  .<> , JK*1 . 1 1 . NTCHA* (K) 

.FOaniT  19A  3,141  _ 

00  aO  Ol.NCH 

*210-50  .  ( CHASAC-LJ* .  U  ,JK«U9I  wHCHAJ*  (XL 
*210  70.N.M.S 

pnanAT  upib.ai _ 


81 

50 


100 


SUT«S 

i£ia-90.ia£Mt.S2JU_ 


r  02  MAT  Ci.-A.2ri0. 5) 

it  i;o£.HT-ax79o  _ao._ioo _ 

C0NTXNU2 
NIMH. 


:rioi.LS.i.£-io»  so  to  i?o 

_J»1— .  _ 


113 

-220 


130 


*210  120.CONC(JI .YCIJ.MM1 .0UM.XC.XC.IST0P 

rOANATUf  U,3,IUl..i5l _ _ 

IPdOTSP.ST.O)  so  TO  130 

_ 


030C02 

-mm 

COO  00* 
-303  845 


000006 


102LTA.2202LT11  .S202LT12  >*2MT0X1 , *2X  .UPS.  Xl.'k  J00  007 

,-nmtnn  TPMT  I  m  I  -rf»fn  Ti-.vi.rutAM  .  i  .»fwijmi...ram1.l  ..-uiai - CO  0-222- 

120)  .T  (311  300  003 

.COMMON-  uxiiau.' «ax  1902)  -x&m-i9.i2)-f.uiaic.n.3iu - ooaoai 


C0MMCN  ACTA (—0 .2) .0UUT2UI31.2} ,*2YL (31.2) .U (31.2) .UPUIS (31.21 .  UNA 
.rrn  t-A-Al  -iijattht  ixi  -ji  -«m  -  »-  .nxi  ri.ni  -?v-.n.-,  r-.<  /t.  ,.  »cr.  n 


000013 

■200 


3XMC.HSI0N  CNU(31 1.01(31)  .02  (31),  03(31 1  ,7*53(31)  ,*CY  ( 311  ,  TX  (  313  .  UH 

.luiuuuixuauuzuuaui)  .xtau-xf  i 

2(31). NMN (12.21 . YNM (31.2) . TPLUS (31,2) 

— ox-Ntmiau-saMccii  jmuuiimu _ 


003310 
-200,111  ■ 


SO  T3  110 
JCSNCaJ  .  . 
0K« C2NC (1) 
.CuCI»C«/CI 


000012 
H2- 
300012 
■100,11.3- 


oaoot- 

-201115- 


300016 
-011017- 
w 00  013 
-0001)19 


300320 
-100,021  - 


300  022 
-300  07.X 


300  02'. 
103121 


000025 

-003122- 


C00C23 
-ICC.  029 


300030 

—n.ir..v 


300032 

-0.00133 


00003- 
caa 135 


300036 

-300137 


300  C  33 
-000139 
3000-0 
-330  0-1 


3000-2 

3000.3 


303  0— 
-310-0.5 


300  0-6 
.3401-7 
000  3-3 
a  n  n  n  -  A 


000350 

-3401=3 


0000:2 
-40015  3 
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00  1*0  <•!• JCONC 

_ 14d._CCt(!C*JUU»CatlC(IU-ZCa - 

00  150  K»l. JCONC 
*»* 

:r(c:ucK,fiMt.s£.o.5)  so  to  ioo 

_  _  151LXQMTINU6 - 

160  NZ*H-Z 

_ S«J3T  *11. _ 

SO  TO  ISO 

_ t70_H«jJt»5 _ _ 

ISO  1*1 

_  190  x£AO  .riiiti.7in.TL/imr) _ 

ZOO  ?OftN*rtSF10.i> 

_ Tin  i  a  i  _sa.  u_  zu _ _ 

:*tn 

_ SO_TO_190 _ 

Z10  CONTINUE 

_ OFWan«<'  LCH < 13 _ 

QV0L*QFLCN*8.33/ <60.  *42. *1 

_ N«£_-t _ 

00  TOO  LL*1 tNPLOT 

_ NOHCAgSILL) 

00  760  NH«1..NC 

_ IPJ  ( N1.  2Q.2  .AMO.  <gLAS.NE.OI  -ANO.  LL.EQ.A3  GQ  TO  750 _ 

IF((NH.£(J.2  .ANO.  <WS1L.N£.0>  .ANO.  11. £0.71  SO  TO  760 

_ :F(11.n£.1I_SO  .to  771 _ 1 - 

vFs»v«t> 

_ . _ UF5  *F9£JU.t 1  *  VF.S _ 

JFSN*UFS»12.*2.3-*.  01 

.  ..  .  3F*5/12. _ 

AS10T«S*N 

»  SI  (37F*ASl.  IT /!<»*. _ 

JX*OI«. 0510207/ (ASLOT*1Z.*S0.1 

_ 3.0_ZZ  3  J»1.N _ 

CNUCJI»J.  609013556-0  J- <7. 17977611E-07 1  * (T( Jl  >  M7. 15490612i-0 

_ 1 _ 3 >  *.tT  ( Jl **Z)  J.  Jfl*77l66£r  111  ».IT  U3T2  11^16^.65377 97au£jlklllTl 

Z  Jl*»<.l 

_ _ ^j(j.N**»*F^«a(Ji»v<ui _ 

J»ATICtJ.N«l*UIJ,H*»/UF5 

220 _ CONTINUE _ 

CO  Z33  K*t.N 

. . .  IF(URATI0(K,«m-,9?».  Zs.a.ZvflvZJO _ 

ZOO  CONTINUE 

Z*»0.  _  JSlTi.ru-l.NHI _ _ _ 

J*«*l 

_ 10_ZJ.C_'*llxN _ 

1*J-" 

.  rz_(.<il«r(i,.Nfii _ 

'JZ(N>*U(L.NH» 

•JR»TI0ZIH*UR»TI0a..,7N» _ 

250  CONTINUE 

_ 30_26  0_.Iii.iJ* _ _ 

T  f I.NN9 •TZ  (II 

_  _ J (I ,NNI*UZ (I)  _  _ 

U«ATI0<I.NN»*UWICZ"(I1 

_ _U*M  II  *.Uf  I  »NH1  •  12.  ?.2 . 5 NJL*. 31 _ 

YN(I3*r(I,.NN»*Z.5fc 

_ 9£Y1I1*UFS»Y( I  .NH1/ (1Z.«CNU (I ) ) _ 

?EYL(I,NNl*ALCGifl  <?EY(I1 ) 

_  260 _ CONTINUE _ 

01(11*1. 

_ _  _OZf  11.*1.£-10 _ 


aaooi* 

•1AOH5'- 
300055 
maB-7 
300  055 
-mans 
300060 
000061 
303062 
3«n<ifei 
00006* 
000065 
000366 
300057 
3  00  Cod 
iiima 
300070 
ono  .in 
000072 
300073 
000071. 
300  C, 7  5, 
000  076 
.0.0007.7 
000078 
300  375 
300060 

.laaasi, 
ooosaz 
-110.063- 
00086* 
-10106  5- 
300336 
300  0*7 
1330;a 
-100014. 
30033d 
COO  036 
300050 
-0  3C.C51. 

0  00  0  52 
_C300i! 

30009* 
.300  395  . 

330096 
-1UC57 
300093 
-100  099 
300100 

_caaiii 

300  132 

COO  113 

COO  13- 
_130  105. 

300  136 
— 100117. . 
000  103 
30010? 
308 113 
_!0.0  111. 

300  112 
—3.11113 


33(11*1. £>10 

._30_Z20.L»tJ< _ 

at  (U*U.-<JNATI0(L.NN1 1 


33<u«j**Tiou.»H«<t.-iuaAria«L,«Ni**i.)  * 

-Ccnlinue _ 

ao  240  L*2.N 

_ .uiujuuwi _ 

continue 

-UUl>U£ixg _ 

30  298  1*1, • 

- *Coato«ccmN2.(«ti - 

*Y(LI*»C(H2,«NI 

..  .ne»nem _ 

CALL  INTEGRAL (N.YY. Qt .QELTA1) 


CALL  INTEGaSL(N,YY,03,0ELTA3l 

_  .,3Q  .iaa  i»i»m _ 

Y0ELTA(I,NN*»r(I.NN1/0ELTA 

_ t  QELI  AULL^HL*  YLL^«U20£LXA1 _ 

Y0ELTA2 1 1 ,NHI*  T  Ct .NH*  /0ELTA2 

_ .  .JHATIflZ  LLkOllLa  iw>Lnir.m,n»n./M _ 

CONTINUE 

_ 

•JA*C*aV0L/A«6A 

_ UAWGJ1*UA\tCil2^i2.3<iiL^aA _ 

?EaiN*<  C2. »MI*<2. *101/12. 

_ MTCtA*L«iAR£A/J.EaXM _ 

HYOIAN«MVOXA*.J0*A 

N12*CEL.TAl/0SL.TAi - 

H32»C£LT13/0eLTA2 

j.£x*uf5**ui/<cnu(ij*ii.x _ 

OENALF*UFS*  CH/2L.*/CNUtl» 

_ --YS3EL  T  A-UFSIOEL  t»/CSKi(  11 _ 

9£C£LTAl*UFS*QEtTAl/CNU (l J 

»e;ELrA2*u(rs*0eLTA2/CNUllL _ 

■2EnvOIA*UAYC*HYOXA/CNU  <  1 J 

.  CFrw«0 ,  .53/  IALOS  <0.  05l*NE*>  ••!.) _ 

CEL  T  AN*0£LT  A  *2.5* 

_ 3EL.TA  li'l*0ELX.Al*2  ,  SL _ 

3ELTA2N«QELTA2»2.3L 

3ELTA3N«0ELTA3*2.5L  .  _ 

*IK»«U) 

<N»rIN*2.3»  _ _ _ 

00  313  1*1. N 

_ TN«Xl*Y«XtNNI/lZ. _ 

«FtI>*«<I) /12. 

_  1£TA(I,NN*  *YN(  I)  rSORT.lUFS/JCNUiiLiFF LD.l >_ 
CONTINUE 

_ 30  323  :*1,N  .  _ _ 

UTAU*UF5*5MT  (CFFW/2. I 

_ auu.T4uii,  (in*  *.(uFs.-.'mj  nnjxiuiau _ 

UPLUS  (I . NN I «U ( I .NH* /U TAU 

_ .  YPLUS.U,HH**'MI,HH**\)T1U/ U2.*CNU.UU _ 

YPLUSLIX.NNI*AL0S13  (YPLUS  (I. (IN)  I 

YMH  (  X  ,NN)  *Y.lX.t  _ 

CONTINUE 

_ IF  (QI.lE.  1  «E“.tfl.)  .3G_LG_1»0 _ — 

CALL  POLY(LANSOA,NCCN,4Y,2) 

_ 30  338  I«1,JC0NC _ 

YLA  N  ( I ,  NNI  *  YC ( X .  NHI  /  LA  NSC  A 
.  PRINT  358 , X  DENT _ 


ooaiK 

JUULL3. 

000120 


300  126 
,,300,121- 
000125 


3  33  1-a 

cnai^r 

300  1.4 
333153 


30315. 
-103-15.5 
303  lie 
-.300157 
303  153 
-100.153 
300  Is 0 
-ItfllaL 
300  152 
300153 
000 15. 
-100,1*3 
303  las 
-lOOliT 
300 164 
-300-11.9 
oooir: 
—13  3 1-~  1. 

303  1*2 
-330  1 ZJ 
300  1.*. 
.130175 
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390  FORMAT  UM1, EON 


rest  number 


.til 


aaoirs 


’(i'll’? 

print  jao.  ufs.ufsh 
»RIJ!T  9  tfl ■ JAVGi UAUfiM 

oocira 

5331T9 

PRINT  J7B.0El.TA,0euriN 

PRTNT  M4..1H  TU.MiriM 

000160 

“HINT  390*aCi.TA2,OcirA2* 

99THT 

0001*2 

1fl3  »  i  * 

PRINT  520.MY0IA  .HYOIAH 

PRINT  W13.H1J 

000134 
mu  »r 

PRINT  -20.H32 
.  _  PRINT  ija.jfy 

000 136 

.  .Af»a««? 

PRINT  £40,R£HAI.F 

PRINT  Sill.  SPHiniA 

000143 

aa.ri  aq 

PRINT  -(.O.REaELTA 
... .  PR  TNT  ..9n.AF.1Fi  Til 

000190 

PRINT  *«3,R£DELTA2 

PRINT  Tlfl  1 

000192 
flan i4i 

P9. INT  47 0*  TUI 
mint  <*9a.  rwm 

COO  194 

’8/11  4C 

PRINT  530.CFFW 

^^UTill  r<5  3A<2Q  QM_  r^ANIt  MH2TE-  5  efl**W-A 

000196 
inn i it 

PRINT  JSO.UTAU 

PRINT  97n.CI  _  .  . 

000196 

M.11  GO 

print  soo.ai 

3R-MT  4in.aFi.au 

000200 

503211 

PRINT  a2a,0R0U 

PRINT  SRfl.Ut 

000  202 

5332*:? 

PRINT  600,3 

PRINT  iia.zu  . .  ... 

0Q9204 

.133219 

PRINT  SMfl.CHCI 

PRINT  4SO.LANROA 

300206 

333232 

ISC  FORMAT  (IX, UH  Xm  ,F9.3,-h  ZN,.F9,3.<.H 

.  370  FGRN4T  (1X.11H  3SLI1  >  .F9.3...H  IN..F9.3...M 

CM.) 

CM.) 

000203 

1  3  2  2  n  a 

ISO  FORMAT  (IX, AIM  3ELTAI.  »  .PR.3,-H  IN.,f9.J.*M 
3R0  FORMAT  (1X.UM  (2ELTA2  *  .P9.3.,M  IN. .F9.J. -M 

c*.  > 

u*-l 

100210 

3332*1 

-00  F09H4T  C IX  •  UH  3CLT43  •  XN..F9.J9** 

-»13  FOP  MAT  tiX.il,**  M12_  a  .F9.31 

CN.  ) 

000212 
_  03322.3 

-20  FORMAT  (IX, UH  H32  »-,F9.J) 

-30  FORMAT  (IX,  U  H  _  REX  «  ^FR.01 

— E Eg 

wKmn 

A»0  FORMAT!  IX,  UMRE0SL7A  *  ,F9.0> 

■.SO  ffC.R HAT  (IX,  11HREQELT A!  .  .FR.01 

000216 

233217 

-SO  FORMAT  ( 1 X,  UhREOELT  A2  ■  ,F9,0» 

-70  FORMAT  (1X.11MH20  TE.1P  *  .F9.3) 

000213 

033219 

-40  for mat  (IX,  UH  T  HAUL  «  ,F9.3> 

-RB  FOR  HA  T  ( 1.X, UH  NU  .  *  .FR.S1 

0  00 220 

3  03  22* 

£00  FORMAT  (IX, UH  UFS  ■  ,FR.2,4H  rTPSEC  ,F9. 

.£10  FORHATUX.* _ UAVGr  *,(3.3^.  FT/SEC*  .79  .3.  «  , 

3  *  3M  N/3cC  l 
«/SiCM 

000222 
_  333223 

520  FORHAT  (IX,  UH  HYOIA  «  ,F9.3,*M  FT.,F9,3,«H 

.  330  FORHAT  (IX.  1AH  REHYQIA  »  . F9.il 

*•  > 

000  224 

0  33  22* 

£1.0  FORMAT!  IX  ,  UH  REHALF  »  ,F9.0) 

-558  FORMAT  ( IX. iAH  .  ...  CF  •  .F9.A1 

300226 
333 22T 

540  FORMAT! IX, A1H  UTAUFU  «  ,F9.*I 
_ 570  formatiix.t.  canc  in  •  ppm»> 

000223 
_ 033229 

540  F0RHAmX,»P0LTF|.3N  «  •,F9.2,»  CC/HIN») 

—590  FORMAT (iXifROUT-UIN  *  •.FR.A.p  FT/SEGM 

300  230 

333231 

500  FOR HAT (IX,  *  SLOT  *  • ,F9«M, •  IN»* 

=10  FORHAT ! IX ■ *  OFLflN  *  *.FP.3.*  3PH»> 

300222 

3  33  2  5*? 

420  FORHATIIX,*  Q»OL»  •,F9.3,»  FTJ/SESM 

£30  format ( IX ,• —  SHALL.  ■  •  , F9,  2 , •  PPM*) 

300234 
-  .  _  .  133239 

eHO  FORMAT  1 IX , •  CA/CI  «  • ,F9.»t 

080  236 

450  FQRMATdX.*  LANBOA  *  •.F9..1 

_ 108237 

1 
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PRINT  JJO.IOEMr 

- -  PRINT. -663 - 

eta  f «n*ri//»  nm.i  teen.)  uift/si  Jt«/si  u/ufs  t/ 

’nn  nt  ret  i.-Y  nr.  arv  »>»»i  ti  UIO  rtru  .. 

Z  ZUN.IVI 

- - - aa  «7J-  ui-n — - - 

670  PRINT  640.  T  1 1  ..INI  ,»H«I>  , U  tl  t.INI .until  ,  UR-  710  d  .HH»  ,  T0ELTA1  ( I  ,» 

- 1 - »u  »aiU»Z.l.UWUZ6tl  II  .a£TltI.NHl.ja£LT06I.»ita»^4iTA 

Z  I).Z(I) 

aJLC  £aUMLtXlX«JJCLIU3X—  -  _ _ _  — -  - - - - — 

PRINT  600 

694-  70S NOTA/O*-* - T  «SUJ - XSZtLti - UOUf  S - tl/71  URmS..  — 

tYPLUS  OEL  UTiU  Y/H4CFH-  LOStYPt.US>  •/ 1 

- ao-.Tio_i«v.Ji - , - — 

PRINT  73fl.r(I,«NI  ,T0£LTA(t.NNI,U*»Tt3<t,«m.URiri07«I,rtm  ,  U 

-  —  1 _ »i  inir.wi.mimr.Mi  .aiiutJiiiT.Mi.nMitJ.mi^m  ua  -LI-4UU - 

TOO  F0RH4TtlX.4F13.31 


PRINT  7Z0 

—  gna«iT ujj. _ y tun. 

00  730  Jal.JCONC 

—  r-tn  jptmt  Tua.rr.LA.HK 

740  PORH4TUX.5FIO.31 

_ SO- 10—770. _ 

750  flNNtLU.HMlaltaLAS 

- sa.to_7*a - 

760  NNN (LL » HHI *<N4LL 

_ _ 5O_TjC_Za0 _ 

770  MNNtU..NNI#N 

- IFtU-^ca^U-.OR^- 

780  CONTINUE 

_ _ cou— eunputuwni - 

CN*0. 

_ ONtllO^ _ 

QUO. 

_ ci»<u - 

RENINO  1 

_ „sa_Ta.iaa _ 

750  CONTINUE 

:r\a.<»9 _ 

900  CONTINUE 

_  £(U ...  - - - 


300236 

-aocjja- 

3002.3 


300  24- 
-400  7-9- 
1002-5 


J 00  2.4 
■000  2-9- 
000250 


0002:9 

-400255- 

103256 


300262 
-440263- 
300  264 


000264 
-40026.9- 
000  270 


ici. oi.rc, og.alu  a.  w.zztasssj. OArAPAci.oi.ro, os.  1200.  no.  mzio  ,  ,s  ooooac 
?-*»<»  &^uim>i*»A6Ma^aAtA«*6t .au.r;>a>,  a»*^m«<iAS»ca»Lia*uu.  w~ ~ ~ 


3cz. *•**,!•  3, v.pcir, master. Lia*ARY-oi»eCTaaY.C2.*»»*.i>  l.xj  ioqooo 

_ niiii _ oooaoi  - 

overlay <e>  oaoooz 

_ .•.n.iwm/i»a»yM.i.{.»ain  -*»i  -  .  nnilJAl- 

OOMHOM/IOVeR/IENa.ZPASS  300004 

_ i£!UU0 _ aaoAot- 

ipass*o  aoaaa6 

_ ia_  oomt  inuc _ oao-a&z- 

CALL  UPOVcR  (56,4MTAG1>  000  300 


003010 

-130011 


100012 


CALL  jrav£*m,4«<TASl) 

.flfl-ia  i,o _ 


20 


:ont:*u£ 

■MO _ 
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_5UBBOmi2i£-£UNFLflT  112) . 


COMMON  YLAM(31,Z1  .TCCU1.21 

.  caniMN_cx.taiui2«cu»£- 


CQNHCN  0ELrA,aELm,GELTA2,MY0IA.HlZ,HJ2,:GENT,LL.NC.MCH,;<Pt.aT.  BE 
_ 10  EL.  7  A  <A£G£L.TAL,  ft£3£J.rA2.,££HlQZA,A£X,U£S..J(ZJI_ 


oaaaoz 

.0  30.3  33 


ooaac<. 

~3f.  flBu 


cannon  zfntiioi  ,.ycas£<i2i  ,nxcnab<i2i  .nychmiizi  .ntchapiizi  ,  oao3is 

.Uflj.rjju _ azmaia 


COMMON  XAX(9,121  ,YAX(9,12I.TITL£(9,121,i;MAaACt9.201  003006 

ITIQI40.21 ,U«ATI07(Jl,2j ,Y ( 31,21 ,rOELTA(31,21 .YOELTA1 ( 31, 2) ,  YOELTA  300007 
.  22  (31,21  ,  YP.LUSL  ( 31x2J _ 000007 


COMMON  CCA  <3 1,2  I 
1IMENSI0N  -PACT  (123  ..KGBIO  1 121  .  XLiagL  1 121  .  <SVH  <  1  2  ■  21  .  LggTGX  (121  ■ 


000006 

;  "rnirinq 

1107(12)  .NO (281 .  NUEELUZ)  ,NLwGX(121  ,NLOGY(12)  (LOOK  12>  ,  NL0G2(  12)  300  009 

_2,NSTM<12.21  ,NXt  M2)  ,NYT  1121  .PX ( 121 . PY 1 12 1 ./ axl < 1 21 ■  iHi  t  121  .  xB2  1 12 _ 300  039- 


3) ,XC1 (201 .XSN0I12) .XGSID (12) ,XNUM<28) ,X0F(12),  XUNIT <12 ) , 7 AXL (1Z ) , 


GIMSNSION  ?1  ((.0.11,21,31  (<,0,11. 2)  ,P(20I  ,  YNUM  (20)  .  IPZ  (12 1 ,12  (12,2) 
(lATA  (NLflGX(I)  ■I»i.  iil/*»a.  t-..»n/ _ 


300009 


CATA  (NL0G1  (II ,  1*1,111/11*3/ 
.GATA.  (MLGGTiZJal9li.nl  211*  XL. 


000010 

■moan 


000012 
■mm  Y 


GAT  A  <NLGG2(I1 .1*1. 111/11*3/ 

JJ.AT  A.  JLGAiax  (IJ.n*lilU  /.i»0 .  loPiQ  L. 

OATA  (LGaiGTdl  ,1*1,111  711*0/ 
JJAtA_iKG5I01IlaX*l  xlU/JJLIl/. 


30001* 

300019 


GATA  (NS  YM (1,11 .1*1,111 /Il*67000333a000a0033/ 

_ GATA  (NSYM(X,2  1,1*1  ,lU/.iil 

13000000000000009/ 

SAT  A  KSZfi  T^»07,n/Xt.  1*4.03/- 


303316 

n n ft  "i  7 


OATA  (XGSIDIIl  .1*1, 111/6*0,3.  ,.*.5/ 
3At*_(X.G.?I9.<  I.Ut23UUi(lli2aJ2 _ 


aocois 

_30JL31i_ 

300019 

_101323_ 


GATA  <FACT(tl, 1*1, 111/11*1./ 

_ flAT A._(KLA3EL( II  xl*l,ll)  211*1/, _ 

GATA  (NLA3EUI1  ,1*1,111/11*23/ 

_ OAIA.  (NIT  <11, 1*1,111/4711.,  13, 11*0 

GATA  (NYT (II ,1*1,111 / 11* 1 1/ 

_ 3AJU_<PJS.<  11x1*1, 11L2>*  3 ,a-2 .5. 3. .2. .**3./ 

OATA  (PT(I), 1*1, 111/11*2, 5/ 


003C21 
3  BO  0  22 


000023 
303  32L 


000025 
.333  026 


000027 

JL0332JL 


_ OATA  (XAXL.II1  ,.I*l.ill/<,*5,,6,ai,x/aJa*iix/. _ 

OATA  (YAXLII1, 1*1,111/11*5./ 

_ OATA  (XUNITCl  ,1*1. 111/5,  •  l.,3,«10x.,3«.xl,  ,3  . x* 

GATA  YUNIT  (11/0. 250/,  (YUNIT  (11,1*2,  111/6*1. ,5. 

_ 3AiA_iXG£<.Ux I-»ixiiL/-**a.  .2.  .6*0./ _ 


000029 

_ 300.03.0— 

300  <3  31 

L1*Z _ 3.80  332  . 

.*0., 23. .1. ,.30.5./  300033 

_ m0.3k_ 


OATA  (Y0F(I1. 1*1. 111/11*3,/ 

_ OATA  .(XEHDII)  ,Z*UU  1./2.V  I.,  i.2x»- 1,.- 2.5  , 

GATA  (YSNO(I) , I«l, 111/2*7, ,*21. 5 .2*7. , -21. 5 ,2*7 ., -21 .3 . 7 ., -12./ 
_ J3ATA  (*81(11  ,1*1,111  /imi«z _ 


OATA  (X32CI  ,1*1,111/11*0,/ 


303  335 


303337 

-iiaojua.- 


OATA  (Y82II),  1*1. ill/U«4. 3/ 
.OATA  «ST*MIt  il  ,  Z*l,  U|  /t  l*2/._ 


303039 
inn  f.,a 


0000.1 

JlJlOJla,2_ 


OATA  (KSYM(I,21  ,I*liH>/J*0,2,2»0/  000043 

_  OAT.  A  .(XC1  (Z  ).*l*.ll5  I  /.  3*1,273/.,  ( XCL(I)  .1*6.131/5*3.779/.  (XCim.IH  3000.. 
11,131/3*5.4  73/,  (XCi(  11,1*16.201/5*7.375/  3000-4 

_ 0  A  t  A—LY.0 1  JXJ-xi*-!  J  1 6 , 3J  2.41.1^252 _ 330  3.3  ■ 

OATA  (YC1(I), 1*2*17, 31/4*1,00/ 

_ OAT  A  (YCKIl  .1*3. 14. 31/ 4*0. 78/ _ 


OATA  (YCKIl  ,I*<>. 19, 31/4*0. 53/ 
_3AT  A _ (TCI ( 1 1  x_I*5j 23  1 iXO>±9-,ZS.t- 


3  00  0-6 
JBOOac. 

3000-4 

_3.00.Lx3 
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OAT  A  <NO«I>  .I»t. 201/0,  .»»I.3»0.5»3.1.3,  1,2*5/ 
.  Qiri-UNUIUUv^t.i 
II  »tl. 151/5*5. 675/.  IXNUnUl  ,1*16,201/3*6. 675/ 


oaoaso 
.  coco  si 
aaoa=i 


■Ptllllggfll - -  — - 

-  -  _ ...  ■  _3amxz 

aoaa?J 

Tnn  n  =1- 

pui«ni) 

000055 

?  _ 

M.iflSi 

P<61»«£X 

000057 

_  ?(7)*flf4yQri  . -  - 

PCil«P60£LTA 

300050 

Pt31a*rn£l_TJ.:L  -  ...  __  - 

P(ia>M£0£LTA2 

300061 

pu  11*051,  IX 

_  _  -- 

P«12>*0£LTA1 

000053 

.  Pd.tlafltt  TAP  . . 

®U*I«H1 2 

000063 

UlflflAA 

p<i6i«c: 

300  067 

.  ?  f  1 7 1  mr 

.inn  /» i  a 

■•lilac* 

300069 

3C1A»*III 

mirn 

■  (20J  *5 

300071 

2Q  l*  ;s«.^P(_aT  _  _  _ 

TQn  n?? 

*C««CAS£(t> 

300073 

.  30  IS  Kat.JiC . 

- 

ISZ*IZ(I..il 

300075 

30  t  0  .1*1.132 

.  ..  T  3-5  3- Si 

piij.i..h*u<j.yi 

000077 

OKJ.l.KIlYU.m 

.  _  ancor? 

Pt(J.2,.H*YQ£LrA<J,*l 

030079 

.  21(J.2.HI*URATIOIJ.ftl  ..  . 

*  .1311  a  A3 

Pl<  J,3.nl*tOSLTAl  U«-H 

303  041 

.  3l(J.3..«l*uRAT13c  J.(*I 

100  34  2 

■If  J.6..H  -Y0ELTA2  (J,9» 

300043 

'31U.4..Y1  »U«ATI0«J.HI  _  . 

^aanAp* 

PUJ.I,.1»*«6YL(J.«» 

000035 

01  f  J.5..Y1  *URA7la«  J.PI 

P1U««.1I  *URATI0 ( J,  W> 

030067 

31tJ.6..H»ltT4lJ.)H 

1Q2Ja« 

P1U,7,M1*YPLUSL(J,M1 

300049 

_  .  'Jl(J.7,*t*UPLUS(Jt',t  -- 

■maaaa 

Pll  J.4."l ■  Y0£L  TA ( J •  HI 

300091 

31IJ.4,MI*(3UUTAU(  J.MI 

aoojiiz... 

»UJ.9.11*Ui»ATIO«J.ll 

033  053 

...  -  -  PltJ.9.21*URATI07«J.l) 

31tJ,9,-*l*YO£LTAIJ,l> 

000095 

nu.u..inccK(j.«i 

31U.13,rtl*YCC(J,«l 

330097 

-  .  Pl(J.lX.«l*CCmj.«) . 

_  _  aflflnfld 

11U.lt,MI*YLAH<JfIH 

000059 

to  co.ntinus  .  . 

_ aaa  tea 

29 


39 


30  23  ;«t,.NCH 
YNUmI)»YCl(Il 
*PIT£(l,73l  HP'.QT 
CO  60  1*1. .'•PLOT _ 

•yc»*ca;i(ii 

_ _ 30  30J«*l,HC _ 

:pi  (Hi«rzd,H) 


.  «NTCHtP«NTCM44(X>  . 

XC*,5*  lxai«;i-*GT*XHTCMA*) 
r  C*t.625  _  .  _  _ 


300  101 
-1103X02. 

330 133 
-0X0X1  >* 


030  105 
-000X36- 


000  137 

JLiaxoa,, 


ooaios 

-00313.0. 
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nrite ti.aoi  ricTin.xGsioti) .rosin ti> 


again 


« 

-AAizzmu  <LA&c.i-iL)-*JiL*azi-Liiuazisunmt*La&JUii  Mica., 

t  «i)  .Ntasztn  .uSRiaxti)  .lgriby  tio 

■  ■JflflUZ- 
aoaii2 

•i  n  u  n  .t 

30  -3  3*1 tNCH 

300 11- 

4iltTT-t.if.4fli  tci 

3RITE  tl.100)  (CHARAC(JK.J) .3**1. 31 ,NCH4«(J) 

000  116 

-ft] 

9  tJJ  .MCT ->i  n  1-11  _  _ 

^11**7 

sa 

CONTINUE 

4*TT*ll.t7«l  tc ASStTl  .NVTIT1  TIT). HUT 

aaa 113 

-  -flrt  *_ia_ 

4RITS  1 1*130)  PXCI)  .Ft  (I)  .  XAXU  I) ,  r  AXU  ( 1 1 ,  XUNIT  1 1>  ,  TUNIT 1 1 )  .  <CF( 

000  120 

i 

n.yoiPftl  _  . 

...  ‘•flfl 

4*in  (liUOt  XcNOtl),  VcXOtXJ  •  XStfX)  «X82IX)  fYdl(X>  tVS2C>  «xC'VC 
MflXT?  tL»n  .HXCuAS  ttl  _ 

aac i2i 

-  flflflLM 

9SITS  (1.100)  (TAX  ((.,!>  ,L*1.9),NrCNAR(I) 

3»rrs  t  moo)  itttl- (L.  :i  .l*i.  4i  .htcmis  tti 

330123 

1  7U 

3  SITE  (1.130)  (KSTNtl.N)  ,.r»l,rtC> 

000  125 

*»rr*  (TTtt.Li  .1*1. nci 

3flilL2fc 

3RITS (1.130)  (IPZ (L) . U*1 ,«C) 

000 IZ7 

(TPHTfLI.l_al.lfll  _  _ 

_.3flftl2A 

00  60  N*1,NC 

300  129 

_ isisizxnrj _ _ _ 

,100.136 

asiTEa.lFNT)  ifiu.i.m  ,j«i,iszj 
tFocsv^c.Mj  ,£Q.ai  go  ro  &a 


_ xhite  ( m-o  i ..  Nsjtfi  <  jjtsXNUn. 

60  CONTINUE 

_y_E!»  N*  r  1 9  is  j _ 

SO  FO^fiATCJFlD.il 

_  «a  ^aSNiro.Fio.ii _ 

130  FOSNAT (9A4.IJ) 

.115  FORN«T(F12.5,.Fl0.5jt5».-_, _ 

1ZC  ?09H»T(3I3.F3.3) 

t3a_e.0RN4Tt.4F.ia.4i) _ 

144  FQ9H6TU5.I3.F3.3I 

.I5.a  FORMTU&H) _ 

160  FQ9N4T(10A6> 

_ ». IT URN _ 
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000131 

000132 


080133 

:aa,  13.14. 


000 133 

..inaxia.. 


000 137 

aaaiia. 


000  139 

oaaito 


ooom 

-oaana. 
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oouble  fUEcisxoN  a. a 

ati.  x  u.»3-t-.a-a  «-t  i .  lafAusuz, 

DIMENSION  X3Ai‘A(31)  ,rOATA(3U 


2 


— aa  - 

*P?*KPl»l 
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OUUWmUiUltUUU  —  -  _ 

30  10  I»1.N  000006 

- CF»K^a-tl~M.  tO-ta - 100-04-9 

rt<I,J|iH(X,J)-rt<I,iC)>tt!K,J>  000  010 

ia-ca!irxj*u£ —  - — - ooaoit- 

;rWN  3QQQ1Z 

_ ilia -  ■  -  .  ..  —  .  n.  _ 04004  v. 
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■agsgautw8D-wiw>t» 


Mimm  at 


■-.aaaoi-. 


OVERLAY  (94» 

— — COMMO  tUiC  VER  XI  £  IIM  P  A$S- 
CONHON/G/PX ,PY 


.  COttMOMA^Y  H/ lACM^SX*  tf.  >- 


CONMON/SIN/XLI.XKIZ*  .XLIN2(1Z> .YLIM1(12>  .YLIM2(12J  .KNL 

..r.MHOit/Lnii^Lar.x.YirL.rAx.MTiTMAa ,  nma»»taajaita»s<»tax— 


3  03002 

-444103- 


000034 

-444045 


3:3306 

-344403- 


SIMSNSION  ISZ<16J.I?Z(16) ,<SYM(lil  303008 

- GinfrXSMIA-XAX4*AaPtAX4»W4ttLS4»UlEHW«A - 444444- 

C...NPL0T  IS  THg  NUMBER  OP  SEPARATE  PLOTS  000010 

C...XGRia,YGRIO  ARE  THg  5310  WIDTHS  303012 

.o^jcLAaai^XMMCAca-Miciiiausamcic  «iuia at-iAmam- una  jrrsas 344443- 


C...CN  THE  PLOT 

-o— -'u-3a£L-  -is.  JMS-Mwaga-or-  labelej— variables— — 

3...IF  <5310  *  !■  A  5310  HILL  BE  PLOTTED 
-C— ..SLOCX  t.S  THg-  mcagsi— PCWga-O*— THi  uaC-X-AXIS- 
C...NLCGI  IS  TM£  LOWEST  POWER  3F  THE  L35  X-AXIS 


C...M,aG2  IS  THE  LOWEST  POWE3  OP  THE  LCG  Y-AXIS 
-C-.— LG3iax.i43iaY— maiCAJS— wheimer-lcg.  colas.  Age  ra  a  aax-wM- 
C...nCASE*MUMBER  OP  FIGURES  0*  THE  PLOT 
4-— aXT.aXTt-A.RE_  NUMBER-  ap_ TICX-MAA<S- 


C...hGT«mEIShT  3P  THg  CHARACTERS  TO  BE  PRINTED  IN  TITLE 

-C...PX.PT  are.  aEGi.sm.xr.  caoRniMAiEs  of  thf.axp.s _ 

C...XAXL. TAXL  ARE  THE  AXIS  LENGTHS  IN  INCHES 
C  .  .UUNII . r UNIT  ARE. THE- 5CALE0— LENGTHE_a£_tHE_AXEE _ 


:...xoF,raP  are  the  jfpssts 

C...XENQ.  TEND.  ARE.  THE  CaaRaZNATES-aF-ENPLE- 


Ca.aX31aX32aY91.YB2  ARE  THE  3OR0ER  LENGTHS  IN  INCHES 
0 . -^XCxYC  ARE _ I_H£_  3EGZNNLSG_C.0GRQIiSATr  E  Of  THE  IU,L£„ 


C...XAX  IS  THE  TITLE  OP  THE  X-AXIS 
.  C  .  —N  XCHA.R.  IS  THE- NUMBER.  0£_  CHARACTERS-  IX-AAX- 


C...YAX  IS  The  TITLE  OP  THE  Y-AXIS 


aoaoiia 

444415- 


300016 

444442- 


000018 

POOCH- 


000020 
■  344421 


000022 
344  425. 


303  02a 
tih  n?c 


303026 
344  027 


000028 

-04042.9 


300  CTO 
.300.0.11 


330032 

_J.40.U1 


3  00  034 


C-^.RYCHAR_  IS.  THE.  NUMBER-OP  CHARACTERS-IN-1AX _ 

C... TITLE  15  THE  TITLE  OP  THE  PLOT 

.. C  ...nTChar  is  thp  .iuhr-o  ag  CHARACTERS  th  thf  ttti  ^ _ 

C...<SYH»1  INDICATES  THAT  SYH80LS  ONQ  uINcS  ARE  TO  3t  PL37TE3 

C...KSYM«2  INDICATES  THAT  ONLY  SYMBOLS.  ARE_ia_3E_PU3.LTEa _ 

C...ISZ  IS  THE  NUMBER  GP  POINTS  TO  3E  READ 

C.a.IPZ  IS  the  NUMBER.  OP  POINTS-TO  3£  PLOTTED _ 

C...IFHT  IS  THE  FORMAT  OF  THE  X,Y  OATA  TO  BE  READ 

_ CAU^.?U3.TSU.a.UJl> _ _ 

REAO(1.2fl>  NPL3T 

...  OQ  110  <LN»l.’.NPLO£ _ 

<ML*<LM 

.  _ REAOd.lOI  FACTa.XGSIfljY.GRU _ 

10  P0RH4T(3P10.3) 

_ ULA-PACTQRtPACTl _ 

RSAO(l.ZO)  <LA BEL  a  NLA EEL  a  <GRI0a  NLOGX.NLOGl ■ NLOGY . NL3G2 , LGRIOX 

_ 1  5RI0Y _ _ 

20  PQRHAT (913) 

-  .  IP(<LAaSL.Ea.  1)  CALL-LA 8£L.(HLA36U _ 

REA0U.301  hCASE.NXT.NYT.hGT 

.  ...30.  FORMAT (31.3, F.5. 31 _ 

REAC(i.-O)  °X ,PY,X AXL  a TAXL. XUNIT.YUN I T,X OF, YOF 

REAOU.-O)  Xa NO .  YgyO .  X31lX8Z,YB L.I..Y.S2 a  F  Cj-YC _ 

40  FORMAT (8FL0.il 

_R £40(1,50*.  (X  A X  (L  I  ,L*l,9H;X.CHAR _ 


_ 300  5.35- 

303  036 
to  a  n-Y7 

300038 

_ _ 344  UB 

3000-0 

_ 0444—3- 

303042 

i an q.3 

333  0  — 

_ loaa-j 

333  046 

_ 1044-2. 

3330- 5 
1P01-A 
L  300050 

_ 040  010 

30SC51 

_ 130052 

303053 

30305- 

000055 

_ 343154 

300  357 
_ 190.155 
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50 


AEAOC1,50>  CYAXCU  >01.91  .NVCMA4 


F0.9NATC9  66.161 
_ 9M0ASE- 


AEAD11.60I  (KSTHCU  ,L*t.N» 

JSAaauii-. imiu  *uaauL  - 


9EAQC1.601  (IPZCU  ,k*t.«* 

_ 60_r  .04*61.11613) - 

PEAOC1.70I  CIFXTCU  .k*1.10) 
_ ZLcOWArtlllAOJ. 


XkINl<KkNI*XOF 
.  <i  rx?  (<i  ii»»uNiT.t3f..aafli 
yunickdimyoF 


_r  uifttiKui  i  »iun  rt*t.cr i. 

:au.  PLar(X82.rai.5> 

CALL  PiaTCXgl.791 -PI 

call  PuoTixai.raz.2) 

_CAki_eioT.txataaa.tZi— 


4.0.04. 


CALL  PkOT(X32.YSi,2) 

frntsarn.go.il  call  eatn i yaxl. yaxl ■  xsarn. reatni. 
:;u  cha«ixc.yc.o..h«t. title. NrexAai 

tPiMLQsg.Mg.a  .an.  .4Lasr.Mg.ai  pall  i 


X , NX CHAP, NCOS  1 ) 
..  XSF.fX  AXL/XUNLX— 


asArrx.M.  ...  it. 


rSP»»AXk/XUNIT 

_ iAkk  »LOTtPXtPTt-JJ _ 

_  CAkk  OFFSTCXOF.YOF.ll 
_ iAkt.  SCAkf.«XSP,TSP_,U_ 


SO 


90 


[FtNkaex.Nc.ai  so  to  so 

CAkk  AXIS  « 0.  ,  a  .  .XAX.N XCHAR.X AXkjJtXT.  1*40 _ 

tP<KkOCY.NC.a>  SO  TO  90 

,5Akk  AXIS  CO.  t-0.  .  YAXtttYCHARtJAXkOtXT .0, .3.1% IS 9/2.  t 

30  iaa  j*i,n 


..1»IS2U>_ 


;z*i?zcji 

.—SEACCl.IFNTI 

9EA0C1.IF9T1 


CU(Iull,I>t..MI._ 
CVCI.JI  ,1*1, Nl 


_<SYX9»<SYX<  J) _ 

IFIKSYM8.Nt.ai  CALL  5YN8CXCIZ.J) 
.IFJICS  JX3.  60  ..2  l_SQ_LQ_l  u_ 


CAkk  PkOTCUCl.jl.Vll.JI  ,31 

.30  tea  i*i, iz 


i 


IFCUCl.JI .kT.XkINICKkfO  .CP.  UCI.JI.C7 .XkIX2(XkX> >  30 

too  .  .  .  _  . _ 

IFCYIt.JI  .kT.YkI.XKXk.XI  .0*.  YCI.JI  .ST.YklXZCKk.XI  I  JO 

_ioa_ 


.100 


__j.ia 


CAkk  PkOT(UII.JI,VCI,JI,2) 
-CONTINUE. 


;ikk  ASSET 
CAkk.  PkOTCXENO.  YENQ.^ 
CAkk  ASSET 

CAU_.£.XPJ_r  CXENOjutEMIU 


•NO 


000059 

'inn  it.i 
000061 

0000*3 
tlnn  Hi. 

300065 

jjxajiai. 

300 C*T 

300069 
inn  in 
000071 
onni77 
000073 
304076 
300075 
3,00377. 
000077 
-3104ZA 
300079 
^nnnxn 
300060 
-1440A1- 
3000)2 
■400063 
aaao6- 
-340035 
300066 
JL0J.44Z.  • 

3000)6 
_1004a3_ 
330390 
■143  331 
300092 
300093 
00009- 
..33.34.95 
303096 
033  331 
303396 
.103.039 
330 133 
-.30413  0 
330131 
TOO  1  Of 
330  102 
.104103 
30313- 
-.314143. 
330  106 
080107 
303130 


- - wwmcsa-uwws - 

COW*ON/G/PX*PT 

OOOC02 

LX*XAXL/<2.*XS»IO>*.0W9 

000004 

LY*YAXL7  I2.»YMI3I*,S<H* 

000006 

YAaPY*YAXL 

000003 

?*f»x5sro.tr-j , --«•  «•?  t,7  if 

CALL  ft.OTIXA.AY  ,31 

000010 

<•-*«,  V  ai>at-t*4.n-VAT^^ 

30  IFlYSaiO.LT.l.i-S)  SO  TO  40 

000012 

CALL  plot  IPX. r A. 2 j 

oaoaio 

40  P»PX 

000013 

v««r O.LT* l«£-fF  fn  Tfl- 

00  SO  I«1,LX 

300013 

_  *mQ+*rratry  .  - 

call  plot(p,pt,3> 

030020 

ffAIrL  9l.aTIP.VA-?>  _ 

IF1I.E0.LX>  SO  TO  eO 

000022 

call  PLoriP.rA.3i 

00002L 

sill.  pi  ar(P.»4.»Y 

30  CONTINUE 

000026 

&fl  CALL  PL0UPX.AT.3J 

_  nan  .a  72 

70  0«»T 

000023 

*rjTG»ra.LT.l.--«l  SO  TO  40 

-  _  aacfl2^_ 

30  30  ;«L.LY 

000030 

o«a*»c»io 

iniifni 

CALL  PLOT (PX.0.33 

•  CG0  0  32 

CALL  PLOT (XA.O.A) 

iflflcn 

IFd.SO.LYI  T£  TURN 

000034 

c*o*r CP :a 

.maa’5 

CALL  PLOT (XA .0 .31 

000036 

...  CALL  PLOT  IPX. 0. 21 

'inn  nr 

30  CONTINUE 

300  013 

30  CONTINUE 

aaofl^a 

AETU2N 

0000.0 

£N0  ...... 

3JUU1-.L 

[  I 

t  1  281 


_ gg<gmigfciggn_  l  f  <tt ms _ 

;...xc,rc  the  beginning  scquoiujtes  zr  the  rfjatinus 

OQG  002 

i...CHi«AC  IS  THE  ALFM48ETIC  CHJRACTEaS 
^....iCm*3  rs  rut  miimmfs  n»  i:mj»sr:n  r*i 

000034 

annii«« 

C...rf*U/E  IS  THE  ViLUC  or  THE  rf»aiASLj 

uaace 

C * . .NO  IS  TH£  NUMBS*  OF  OICITS  to  THC  RIGHT  Of  THi  Q2CINAL  POINT 

CMAftJ.fi  ifti 

000003 

oo  *o  i*i  *n 

9*&o<i.iai  jtcf  re.  <num-jhium 

oooata 

motile 

to  FORMAT < UP  10, 5  > 

_  _  R£JQ(1*2Q)  _(  CMAftAfi  1J1  >~Jsl  «  31  »  .'4fiMAft_ 

000012 

20  FMM4TI94S.I8) 

JfiDll.XO)  rfJLUF.MGT.MQ 

000014 

30  FOR HAT (F 12, 5«FlO  «5* 131 

,  -ALL.  CHAftltC.ffi.fl- .HCTaCHAflAC-.MCMAftl  _ 

000  Ole 

40  CALL  NUMPVT  (XMUM,YWM,0#tMCT*  VALUC'NO) 

RETURN  _  ___  __ 

000019 

zNO 

aoau2o 

50000 2 


suaaauTTti*  tYnanm  tx.ji 


..NSYM  IS  IHS  SYM30L 

^.jtsiMt.is-iHS.  mTaavflu.^r-MHi;R-S3ta«aiA-*af.  to  u  »uaa 

..M  IS  THg  H*ISHT  Of  TH6  SYH90US 


_u£.iBoa/  si  kuso.w£5- 


COBNON/SZN/XLXHl  ( 12  >  •  Xl.I'l?  ( 12 )  •  YI.XM1  (12)  .  fLIfl2(  121  ,X«C 

arii‘1  m  -  i  n  i  >K»».ftTiif  .n  _ -  -  .  - 


10  aoi»4T(A5.I3.AJ. 31 

-10-20  .I»l»I2*XSXNr - 

lA(UCI.J)  .LT.  XLI.Bi  (KMU 

. . . iKuxx.>jtu.r-rum.mwu- 


CAU.  VeCTOA(U<I.JI .Y1I.JI ,l.l.-l,NSYBI 
-23LXQM  tI0Ul£ _ 


.30.  li<I. JI.5T.XtlB2Um.il  50  TO  20 

.c«.  ■  /  u  ■  a^iuumium  .50-10-20 


aeruoH 

- *M - 


il  00.501 
300004 
300005 
300306 
3,300, OX. 
300004 
300.00  » 
300010 
300011. 
000012 
mao  :o 
oooou 
300X1.5 
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ag^gnnrwryn  i  rtrrur. _ _ _ 

1«C31 

CONHON/LCS/NLOGY,YAXL,YAX,.'IYCHA«,.XLOS2,LOaiOX,LS«IOY 

300002 

_  arntMSiflu  r rfUXHL111 

QINCNSION  XLOGU03  ,IX<1»  ,XAX(9t 

00000- 

X-»AX1S_PQM-B<  rtr  IQ 

PX4X»*X 

000006 

PYAXmPY 

tMNcoGx.za.ai  so  ;o  ro 

000006 

XAXPvPY*XAXL 

CALL  PLOT  IPX. PY,  3) 

000010 

CALL.  PLOT  13CAXP.AY.2I 

000012 

pxs«*x 

PXaPX-.l 

00001- 

?Y*PY-.2§ 

aaiifli « 

XlaXfcXC*l.£-4 

000016 

>ILaG*NL,aGf*?*^2 

a.tflfli  r 

XNLOGX-NLOGX-NLOGA 

000016 

x2*xiyxNLaGx  _  *  _ 

PXZ«PX*.2S 

000020 

PY2wPY*.fl7  .  _  _ 

iPiNLaGi.-r.ia)  .'it*t 

000022 

rPfNL0Gl.S2.10l  NL«2 

nnn  n?i 

CALL  CHARfPX,PY,0...10,2HL0.21 

00002L 

3 ALL  CH»«<PX2.PY2.a...117.NLQG.NLI 

.inn  up* 

NC3GX1«NL3G|.»1 

300025 

00  la  TsNLOCXl  ■  ,‘iLQCX 

inn  a 27 

tru.LT.tai  ic«t 

300  026 

. .  .  ird.cs.iai  .il«2  . 

0flfl02q 

000030 

PX»PX*XZ 

inn  n  ti 

*X2*PX»«2S 

300  0  32 

CALL  CHAP  IPX.PY.a...ta.2Min.2> 

inn  nil 

10  CALL  CHAR<PX2,PYZ,a.,.07,IX,NLJ 

03003- 

0»a w  x-axt*  rilX_»APKS  UNLZSS  .GRlfl  rs  ro  3.  oraum 

ananas 

PY*PY».25  1 

303036 

rpiLCSiQx.za.ii  so  to  ta 

300  037 

PYt»PY-.07S 

000  336 

?Y3»PY-.l 

000033 

00  $3  I*NL0GX1 ,  NLOCX 

0003-0 

IQ  «£  Jrl.lQ 

anca-t 

XJ*JY1 

0000.2 

CALL  PLOTfPXX.PY.il 

nna  i«.i 

IPfJ.Ns.l  , ANQ •  J.nE. 13)  SO  TO  20 

3QQQ~« 

_  CALL  PLOTfPXX.PYi.pl 

SO  TO  30 

0000.6 

20  CALL  PL  OT  t  P  XX .PY 1 1 

ma  c-,7 

30  XL3GI J»«4LOSIO(XJI 

aoso-4 

rPtj.PO.isi  r.a  rn  so 

jaflfl-s 

LO  PXX«PX3»XLOOf JI*X2 

0000:0 

SO  PX3«PXX 

jaaa?i 

C., .LABEL  X-AXIS 

0000:2 

.  AO  PXaPXAX 

flflllflvl 

PY»PYAX 

000  03*. 

NXCMARaTAflS  f  NXOHAPI  _ 

5« 

XNXPNXCMAR 

000  0  36 

PXS».  S*fXAXL-(.  1-»XNX1)  *PS_ 

anflfl?7 

pY5*®Y-.i.J 

900054 

CALL  CM A&(PX5«PY4.i2*..l  k.XiX.MXCHl^I 
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5 


1  0  iFlNLOGY.iQ.OI  SO  TO  140 

- - YAXP*PY*YAXL _ 

CALL  PLQT(PX,PY,31 

C... PRINT  9-AXIS  POWERS  OR  19 

- irr.ua  r _ 

PYJ«PY 

.  ...  aimarm,uij  _ _ _ 

YlaYAXL»l.S-9 

- . .  . .  . 

YNL0SV«NL0Cr-NL0G2 

- YZaYtXtNLOCS - 

PX2aPx* .  22 

-  - PYZ.aPYa.47. - 

if/nlogz.lt.iqi  nl«l 
- irt>u.  nr.i  mnf  ... 

call  cma«<fx,py.o. ,.ia,2Nta.2» 

- - CALL-C*Afi.taxZ.ax4.a...a.7..NLS6.NU - 

LL0GYl*NLaS2*t 

- - 30-40— i»NLaCY.l_NLOC.X _ 

IFCI.LT.103  ,yl*1 

■ - LT.1I.CE. . ta»— :il*2 _ 

IX<U«I»2“»2 

_ PtsexAtz _ 

PY2*PY* .97 

. -_CAU_.CMARiax.PYta»..ia»2H14*Zl - 

60  CALL  CMAP«PX2,PY2.a...aT.IX.J*LJ 

.  ts  -ticx-aARKS.  aua.  fiain..:.s-.:a-.nE  j&uul 

PXaPX*..  a 

—  LKtL3Riot.£a.iL  .Ga_ra.  ua _ 

®xt«PX- . 375 

.  .  PXJaPn-.l _ _ _ 

CO  120  I»NLCGY1 . NLCGY 

_ -Id  -141.  J-1U.3 _ 

YJ»J*l 

- - 0ALL-2U3TIPX.PYY..JJ- _ 

IFU.NE.l  .A NO.  J.NE.13)  GO  *0  20 

_  - _ CALU  PL3TtPX3aPYY.2J_ _ 

iC  TO  109 

- 90 - CAU_-Pua*t?Xl.PYXi2J _ 

100  YLQGl J) aALCGlO (YJI 

IFCJ.cO.iai  GO  73  120 _ _  _ _ 

110  aYY»PYJ*YLCG<J)*Y2 

120  PYIaPYY  .  .  _ _ 

C.  •  •LABEL  Y-AXIS 

_  130  PXaPXAX _  . 

PYwPYAX 

_ _ NYCNAR*IA8S  (NYCHAR1 _ 

YNY  *4 YCHAR 

_  .  -  PXf*PX-.J  ...  _ _ _ _ 

PYO*. 5* (YAXL-«. 1»»YNYII *PY 

- CALL  CHAR  <FX3. PY8. 3.  1410  9/2.. .  U.  YA 1.  .YYCH  44 1 

140  CFtLGRIOX.NE.ll  sa  73  179 

C...0RAW  X-GRI0  _  _  _ 

PVAX«PY»YAXL 

.  ...  PY5»PX  _  _ _ _ 

PYMPY 

_ ?X7«PX _ 

PYTaPY 

.  CALL  PLOT (PT*,pYS,J|  _  _ 

4XL0S*4LCGX-YL0C1 

30  113  :«1,NXL0G  __ _ _ _ 


300050 
909911 
003062 
40A45J- 
900064 
040045- 
300006 
4444*3- 
000066 
4044*5- 
300070 
400471- 
000072 
440471- 
0  01 17- 
440478— 
300076 
440477- 
000078 
444475- 
000030 
4/13  3.1  .. 
000062 
aaaa«3- 

300  06* 
440.465- 
0  03  0  66 
OC.TCi, 7... 
000066 
■443.16  6. 
000  0  30 

100091  . 

000032 
043.393 
00009- 
inn  rati 
000096 
.040  437- 
000096 
040  0  99 
300100 
.inn  t  n  < 
000102 
.111443.. 
00010L 
.340,1,06 
004106 
-344447— 
ooaias 

.140103 
000 110 
,  840  1 :  ,1 
000  112 
-114143- 
30011* 
401145- 
000116 
-111  ILL. 
000 116 
411443. 
000 120 
-100424  . 
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XJ«1. 

_ 30-150-  jmV*. 

XJ*KJ*1. 


CAUL  PL0TIPX6. PT4, 3 1 


XJ*XJ*1. 

_ CEu.ia.ii  G0.ra.i4i_ 

PX6«PX7»x2»ALOGiO  tXJ> 

_ CALL  PLaTlPX6.PV.lX.31 

CAUL  PLOT tPX6, AYS. 2) 

_t5fl _ CONTINUE _ 

160  PX7»PX7»X2 

.1/0 _ IEIL5  RXCY.nE^JJ _ rULURN _ 

...ORAN  T-GPID 

_ 2JU  X»PX»XAXL . — 

PX»«PX 

_ SMjiBX _ 

CALL  PL0TIPX6.PVS.31 

_ NIL0GPNL0.G.I=«tLJG2 _ 

00  UO  I»1,NYLJG 


30  ISO  J*l,5 

_ 1  .  JSEJL*  L» _ 

PY6»PY7»rZ»ALOGlO tYJI 

_ CALL_EL0TtPjt6.PY6«3.1_ 

CALL  PL0rtPXAX.PY6.21 

_ tj»  r.4*u _ 

xrtj.e0.51  go  ro  100 

_ ._PY4«PY7*»2*_»L0Gia<YJl. 

CALL  PL0rtPXAX.Pt6.31 

_ _ CALL  PL0TtPX4.PY4.tl_ 

100  CONTINUE 

190  PYT.PY7 *t? _ 

RETURN 

_ ENO _ 


0001.6 
300  l^Z. 
000  1.0 
can.i.i 

3031:0 
ana  m 
300152 


00015. 
JUUUii. 
000  156 
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MYUIH-  0.071  M  KtOll  1HZ-  3377 .  H3Z-  1.670  hi  01  -  O.OZOOO  IN 


APPENDIX  D 


HYDRODYNAMIC  PARAMETER  TABULATIONS 


TABLE  23 


BOUNDARY  LAYER  PARAMETERS  FOR  WATER  INJECTED  AT  =  .6148  FT/SEC 


TABLE  25 


305 


18.00  2.817  365801  .200  .032  .019  13510  48822  7682  4644  1.654  1.805 


18.00  2.798  368596  .240  .032  .020  13706  58975  7886  4891  1.612  1.832 


BOUNDARY  LAYER  PARAMETERS  FOR  500  WPPM  POLYMER  INJECTED  AT  V.,  -  .0766  FT/SEC 


18.00  2.786  346175  .199  .042  .023  12928  45972  9627  5314 


BOUNDARY  LAYER  PARAMETERS  FOR  500  WPPM  POLYMER  INJECTED  AT  V.  =  .1532  FT/SEC 


18.00  2.760  410928  .219  .034  .021  15491  60105  9276  5793  1.601  1.814 


18.00 


BOUNDARY  LAYER  PARAMETERS  FOR  800  WPPM  POLYMER  INJECTED  AT  V,  -  .0383  FT/SEC 


TABLE  40 


18.00  2.822  371758  .236  .036  .022  13706  58589  8986  5506  1.632  1.795 


APPENDIX  E 

t 

VELOCITY  PROFILE  SAMPLES 


321 


322 


I 

I 


NOTE:  +10  -  add  .010  Inches  to  y  value;  +15  -  add  .015  inches  to  y  value 


VELOCITY  PROFILES  FOR  200  WPPM  POLYMER  INJECTED  AT  V.  -  .0383  FT/SEC 
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+10  -  Add  .010  Inches  to  y  value;  -2  -  subtract  .002  Inches  from  y  value 


NOTE:  +10  -  add  .010  Inches  to  y  value 
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tSCHE0,C0R£«L6,SCR»L,  CLASS-E 

sftnu(l.x.s) _ _ _ _ 

PROGRAM  SKINFRIC 

_ REAL.  X . L • LAMBOA  « INTGN01 « INTGN02 . INT  GRL 1 . INT GRL2 _ 

DIMENSION  LAMBDA (2561,* (256) ,CMX (256) , YINTG1 ( 256 ) , TINTG2 (256* 

_ BJ"ENSION_ri(Sai  .T2(S0) .OREX  (50) ,OIFINT(50>  .CF(256>  ,0ELTA(256) 

OIMENSI ON  REX (256) *  REXL (256) , WO  (256) , OTAO(256l , 0ELTA8L (256 ) 

_ OIMENSION  CFT  (256) _ _ _ 

COMMON  K, B, VOS, ALPHA, GAMMA, U, RE, CX,XCT , OELT , CXL 
_ l_FORMATi«F5.  0. 3F1 0. 0_.  £10. Q) _ 

2  FORMAT  ( 110,  3F10  •  0, 6F5,*Q  ) 

3  FORMAT (1H1,9X.  >CI»*  .F7.2,/,  10X,  > XCT-' , F6. 3, / ,  10X,  *V0S»%F5.  3, 

A  lOX.  'ALPHA-SF5.  3i/.'iOX,  'GAMMA*'  ,F5.  3./.10X,  *8«*,F5.3, 

B _ IPX, «K«‘ .F5.3,/, 10X, ‘HU-* .£12. 5./. IPX. «TEMP«« ,F6.1,/, 

C  lOX.'REXT-'.EIZ.S./.lOX, *XTRANS-‘ ,F6.3) 

_  L  FORMAT  (•  NO  ROOT  IN  GIVEN  RANGE'*  _ _ _ _ 

5  FORMAT)  9X~,  •X'.tX,  *REX>  ,OX,  'REXL'.OX,  'Dm  •  ,5X  ,  'DELf  A  ' .  AX  ,  'C?  * ", 

X  LX, ■LAM90A*,LX, <V S/VOS* ,6X, 'UTAU'.SX, •oeltabl’.tx, 'CF 

6  FORMAT(F10.2.ZF12.2.Fia.2.F10.3.Fl0.5,Fia.3,2'F10.L,F12.6.Fia. 

7  FOR MAT(1H1,18X, *X« ,11X, « INT1-INT2 * , 1SX . 1  ORE X  * ) _ 

a  Wma'I ('4x,*(«,^S.i, » , *,F6.2.  •) ’TSSToTT) 

C  INPUTS  " 

RTonTTnu?  “ 

REAP  l.U, VOS, L.‘CI,ALPHA, GAMMA, B.X,T,XCT, CXL, REXT _ 

£f(u,eo«o.oigo  to  99  11 

READ  2.NITER.EPS.RES,E0REX,XB,XF,XINC.STEPX,XCTL.XCTH _ 

_ .INITIALIZATIONS _ 

V»3.60901355c-5*T«(..7.17973611S-T»T«(7.15a9  0612E-9* _ 

I  T*(-3.$047?166£-11*T»6.6237797SE-1L) 1 ) 

N STEPS- ( X  F- X  B) /STEP**!.  1 _ 

UV»U/(V-1Z.  I 

XTR  ANS-  RE  XT/UV _ 

UV5-U/VOS 

NPT-XINC/STEPX+1.1 

- sjwsjraMj - - 

AL- ALPHA 

TFI NPT, EorZ*(NPT  72V) NPT-NPTaI 
OELX-  (NPT-O-STEPX 
MINT- (XF-XB) /OELX*.  OOOOV 
REX0-UV-XB 

p  RmrioTS -  ■  — 

OQ  15  I-t.NSTEPS 
*(I)«XB*(I-1)-STEP* 

5NX(I)«CNJCI,X(I)jlL) _ _ 

15  '  CONTINUE 
10  CONTINUE 

5’iiiMA-o.a 

ALPHA— 1.0  _  _ 

JS-1 

9IGX-XB  __  _ _ _ _ 

INTCRL1-0.0 

_ INTGRL2-0 . 0 _ 

_  GENERATE*  TABLE  OF  X.CN, LAMBOA  ANO  INTEGRANDS _ 
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_  00  30  1*1 .NSTEPS _ 

CX*CMXtZ> 

rex<i»«uv*xci> 

IF<X  I  I)'.L£.XCT)  GO  TO  27 
R£XL<H*UV»  (X(U-XCT* 

?r«s??r(T5  ' 

GAMHAaGH 

ALPHAaAL 

_ XL* 5 . _ 

F1*F(XL» 

oo  zs  j*6,4oa 

xk»  J 
F2*F«XHI 

tftifi  <f27Ct  .‘aTa  >w  ra~?s - 

F1*F2 

20  CONTINUE 

CO  TO  99 
2T  CONTINUE 

CALL  FALPOSC XL,  XM. NITER. RES. EPS. LANBOACH  * 

- SO  T(T  2*9 - - 

27  CONTINUE 

- RtxnTTSRBrm - 

L  ANBOA  <  I  >  *1 . 735525336*REX < I) *• . 25 
2B  CONTINUE  ~ 

RE*REXL(II 

y  INTC2  <  r«  -IMTSN02  JLANBOAT ITT  STOTT - 

7 1 NT Cl ( II *INTGN01 (LANBOA(II.X(I)I 

»vooj-uvs7c*Hwrm 

CF»I)  *2./  (LAMBDA  O1»LANB0A<Il  » 

DELTA (II*OELT 
UTAU  «I1 aVNOAI 1 *VQS 

- OELTA8LTTT*U*OECT*12i  /UTAIT(T) - 

30  CONTINUE _ 

C  PREFORN  INTEGRATION 


00  60  1*1, MINT 

- TIGX*BTCX*OELX - 

OR£XCI)*BICX*UV-REXO 

- TETTS*HPT-1 - 

N*0 

- 00  55  JPJS.JE - 

N*N-*1 

- - -  —  ■  ruHr*rlNTGitji - 

T2(NI*rINTS2CJI 

- 55 - CONTINUE - 

JS*Jc 

- “CALL  SIMPSON ( NPT,-5rEPX.n,T2.TIHirriN2T 

INTCRL1*INTCRL1*T1N1 

- TNTCRC2*TNTSRLT*YTN2 - 

0IFINT«II*tNTGRLl-INTSAL2 

- 50 - CONTINUE - - — - - 

RELER*<OIFINT  ININTI  »PREOIF)  /OIFINT  I MINT) 

- IFf  ABSTRELERSTLTIEOREX)  CO  TO  75 - 

0EL*01F1NT ININTT-OREX  TNINTT 

- BrRCOTF*  OIFINTT  KIND - 

IF(OEL.GT.a.O)GO  TO  65 

- ,CTM,XCT  "  - 

XCT*.5*<XCTL*XCTM) 

- oo  to  in - 


65  CONTINUE 
xCTU^Tg? 

XCT*.S*IXCTL*XCTH> 

_____ 

75  CONTINUE 

PRINT  3 ,CI,XCT , VOS, ALPHA , GANNA , B . K . V, T , REXT , XTR ANS 
PRINT  5 

00  79  1*1 ,NSTEPS 

IFIXIU.CT. XCT)  SO  TO  76 _ 

CFTII)»CFII) 

SO  TO  78 

76  iONTINUE 

IF IX ( I) ,GT. XTR ANS)  SO  TO  77 
GN*  IXIIY-xCTV/IXTRANS-XCT) 

_ XLL»1,73SS25336»R£XIH**«25 _ 

CFL*2./IXLL*XLL) 

_ CFTII)*  ll.»GH)«CFL*6H*CFII) _ 

SQ  TQ  TS 

77  CONTINUE 
CFTli)*CFUI 

_78 _ CONTINUE _ 

PRINT"  6,  XII)  ,R£XII>  .REXL(I)  ,CNX  I II  ,  OELTA  (I)  ,CFI  I)  ,LAH80A  (U  , 

X  VVOII) ,UTAUII) ,0ELTA8LII),CFTII) 

"79  CONTINUE 
PRINT  7 
8IGX*XB 

00  80  I*1,NINT 
8IGX*8IGX*0ELX 

PRINT  6.XB.8IGX.0IFINTII)  ,0REXII) _ 

__  ;QNTtNU£ 

GO  TO  9 
99  CONTINUE 
ENO 

SU8R0UfINg_FK>0SIXI..XH.N.R'.1frx3l 

SUBROUTINE  F6LP0S  USES  THE  nETHOQ  OF  FALSE  POSITION  TO  FI  NO  THE  »00T  OF 
A  NONLINEAR  EQUATION  .  THE  USER  SUPPLIES' f'HE  NONLINEAR  EQUATION  IN  THE 
FORM  OF  A  FUNCTION  ROUTINE  FIX)*0.0 
XL  *  LOW"  ESTIMATE  OF  ROOT  ~ 

XH  *  HIGH  ESTIMATE  OF  ROOT 
N  «  NUMBER  OF  ITERATIONS 
R  »  RE5I0UAL.ABSIFIANS))  ,LT.  Rl 

E  *  TOLERANCE  BETWEEN  SUCCESSIVE  APPROXIMATIONS  .  ABSIX1-»J)  .LT . E 
BOTH  R  A  NO  i  NAVE  TO  9E  SATISFIED  FOR  CONVERGENCE 

XJ_*  ROOT~{r^THE "nORlINE SJTEXHJJrTTEN - 

1  FORMAT  I '  SOLUTION  OOES  NOT  CONVERGE  IN', 15,'  ITERATIONS') 

xflXL  - - - - 

X2*XH 

00  10  1*1, N  - -  - 

X21*X2-X1 

x3'»xT*F  f  xiTTl  FTx  i  r=rrx2TT>  xTl - 

7*F|X3) 

IFTISST2  rTCTT  5T5'<T"T  8"7 - 

IF  I ABS 1X21) ,LT.E)RETURN 

f  rk  i  t  t  .  t  i  fin  r  '  ■  '  1 
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DIMENSION  Yl«l>,  Y2<  1! _ _ 

rlNTl»Tl«ll+riCN»*i;*YHN-ll 
YINT7-Y2  (1)*Y2(N)*4.«Y2(N-1) 
NiN-1 


t 

l 

r 

i 


00  ID  1*2, M, 2 

ftNf  l*Yl"NTl*4.»Yl  (t J*2.»Yl(I*t) 
YINT2*YINT2*4.«Y2(I)*2.»Y2(I+1> 
i'0  CONTINUE  - 
YINTl*H*YINTl/3. 

YINT2*H»Y'lNT2/3.  _ 
return 

gjjjj 

FUNCTION  CW(CI.X.L) 

re»l~l 

XL*  X/L 

x*.  oi»(fi*.oooi 

GO  TO  (100* 20 0,990,400, 500, 990 ,990,600) ,K 
^  CONTINUE 

CN*CI*< .97137*XL*(1.4495*XL* (-6.403* 3 .96539* XL) I ) 

IF ( XL . GT . . OS ) CM* . 007*Cl 
SO  TO  900 
200  CONTINUE  ~ 

_ CM*CI* ( .96372*XL*( 1.560 1*XL* (-5*  045574 3* 17097*XL) ) I 

fF(XL.St. ,671CM*.007*CI 
SO  TO  900 
400'TONTINUE 

CW*CI*( .9723S*XL»«1.2537*XL*(-4.6649*2.3152»XL> > ) 
fF(XL.St..?fllCM«;oi4»d 
SO  TO  900 

- j5¥-i6lit;irag - 

CW-CI»( 1. 00124* XL»(-.177746*XL* (.535740-1. 6 046*XL> ) ) 
TF( XL.GT. *93fCM* ,006*CI  — 

SC  TO  900 

"art  CONTINUE  “  ~~ 

CM»CI*).99383*XL»tl.0221S*XL*t-5. 9056*3. 9775»XL) )l 

- rrrxL.ST;;a9icu5:82s*ci - 

900  CONTINUE 

- iFTcM;sT.cr70R.x.LT.».>cu.ci - 

RETURN 

- 99<T  CONTINUE - 

STOP 

&NU 

FUNCTION  F (LAMBOA) 

- COMMON  KVHTVOSTTUPMITGirHHa .TTTREX , CT , X C  T , DEL  7  A , CXL - 

REAL  LANSOA.X 

- ULV5ir/L'Ah90A7V(JS - 

IF  (CX.GT.  CXL)  GO  TO  20 

- XU*7*  I L  PW  I*CX*  *  G  A  MH 1 - 

GO  TO  30 

20  CONTINUE - 

XX*K*  »L  PHA-CXL** GANNA 

30  CONTINUE  ”  —  "  "  ■” 

f*lAM80A*LAN80A— 4*  395604396* ALOG( • Q6*REX*UL V**XK)**2 

- ,ETqBH - _______ - - 

ENO 

- ■FUNcnax'TNTCNorn;  inbda  ;x  i - 

_ REAL  INTGNO 1  .LAN80A  ,  K 

YY*U/LAN80A/V0S 

- IC7STCSCT771 - 

_ IF  (CX.GT,  CXL)  SO  TO  20 
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za 

GO  TO  JO 

CONTINUE 

ZaALPHA  »C  XL  **GA  NNA 

30 

CONTINUE 

IPIX.GT.XCTJGO  TO  10 

OELTAaS.»X/SORT(REXI 
.  SO  TO  IS 

10 

CONTINUE 

0El,TAa£xP<Ka<LAH8OA-B>»»VV”»-<»Z) 

IS 

CONTINUE 

AOELaALOGIOELTAI 

»0L1»4DE1-*T;  '  —  ' 

PZ>Z./K*AOI.l*(B-l./On./IK*KI*AOEi.*M)eL»S*B 
_  INTGNOlaOELTA*  (Pl*fZ> 

RETURN  ----- 
END 

FUNCTION  INTGNOZ (LAMBOA.XI  ' 

_ REAL  INTGNOZ. LAHBOA.K 

con  NON  K . 8 , VOS . ALPHA , GA  H*A , U.RE  X ,CX .  XCTVOETT 77C  XL 
VVau/LAH80A/V0S 
ALVaALOGIVVl' 

IP  tCX.GT.  CXL>  GO  TO  ZO _ 

ZaALPHA»CX»»4AHHA 
GO  TO  30 

za  continue 

Z«ALPHA»CXL»*GAHHA 
JO  CONTINUE 

IPIX.GT.XCTIGO  TO  10 _ 

OEtTAaS • *X/SQRT (REX) 

GO  TO  IS 
10  CONTINUE 

Cr,Tl»CXP«K-;LANBOA-BM  *wa««-K»Z» 

IS  CONTINUE 

A0EL»AL06>0ELTA» 

AOLlaAOEL-l. 

Pl*OELT  A*LAHBOA*GANNAaZ/CX*ALV 
FZ«  1.  ♦  1 .  /  «•<>•  AOL  1*8-1./X 
INTGNOZaPt»PZ 
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